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The  research  work  in  this  report  was  performed  by  Chance  '/ought 
Corporation,  Dallas,  Texas,  for  the  Vehicle  Dynamics  Division,  AF  Flight 
Dynamics  Laooratory,  Wright-Patterson  Air  Force  Ease,  Ohio,  under  Contract 
AF33(6l6)-7£24.  This  research  is  part  of  a  continuing  effort  to  provide 
a  more  rational  and  reliable  theoretical  method  for  establishing  design 
criteria  in  the  specific  area  of  shock  and  impact  for  flight  vehicles  orb 
is  part  of  the  Research  and  Technology  Division,  Air  Force  Systems  Coimar.d''s 
exploratory  development  program.  The  Department  of  Defense  Program  Element 
number  is  6.24.05.33.4,  "Aircraft  Flight  Dynamics".  This  work  was  performed 
under  Project  ho.  1367,  "Structural  Design  Criteria"  and  Task  i!o.  136706, 
"Prediction  and  Measurement  of  Structural  Dynamic  Loads  Including  Fatigue 
Aspects".  Mr.  V/.  P.  Dunn  of  the  AF  Flight  Dynamics  Laboratory'  was  the 
Project  Engineer.  The  research  was  conducted  from  2  February  1961  to 
30  June  1962  by  the  Structures  Section  of  the  Aero  Division  of  Chance 
Vought  Corporation. 

The  report  is  Volume  I  of  a  two  volume  report  whicn  presents  the 
formulation  of  the  equations  of  motion  and  defines  in  detail  many  applied 
forces  of  present  and  future  vehicles  necessary  for  solution  of  the  equa¬ 
tions  of  motion  formulation,  Volume  II  -  Examples,  presents  illustrative 
examples  with  accompanying  numerical  solutions  so  that  the  procedure  for 
manipulating  the  equations  of  motion  formulated  in  Volume  I  can  be  used 
as  a  guide  when  the  method  is  utilized. 

Volume  I  contains  Sections  1  through  3*  Appendixes  A  through  D,  and 
the  Bibliography.  Volume  II  contains  Sections  4  and  5. 
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ABSTRACT 


A  rational  method  for  predicting  alighting  gear  loads  during  landing 
impact  is  discussed.  The  equations  describing  the  motions  of  a  vehicle  dur¬ 
ing  landing  impact  are  developed  for  an  arbitrary  vehicle  configuration. 

The  method  is  of  sufficient  generality  and  accuracy  to  allow  the  formulation 
of  alighting  gear  dynamics  problems  in  flight  vehicles  including  v/STOL  air¬ 
craft,  high  gross  weight  logistic  vehicles,  recoverable  booster  vehicles, 
advanced  tactical  and  defense  vehicles  operating  out  of  remote  areas,  anu 
lunar  vehicles.  It  allows  for  the  effects  of  varying  coefficients  of  friction 
and  damping,  combinations  of  initial  conditions  of  pitch,  roll,  yaw  angles 
and  rates,  vertical,  longitudinal  and  lateral  motion,  slippage  of  the  gear 
relative  to  the  alighting  surface,  flexible  alighting  gear  and  vehicle  struc¬ 
ture,  simultaneously  applied  triaxial  ground  loads,  and  various  types  and 
number  of  alighting  elements.  The  formulation  is  intended  for  the  cantilev¬ 
ered  type  of  gear,  although  the  articulated  type  may  be  handled  through  some 
extensions  of  the  formulation.  A  survey  of  the  various  types  of  forces  which 
occur  during  landing  impact  is  made,  and  the  manner  in  which  these  forces  en¬ 
ter  the  equations  of  motion  is  described. 

The  general  equations  nay  be  reduced  for  a  particular  problem  by  imposing 
the  vehicle  configuration  and  any  simplifying  assumptions  directly  on  the 
equations.  Several  illustrative  examples  with  accompanying  numerical  solu¬ 
tions  are  provided  in  Volume  II,  " Examples". 

The  report  may  be  used  as  a  guide  in  the  formulation  of  a  landing  impact 
problem. 


PUBLICATION  REVIEW 


This  report  has  been  reviewed  and  is  approved. 
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The  unit  vector  along  the  axis  of  the  body  coordinate 
system. 

The  unit  vector  along  the  of  axis  of  a  component  coord¬ 
inate  system. 

The  unit  vector  along  the  axis  of  the  articulated 
gear  coordinate  system.  * 


Paragraph  3*2. 5*2  -  The  tire  cornering  coefficient. 

jl  l  The  spring  rates  in  the  support  members  of  contacting 

elements  a.  jJtr. 


fCa  The  piston  bottoming  spring  rate. 

The  torsional  spring  rate  of  the  cantilevered  gear 
“  about  the  gear  centerline. 


^81 )  ^00  The  iower  and  uP-per  Piston  bottoming  spring  rates. 

iv  iv  The  lower  piston  bottoming  spring  rates  of  the  main  and 

'BM  ^  nose  gears. 


[KnJ 


The  stiffness  matrix  defining  the  normal  modes  of  vibra¬ 
tion. 


Kr 

k5N 

^TH  ) 

^XM  j  ^XN 


The  spring  .-.-nzc  of  the  relief  valve  spring. 

Metering  function  for  the  return  stroke. 

The  tire  spring  rates  of  the  main  and  nose  gear  tires. 

The  main  and  nose  gear  fore  and  aft  spring  rates  at  the 
axles. 


K'xmi  j  The  main  gear  fore  and  aft  axle  spring  rate  polynomial 

coefficients. 


K  -x',  .Kx'a  kj,'j  The  cantilevered  gear  fore  and  aft  axle  spring  rat' 

1  '  polynomial  coefficients. 

K  u'j  .  k*  iy'i  ,  Ku'y  The  cantilevered  gear  lateral  axle  spring  rate  poly- 
*  '  nomial  coefficients. 

|{xy  The-  cantilevered  gear  axle  spring  rater  are  con- 

’  *  1  sidered  as  functions  of  stroke. 


xiii 


ASD- TDR-62  555 
Volume  I 


mm 


a. 

L 

l 

u\ 

U 

i 


4 

4 

h 

4 

4  )  *-N 

4 


*Lr 

i  j 

Aun )  ^i/n 

L, 

Lj. 

CL.  ] 


fhe  stiffness  matrices  defining  the  potential  energy  due 
to  elastic  deformations  of  a  component  along  the  x', 
and  ^  axes  respectively. 

The  stiffness  matrices  arising  in  the  type  of  structure 
in  which  an  applied  force  in  one  direction  causes  a  de¬ 
flection  at  right  angles  to  that  direction.  These  ma¬ 
trices  are  not  retained  in  the  equations  of  motion 
except  in  the  appendix,  as  they  do  r.ct  occur  for  normal 
structures. 

The  position  vector  from  the  origin  of  the  body  coordinate 
system  to  the  undeflected  position  of  an  elemental  mass. 

The  Lagrangian;  kinetic  energy  minus  potential  energy. 

The  separation  of  the  spike  point  from  the  spike  center 
of  mass. 

The  distances  to  a  set  of  panel  points  from  the  u/  axis 
about  which  a  component  may  rotate  as  a  rigid  body. 

The  position  vector  of  the  axle  in  the  undeflected  vehicle. 

The  distance  to  the  center  of -mass  of  a  component  which 
may  rotate  as  a  rigid  body  about  a  line,  measured  from 
the  line. 

The  axle  offset  length. 

The  bearing  separation  at  zero  stroke. 

The  distance  from  the  cantilevered  gear  centerline  to 
the  support  strut  connection  point. 

The  length  of  the  j-th  bay  in  a  simple  beam. 

The  separations  of  the  upper  bearing  and  the  axle  of  the 
main  and  nose  gears. 

The  separation  of  the  upper  bearing  and  the  axle  of  the 
cantilevered  gear. 

A  vector  collinear  with  the  thrust  vector. 

Initial  bearing  separation  in  the  main  and  nose  gears. 

The  distance  from  the  uppor  bearing  to  the  contact  pal . 

The  distance  from  the  lower  bearing  to  the  contact  pad 
before  stroking  occurs. 

The  transformation  from  panel  px>int  loads  and  reaction 
forces  to  panel  point  shears  and  moments. 
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Matrices  formed  by  partitioning  CL,] 

ffhe  transf 01  .nation  from  panel  point  loads  to  suppe.-.-t 
reactions. 

The  transformation  from  panel  point  loads  to  panel  point 
shears  and  moments. 

The  vehicle  mass. 

The  mass  of  a  component  of  the  vehicle  which  may  move  as 
a  rigid  body  relative  to  the  body  ares  C& -  motion). 

The  mass  of  thei-th  component. 

The  bending  moment  in  the  J-th  bay. 

The  masses  of  the  main  and  nose  gears. 

The  mass  of  the  wheel,  tire,  brakes,  and  piston  in  the 
cantilevered  gear. 

The  mass  of  the  wheel,  t-'re,  and  brakes  in  the  canti¬ 
levered  gear. 

The  unsprung  mass. 

The  masses  supported  by  weightless  flexible  beams. 

The  exponent  which  indicates  the  exact  polytropic  nature 
of  the  pneumatic  compression. 

Spinup  moments  on  the  main  and  nose  gear  tires. 

The  sum  of  the  applied  moments  and  restraining  (scissors) 
spring  moment  c a  the  torsional  motion  of  the  cantilevered 
gear  unsprung  mass. 

The  spinup  moment. 

The  moment  on  the  cantilevered  gear  axle  due  to  the  force 
C?^'g>  which  acts  normal  to  the  wheel  plane. 

The  components  of  the  total  applied  moment  on  the  vehi "If . 

The  components  of  the  applied  moment  on  the  vehicle  due 
to  the  ground  force  on  a  contacting  element. 


The  moment  on  the  spike  about  its  center  of  mass  due  to 
h?6  ground  forces. 


jNp  )N^ 


Components  of  the  total  applied  moment  on  the  vehicle. 

These  are  non-orthogonal  components,  as  they  act  about 
t.he  2  axis,  the  line  of  modes,  and  the  x  axis  respectively. 
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The  bogie,  ski,  or  skid  restraining  moments. 

The  total  displacement  of  an  elemental  mess  ii.  the 
vehicle. 
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The  symbolic  form  for  the  panel  point  displacements. 

The  axle  total  displacement. 

Pressure  in  the  airchamber  of  the  pneumatic  spring. 

The  undeflected  uire  pressure. 

The  critical  pressure  in  a  gas  chamber  with  relief  valve. 
The  tire  wall  equivalent  pressure. 

The  average  tire  contacting  pressure. 

Tire  pressure  rise  on  deflection. 

The  total  displacement  of  an  elemental  mass  in  the  i-th 
component  of  the  vehicle. 

The  displacement  due  to  elasticity  of  an  elemental  mass 
in  the  i-th  component  of  the  vehicle. 

The  displacement  of  an  elemental  mass  in  the  l-th  com¬ 
ponent  due  to  displacement  of  the  j-th  component,  to 
which  tne  i-th  component  is  affixed. 

The  displacement  of  the  lower  bearing  in  the  cantilevered 
gear,  at  right  angles  to  the  gear  centerline. 

The  components  of  PL 

The  displacement  of  the  upper  bearing  in  the  cantilevered 
gear,  at  right  angles  to  the  gear  centerline. 

The  components  of  Py. 

The  vector  displacement  of  the  instantaneous  center  of 
maos  of  the  vehicle  from  the  origin  of  the  body  coordinate 
system. 

The  components  of  P,'  . 

The  panel  point  displacements  parallel  to  the  component 
coordinate  axes. 


a 


The  panel  point  displacements  due  to  elasti:ity. 

The  panel  point  displacements  in  the  i-th  component  due 
to  displacement  of  the  j-th  component,  to  which  the  i-th 
component  is  affixed. 
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Components  of  the  displacement  of  panel  point  i.  on  the 
piston  due  to  displacement  of  the  cylinder. 


P*n,P: 


The  axle  fore  and  aft  displacements  of  the  main  and  nose 
gears. 


P.»P* 

0 

w 


Pressures  in  the  two  regions  of  a  liquid  spring. 

The  total  applied  force  on  the  vehicle. 

A  column  vector  of  modal  coordinates. 

A  column  vector  of  modal  coordinates  with  the  higher 
modes  deleted. 
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QP 
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The  aerodynamic  force  on  the  vehicle . 

The  generalized  force  associated  with  the  generalized 
coordinate 

A  generalized  coordinate. 

The  parachute  fores  on  the  vehicle. 

The  thrust  vector  and  its  magnitude. 

The  gravitational  force  on  the  vehicle. 

The  components  of  the  total  vehicle  force  along  the  axes 
of  the  ground  coordinate  system. 
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$c  $8 

Qtj  &&■ 
Qxgn  ,0i 


Q^'gp 


The  components  of  the  total  vehicle  force  along  the  axes 
of  the  body  coordinate  system. 

The  panel  point  applied  force  components. 

The  panel  point  load  components;  the  difference  between 
the  applied  forces  and  the  inertial  reactions. 

The  components  of  the  ground  force  on  the  vehicle. 

Two  components  of  the  ground  force  on  the  contacting 
elements  a,  Jr, 

Two  components  of  the  ground  force  on  the  main  and  nose 
gear  tires. 

Components  of  the  total  applied  force  on  the  piston. 

The  ground  force  acting  along  the  piston  ^  -axis. 

The  aerodynamic  lift. 
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The  total  stroking  force;  the  sum  of  the  hydraulic, 
pneumatic,  bearing  friction,  and  bottoming  forces. 

The  total  stroking  forces  in  the  main  and  nose  gears. 

The  components  of  the  gravitational  force. 

The  position  vector  of  the  origin  of  the  body  coordinate 
system. 

The  position  vector  of  an  elemental  mass  in  the  vehicle 
relative  to  the  ground  coordinate  system. 

The  transformation  from  the  Eulerian  angle  time  deriv¬ 
atives  to  the  components  of  the  angular  velocity  expressed 
in  the  body  coordinate  system. 

Tile  tire  undeflected  radius. 

On  the  cylinder;  components  of  the  support  reaction  along 
and  at  right  angles  to  the  gear  centerline.  On  the  piston, 
the  bearing  forces. 

The  effective  drag  arta  of  the  parachute.  The  symbols 
may  be  separately  defined  as  the  drag  coefficient  and 
area,  but  no  standard  definition  holds  (see  Ref.  5)* 

The  effective  area  is  generally  measured  for  each  para¬ 
chute. 

The  stroke  of  the  piston;  the  displacement  of  the  piston 
relative  to  the  cylinder,  positive  as  the  piston  enters 
the  cylinder. 

The  component  of  the  ground  force  on  a  contacting  element 
parallel  to  the  ground  Y  axis. 

The  kinetic  energy  of  the  vehicle. 

Time 

The  matrix  relating  the  displacements  at  the  bearing 
points  to  the  displacements  of  the  panel  points  on  the 
cylinder.  The  elements  of  this  matrix  are  determined  by 
the  stroke,  the  position  of  the  panel  and  the  inter-:  - 
lation  scheme  used  to  relate  the  cylinder  elastic  dis¬ 
placement  to  the  panel  point  displacements. 

Tlie  matrix  yielding  the  panel  point  displacements  of  the 
i-th  component  as  a  rigid  body  due  to  displacement  of  tlie 
j-th  component,  to  which  +!’.e  i-th  component  is  affixed. 

The  matrix  yielding  the  panel  point  displacements  of  the 
j-th  component  as  a  rigid  body  due  to  displacement  of  the 
fuselage,  to  which  the  j-th  component  is  affixed. 
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The  general  geometrical  transformation  between  shears 
and  moments,  and  the  panel  point  applied  forces  and  re¬ 
actions. 


The  matrix  yielding  the  panel  point  displacements  of  the 
piston  as  a  rigid  body  due  to  displacements  of  the  cylin¬ 
der  (or  the  bearing  points). 

The  spinup  tine. 

The  matrix  yielding  the  panel  point  displacements  o:.  the 
tail  as  a  rigid  body  due  to  fuselage  panel  point  dis- 
pla cements. 


[TvH3 


The  matrix  yielding  the  panel  point  displacements  of  the 
vertical  beaa  as  a  rigid  body  due  to  panel  point  dis¬ 
placements  of  the  horizontal  beam. 


[TWf] 


The  matrix  yielding  the  panel  point  displacements  of  the 
wing  as  a  rigid  body  due  to  panel  point  displacements  cf 
the  fuselage. 


u 

The  potential  energy  d  e  to 

elastic  deformations. 

V 

Volume . 

r 

The  velocity  of  propagation 

cf  a  compression  wave. 

ij- 

The  velocity  of  the  vehicle 

relative  to  the  atmosphere. 

IV 

The  velocity  of  the  metered 
Bernoulli's  relation. 

hydraulic  cil  obtained  free. 

[nr] 

Defined  by  2g.  2.7-12. 

V8  The  velocity  of  a  point  3  in  the  vehicle,  relative  tc 

Uie  ground. 


Vo 


The  component  of  the  ground  force  on  a  contacting  ele¬ 
ment  parallel  to  the  l  axis. 


)  V4N 
Vi 
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The  ground  force  \/s  oa  the  main  and  nose  gear  tires. 

The  par.ei  point  shear  load  at  the  j-th  bay. 

Tr.e  components  of  the  ski  or  skid  axle  velocity  at  ng.tt 
angles  tc  tne  ej.eo.cnt. 


VM 

V, 


Hfp 


The  volume  of  oil  metered  tnreugr.  the  orifice. 

The  initial  '/slue  of  the  volume. 

The  velocity  of  the  parachute  attachment  point 
tc  the  atmosphere. 


XIX 


The  speed  of  the  tire  footprint  relative  to  the  ground. 


The  components  in  the  body  coordinate  systen  of  the 
vehicle  velocity. 


V  \l  \t  The  components  of  the  velocity  of  a  point  B  on  the  body 
*8'  "5®  relative  to  the  ground. 


The  components  in  the  body  coordinate  system  of  the 
vind  velocity. 


The  initial  volumes  of  the  two  regions  of  the  liquid 
spring . 


3Vit  aV  3 \J  The  rates  at  which  the  volume  V,  of  the  liquid  spring 

— ^  •>  mr-pf  cylinder  changes  with  pressure  and  stroke  due  to  cylin- 

a  |  3A  o  r,  jjer  expansion  and  seal  compression. 


UT  The  tire  width. 


WD  W  WN  The  components  of  the  ground  force  cn  a  ski  or  skid  in 
1  1  the  coordinate  system  of  the  element. 

VvL  The  component  of  the  ground  force  on  a  tire,  in  the 

"  ground  plane  and  normal  to  the  line  of  intersection  of 

the  wheel  plane  and  the  ground  plane. 
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The  spinup  force;  the  component  of  the  ground  force  on 
a  tire,  in  the  ground  plane  and  parallel  to  the  line  of 
intersection  of  the  wheel  plane  and  ground  plane. 

A  component  of  the  position  vector  |R  in  the  ground 
coordinate  system. 

A  component  of  the  position  vector  IL  in  the  body  coord¬ 
inate  system. 

A  component  of  the  position  vector  U.  in  a  component  co¬ 
ordinate  systen. 

A  component  of  the  position  vector  IL  to  the  center  cf 
mass  of  a  vehicle  component,  expressed  in  the  component 
coordinate  system. 

The  panel  point  positions  along  the  x^&xis. 

The  component  of  the  vehicle  velocity  along  the  y  -axis. 


X,  The  component  of  the  position  vector  to  the  contacting 

element  'a'  along  the  x-axis. 

?  ,  3:..,  Tne  components  of  the  position  vector  U  to  the  nose  and 

fin  )  _  ,  , 

rati  1:1  gear  axxec  along  the  x-axis. 
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The  component  of  til c  ve lociay  of  a  point  3  in  the  vehicle, 
relative  to  the  ground  and  along  the  %~&xis  ■ 

The  cor.pcr.ent  of  the  velocity  of  the  trailing  end  of  a 
ski  or  skid,  relative  to  the  ground  and  along  the  X -axis . 

The  initial  value  of  the  component  of  fR  along  the  X-axis. 

The  component  of  the  pad  velocity  relative  to  the  ground 
along  the  X-axis. 

The  panel  point  coordinates  of  the  wing  panel  points. 

The  coordinate  of  the  ving  root. 

A  component  of  the  position  vector  fR  in  the  ground  co¬ 
ordinate  system. 

A  otcpor.er.t  of  the  position  vector  |Lir.  the  cody  coord¬ 
inate  oyster. 

A  oor.pcr.er.t  of  the  position  vector  LL  ir.  a  ccrpcr.er.t  co- 
srci^.^  «c  sys 

A  component  of  the  position  vector  iLto  the  center  of 
rass  of  a  vehicle  component,  expressed  ir.  the  component 
coordinate  system. 

The  pane1  point  positions  along  the  <^-axis. 

The  component  of  the  vehicle  velocity  along  the  Y-axis. 

The  component  of  the  velocity  of  a  point  3  ir.  the  vehicle, 
relative  to  the  ground  and  along  the  Y-axis. 

The  iST.pcre.vt  of  the  velocity  of  the  trailing  end  of  a 
:.<i  .r  skid,  relative  to  the  ground  and  along  the  Y-axis. 

T..=  t-.-ttsl  value  of  the  component  of  IR  along  the  Y-axis. 

The  component  of  the  pad  velocity  relative  to  the  ground 
along  the  Y-axis. 

The  distance  free  the  ground  to  the  origin  of  the  tody 
coordinate  system,  normal  to  the  ground. 

A  component  of  the  position  vector  !L  ir.  the  fcedy  coord- 

'..12*.^  3‘‘ 3  « 

A  : impotent  of  the  position  vector  IL  in  a  component  co¬ 
ordinate  system. 

A  component  of  the  position  vector  !Lto  the  center  of 
mats  of  a  vehicle  component,  expressed  in  the  component 
coordinate  system. 
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The  panel  point  positions  along  the^'-axis. 

The  component  of  the  vehicle  velocity  aortal  to  the 
ground  plane. 

The  axle  height  above  the  ground . 

The  rain  and  nose  gear  axle  coordinates  along  the  -j  -axis. 

The  height  of  point  B  above  t.-ie  ground. 

The  velocity  of  the  trailing  end  of  a  ski  or  skid,  rel¬ 
ative  to  the  ground  and  along  the  2 -axis. 

The  pad  velocity  relative  to  tr.e  ground  along  the  1  -axis. 

The  initial  value  of  the  component  of  IR  along  the  2 -axis. 

The  distance  along  the  spike  measured  from  the  point. 

Paragraph  j.2.-»  -  The  vehicle  angle  of  attack. 

Paragraph  -.3. 5*1  -  The  angle  at  vnic.t  the  cantilevered 
gear  shock  strut  is  canted  forward  from  the  ^-axis. 

Paragraph  -,2o.2  -  The  angle  of  rotation  of  the  artic- 

T he  angle  betveer.  tr.e  -T.ee  1  axle  and  the  ground  plane. 
Paragraph  3*2.-  -  The  vehicle  sideslip  angle. 


raragrap.c  >r. 


-  The  angle  through  vhich  the  vheel 

axle  is  rotated  accut  tr.e  gear  centerline  due  to  tor- 

-  .  ... 
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P 

[n] 

Ljp  ] 

^A 


.e  ang_e  t.trougn  vhicr.  a 
i  rear  centerline  due  to 


ski  or 
torsional 


Tr.e  spt.-ce  apex  angr 


The  transformation  from  the  ground  coordinate  syst.r. 
to  tne  scdy  ocorcinate  systec.. 


The  transformation  from  a  component  coordinate  system 
to  tr.e  body  coordinate  system. 

The  angle  cer-eer.  the  aiiitulated  gear  tension  strut 
centerline  and  the  vheel  axle. 


The  angle  betveen  the  support  strut  and  the 
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y'jr .  The  cosines  of  the  angles  between  the  thrust  vector  and 

°  '  '  the  axes  of  the  body  coordinate  system. 

£  Paragraph  3*2.**  -  A  control  surface  deflection. 

^  The  tire  deflection. 

£  A  roll  control  surface  deflection,  such  as  aileron  or 

A  spoiler. 

^  i,  The  tire  deflection  at  bottoming. 

^  A  pitch  control  surface  deflection. 

^  A  flap  deflection. 

The  displacement  of  an  elemental  mass  from  its  undeflected 
1  position  defining  component  rigid  body  displacement  with 

respect  to  the  body  axes,  expressed  in  the  component  co¬ 

ordinate  system  in  which  the  motion  is  most  easily  de¬ 
scribed. 


jAyi 
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A  ruduer  deflection. 

The  main  and  nose  gear  tire  deflections. 

The  panel  point  displacement  components  due  to  ^-motion 
of  the  component. 

Tte  rigid  body  displacement  along  the  stroking  axis 
(A  -motion  along  a  line). 

The  displacements  alo:ig  the  stroking  axes  of  the  con- 
tactiiig  elements  a,  b. 

The  displacements  along  the  stroking  axes  of  the  main 
and  nose  gear  pistons. 

The  fraction  of  critical  damping  of  the  mode  of  fre¬ 
quency  co. 

The  angle  defining  component  rigid  body  motion  about  a 
line. 

The  angular  velocity  of  a  component  rotating  as  a  rigid 
body  about  a  line. 

Wheel  spinup  angular  velocity. 

The  main  and  nose  gear  wheel  angular  velocities  at  spin- 
up. 

One  of  the  Euler  angles  defined  in  Paragraph  2.7. 
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Qt  The  initial  value  of  8. 

2  The  angle  between  the  component  of  axle  velocity  parallel 

to  the  ground  and  the  line  of  intersection  of  the  wheel 
plane  and  the  ground  plane. 

2  _  An  eigenvalue  of  the  fuselage  homogeneous  panel  point 

F  equations. 

2  ■  An  eigenvalue  corresponding  to  the  i-th  mode  of  elastic 

1  vibration. 


M- 

Sb 


The  coefficient  of  friction  between  a  contacting  element 
and  the  ground-. 

A  bearing  coefficient  of  friction. 

The  lower  and  upper  bearing  coefficients  of  friction. 

The  lower  and  upper  liquid  spring  seal  coefficients  of 
friction. 


A  dimensionless  local  bay  coordinate  used  to  non-dimen- 
sionalize  some  forms  i.i  the  interpolation  schemes. 


P 
P  H 

\r\ 

* 

\<p\ 

[0] 

L£] 

<P, 

L> 

Si 


The  vehicle  mass  density;  or,  the  atmospheric  density. 

The  hydraulic  fluid  density. 

A  column  ratrix  of  undetermined  multipliers. 

The  undetermined  multiplier  associated  with  the  constraint 
relation  F^.. 

One  of  the  Euler  angles  defined  in  Paragraph  2.7* 

A  column  matrix  defining  a  node  shape;  the  elements  are 
proportional  to  the  actual  panel  point  displacements  in 
the  mode. 

A  square  matrix  of  modal  columns  {$} 

A  rectangular  matrix  of  modal  columns  in  which  the  higher 
modes  are  deleted. 

The  initial  value  of  (f>- 

One  of  the  Euler  angles  defined  in  Paragraph  2.7- 

The  initial  value  of  "'ft 

The  frequency  of  vibration  of  a  mode. 

The  vehicle  angular  velocity. 
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ilx-  jSln'  jft-y 


defined  by  Bq.  2.7-H* 

The  components  of  JL  in  the  body  coordinate  system. 

The  components  of  Jl  in  a  component  coordinate  system. 


Subscripts 
A  Axle;  axis 

a  A  point  in  the  vehicle;  a  contacting  element 
B  A  point  in  the  vehicle 
b  A  contacting  element 
c  Cylinder 

e  Ski  or  ski  trailing  end 
F  Fuselage 
f  Final  value 
G  Ground 
H  Horizontal  beam 

i  One  of  ‘.he  vehicle  components,  panel  points,  modal  coordinates,  etc. 

This  subscript  must  be  determined  in  context. 

j  One  of  the  vehicle  components,  panel  points,  local  bays,  etc.  This  sub¬ 
script  must  be  determined  in  context. 

L  Lower  beariig  or  seal 

M  Main  gear 

K  ‘lose  gear 

sn.  A  contacting  element 

o  The  initial  value 

p  Pad,  or  piston 

S  Return  stroke 

T  Tire 

U  Upper  bearing  or  seal 
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V  Vertical  beam 
W  Gravitational,  wind,  wing 

H  ,  A  Components  along  the  x,  and  axes;  or,  in  a  matrix, 

"  referring  to  the  component  of  displacement  in  the  corres¬ 

ponding  kinetic  energy  term. 

x>  )<1t )  "V  Components  along  the  and  c!  axes;  or,  in  a  matrix 

"  6  referring  to  the  component  of  displacement  in  the  corres¬ 

ponding  kinetic  energy  term. 
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SECTION  1 
INTRODUCTION 


The  procedures  accepted  by  the  aircraft  industry  for  determining  the 
design  loads  for  alighting  gear  were  until  recent  years  based  on  the  assump¬ 
tion  that  an  adequate  design  could  be  obtained  from  consideration  of  a  few 
discrete  conditions.  These  conditions  and  the  procedures  for  determining 
the  associated  loads  were  established  in  the  specifications.  The  great  in¬ 
creases  in  airplane  performance  subsequent  to  World  War  II,  however,  cauaed 
the  reinvestigation  of  many  well  established  criteria.  The  introduction  of 
high  speed  computing  equipment  about  this  same  time  made  possible  the  solu¬ 
tion  of  problems  that  previously  were  too  long  end  complicated  to  permit 
economical  solution  on  a  large  scale. 

The  aircraft  alighting  gear  was  recognized  as  one  design  area  that  could 
benefit  from  the  performance  of  more  detailed  analyses  and  the  rationalization 
of  procedures  for  determining  design  loads.  The  first  efforts  in  this  direc¬ 
tion  led  to  the  consideration  of  mathematical  models  which  had  one  or  two 
degrees  of  rigid  body  freedom.  These  efforts  eventually  led  to  a  multiple 
degree  of  freedom  analysis  based  on  a  rigid  airplane  free  to  translate  hori¬ 
zontally  and  vertically  and  to  pitch.  The  motion  of  the  airplane  was  defined 
in  the  ground  reference  system.  The  equations  of  motion  were  simple,  although 
non-linear,  equations. 

Comparison  of  the  numerical  data  obtained  from  solving  these  equations 
with  the  data  obtained  from  airplane  flight  and  drop  test  programs  indicated 
that  the  equatiors  could  be  made  to  yield  more  accurate  answers  provided 
additional  physical  mechanical  details  were  considered  in  the  analytic  form¬ 
ulation.  It  eventually  became  apparent  that  considerable  simplicity  could 
be  gained  by  using  the  body  axis  system  as  the  primary  coordinate  system  in 
lieu  of  the  ground  reference  system.  Subsequently,  the  equations  of  motion 
were  written  in  the  body  reference  system  and  due  to  the  pitch  rotation  of 
the  airplane  a  Coriolis  force  term  appeared  on  the  left  tend  side  of  the 
equation. 

The  addition  of  other  degrees  of  freedom  for  notion  of  the  rigid  vehicle 
and  the  effects  of  structural  deflections  and  other  pertine.it  considerations 
caused  the  continued  build-up  of  the  equations  of  motion  and  the  concomitant 
fen  >  and  geometric  expressions.  By  this  process,  which  might  be  called  a 
building  block  procedure  (that  is,  a  procedure  which  is  built  up  term  for 
term  r.s  the  necessity  for  each  term  is  recognized),  the  procedure  for  d.. iorm- 
ir.ing  alighting  gear  loads  was  extended  to  a  multiple  degree  of  freedom  anal¬ 
ysis  which  included  not  only  the  six  rigid  body  degrees  of  freedom,  but  also 
contained  degrees  of  freedom  for  a  flexible  airframe.  The  buiiding  block 
approach  has  provided  a  procedure  that  is  very  adequate  for  the  determination 
of  alighting  gear  loads  for  conventional  airplanes. 
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The  advent  of  flight  vehicles  that  cannot  be  classed  as  conventional 
airplanes  has,  however,  created  new  problems.  It  is  no  longer  possible  in 
all  cases  to  define  for  an  unknown  configuration  the  parameters  important  to 
the  proper  solution.  Consequently,  at  the  beginning  of  this  program,  it  was 
decided  that  providing  a  rational  method  that  would  fulfill  the  stated  re¬ 
quirements  would  necessitate  a  departure  from  the  building  block  approach. 

A  development  program  was  therefore  established  on  the  basis  of  a  completely 
general  approach  to  the  problem. 

In  an  effort  to  provide  a  broad  base  of  understanding  of  the  program, 
and  of  this  report,  the  salient  ideas  and  concepts  involved  in  the  program 
and  the  general  format  of  this  report,  are  presented  in  the  following  para¬ 
graphs.  The  ideas  and  concepts  will  be  discussed  with  the  intention  of  pro¬ 
viding  the  reader  with  a  word  picture  of  the  processes  being  employed.  To 
provide  a  rational  method  for  predicting  alighting  gear  loads  that  would 
attain  the  wide  applicability  desired  has  necessitated  the  utilization  of 
numerous  mathematical  concepts  and  tools.  As  is  frequently  the  case  in  math¬ 
ematical  developments,  the  derivation  of  a  very  general  result  requires  great¬ 
er  effort  than  the  development  of  a  specific  result. 

Since  the  problem  being  solved  is  a  complex  problem  in  structural  dynamics, 
it  is  obvious  that  the  development  of  the  method  will  require  employment  of 
the  concepts  and  tools  of  dynamics.  However,  the  method  is  presented  in  a 
concise,  easy  to  follow,  didactic  manner  so  that  the  average  dynamics  engineer 
can  readily  apply  this  method  to  his  specific  problem.  Consequently,  the  der¬ 
ivation  of  the  equations  of  motion  are  presented  under  the  assumption  that  the 
careful  reviewer  of  the  mathematical  development  has  a  working  knowledge  of 
the  Lagrangian  equations  of  motion,  vector  calculus,  matrices,  and  structural 
analysis  methods.  The  presentation  of  the  procedures  involved  in  the  use  of 
the  method  will  not,  however,  require  detailed  knowledge  of  these  concepts, 
and  the'aim  of  the  presentation  will  be  to  facilitate  use  of  the  method  by 
the  engineer.  As  a  preview  to  tae  technical  presentation,  the  highlights  of 
the  development  of  the  method  are  presented  here  in  the  introduction. 

Several  iterations  were  necessary  to  the  evaluation  of  the  equations  of 
motion  in  their  final  form.  After  the  first  formulation  of  the  general  equa¬ 
tions,  it  became  evident  that  a  simplification  of  expression  was  necessary  to 
allow  the  equations  to  be  presented  in  a  more  concise  manner.  The  equations 
were  then  rewritten  in  matrix  form,  which  provided  a  means  of  collecting  terms 
in  a  manner  that  would  allow  the  equations  of  motion  to  be  viewed  as  a  set  of 
equations  rather  than  as  individual  equations.  As  work  on  the  project  pro¬ 
gressed  to  the  consideration  of  some  of  the  detailed  problems  concerned  with 
repr-s^-.tation  of  a  complex  structure  consisting  of  several  components,  it 
was  found  that  the  inertial  characteristics  could  be  handled  more  readily  if 
the  characteristics  of  each  component  were  expressed  in  its  own  coordinate 
system.  Consequently,  the  equations  governing  each  component  motion  relative 
to  the  vehicle  were  expressed  explicitly,  and  the  general  equations  of  motion 
were  again  rewritten.  This  form  provides  generality  and  also  is  readily  a- 
daptable  to  numerical  analysis. 

To  provide  a  mathematical  model  of  the  vehicle  that  is  of  sufficient 
generality,  the  position  of  every  point  in  the  vehicle  is  defined  :n  terms 
of  the  position  of  a  discrete  sot  of  points  in  the  vehicle,  called  panel 
points.  Mass  is  considered  to  be  distributed  throughout  the  vehicle; 


2 


distribute-?  mass  inertial  characteristics  are  considered,  but  concentrated 
mass  points  may  later  be  used  for  a  particular  model.  The  panel  points  are 
allowed  only  to  translate  with  respect  to  each  other  along  three  mutually 
perpendicular  axes.  The  mathematical  model  consists  of  N  components,  such 
as  wing,  tails,  landing  gears,  which  are  elastic  when  the  panel  points  on  a 
component  are  allowed  relative  motions,  or  which  may  be  made  rigid  by  allow¬ 
ing  no  relative  motions.  Rigid  body  displacement  of  a  set  of  panel  points 
with  respect  to  a  component  coordinate  system  is  also  allowed  in  order  to 
account  for  motions  such  as  those  of  control  surfaces  and  gear  stroking. 

With  the  basic  mathematical  model  thus  establisheu  it  is  possible  to 
proceed  with  the  establishment  of  a  workable  notation,  to  define  the  coord¬ 
inate  systems  to  be  employed,  and  to  establish  the  procedure  required  to 
maintain  knowledge  of  the  position  of  the  panel  points.  The  right  hand  rule 
is  used  in  defining  the  coordinate  systems.  The  ground  coordinate  system  is 
considered  to  be  an  inertial  frame  of  reference.  The  body  coordinate  system 
is  defined,  as  are  the  component  coordinate  systems,  and  the  transformations 
from  one  system  to  another  are  established;  hence,  the  wing  coordinates  of  a 
point  on  the  wing  can  be  transformed  to  body  coordinates  or  to  the  inertial 
frame  by  these  transformations.  Maintaining  knowledge  of  the  position  of  the 
panel  points  is  accomplished  by  defining  position  vectors  such  that  the  posi¬ 
tion  of  every  panel  point  relative  to  the  inertial  frame  is  known. 

To  provide  the  desired  generality  in  the  most  direct  manner,  the  Lagrang- 
ian  equations  of  motion  are  employed  as  the  foundation  of  the  development. 

The  use  of  the  Lagrangian  equations  require:,  that  the  potential  and  kinetic 
energy  of  the  system  be  defined  in  terms  of  the  ground  system  coordinates, 
but  in  general  this  is  at  best  a  very  difficult  task.  However,  if  use  is 
made  of  the  body  coordinate  system,  the  kinetic  and  potential  energy  of  the 
system  car.  be  readily  defined.  It  would  be  desirable,  then,  to  express  the 
lagrangian  equations  in  terms  of  body  coordinate  system  variables.  This  can 
be  done  by  writing  a  set  of  Lagrangian  equations  to  defire  the  translational 
and  rotational  motion  of  the  body  coordinate  system  and  another  set  to  define 
the  location  of  all  particles  in  the  body  with  respect  to  the  origin  of  the 
body  coordinates.  The  rotational  position  is  defined  in  terms  of  Eulerian 
angles.  These  angles  are  then  so  chosen  that  a' transformation  from  body 
coordinates  to  the  inertial  frame  is  always  possible.  It  is  therefore  poss¬ 
ible  to  transform  the  inertial  velocities  to  body  velocities. 

Since  the  kinetic  energy  of  the  body  is  written  m  terms  of  the  body 
velocities,  the  elastic  displacements  of  a  discrete  number  of  panel  points 
from  their  equilibrium  position,  and  the  velocities  of  these  panel  points  in 
their  components,  the  requirement  is  to  express  the  Lagrangian  equations  in 
the  terms  used  for  the  kinetic  energy.  The  development  of  these  equations 
was  performed  for  this  program.  The  results  of  transforming  the  Lagrangian 
equations  are  shown  as  "modified  Lagrangian  equations".  The  general  equations 
of  motion  are  then  derived  by  substituting  the  equations  for  the  kinetic  and 
potential  energy  into  the  "modified  lagrangian  equations"  and  performing  the 
indicated  operations.  The  resulting  equations  are  long  and  contain  numerous 
terms  that  by  inspection  can  be  seen  to  be  of  no  importance  to  the  alighting 
loads  problem.  In  the  general  presentation  in  Appendix  C  all  term-;  are  .shown, 
but  for  the  purpose  of  this  report,  the  unimportant  terras  have  beer  deleted 
in  Section  2,  where  the  resulting  equations  of  motion  governing  landing  im¬ 
pact  are  shown.  Equations  2.8-3  and  2.8-4  are  referred  to  as  "rigid  body" 
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equations  since  the  principal  terras  are  the  inertial  terms  of  the  undeflected 
vehicle;  the  remaining  terms  on  the  left  hand  side  are  inertial  coupling  terras 
due  to  elastic  deformations  and  component  motions.  The  terms  on  the  right 
hand  side  of  Equations  2.8-3  and  2.8-4  are  the  external  forces  and  moments, 
respectively,  written  in  the  body  coordinate  system.  Equation  2.8-5  without 
the  "rigid  body"  coupling  tens  are  the  ordinary  panel  point  equations  of  a 
restrained  body,  with  the  panel  point  forces  and  constraint  forces  appearing 
on  the  right  hand  side  of  the  equation.  These  three  equations  must  be  solved 
simultaneously  with  any  existing  constraint  equations  in  order  to  obtain  a 
solution. 

It  is  intended  that  the  formulation  presented  here  be  of  sufficient  gen¬ 
erality  to  fulfill  the  requirements  of  any  particular  problem  concerning 
landing  impact.  It  is  also  intended  that  the  manner  of  presentation  of  tht 
method  be  of  such  clarity  to  permit  utilization  of  the  method  by  the  practic¬ 
ing  engineer. 

The  requirement  for  generality  has  caused  some  complexity  in  the  formu¬ 
lation,  which  has  resulted  in  a  considerable  amount  of  discussion  to  explain 
fully  the  myriad  of  details  covered  by  the  program.  The  intent  of  this  lengthy 
discussion  is  to  allow  the  user  of  the  method  the  freedom  of  choosing  the 
particular  segments  applicable  to  his  problem.  By  proper  choice  of  terms 
from  the  general  expressions,  consistent  with  any  simplifying  assumptions  that 
have  been  made,  the  equations  of  motion  for  a  particular  problem  can  be  ob¬ 
tained  from  this  report.  The  utility  of  this  procedure  depends  upon  whether 
or  not  significant  time  can  be  saved  in  formulating  the  desired  equations  of 
motion  by  using  this  report,  as  contrasted  to  the  time  required  to  formulate  . 
the  equations  independently. 

To  assist  in  attaining  the  desired  utility,  it  has  been  necessary  to 
temper  somewhat  the  desire  for  generality  of  the  final  expressions  describing 
the  mathematical  model.  Therefore,  for  the  development  in  the  main  body  of 
the  report  certain  assumptions  have  been  made  which  eliminate  some  terms 
whose  effect  on  alighting  gear  loads  or  vehicle  motion  is  negligible  for  any 
known  or  postulated  vehicle  configuration. 

Section  2  of  this  report  is  devoted  to  the  development  of  the  equations 
of  motion.  This  section  presents  a  somewhat  detailed  expla-.iation  of  the 
various  mathematical  concepts  employed  in  the  development  of  the  method.  It 
has  not  been  possible  due  to  space  requirements  to  include  in  Section  2  all 
the  algebraic  manipulations  required  in  the  development.  As  noted  above, 
certain  simplifying  assumptions  made  in  Section  2  have  resulted  in  the  omission 
of  v^rtain  terms  in  the  final  equations.  For  the  benefit  of  those  dynamicists 
who  are  interested  in  the  complete  development  of  the  equations,  Appendices 
A,  B,  and  C  are  included.  In  these  appendices  the  entire  development  -'.s  pre¬ 
sented  without  recourse  to  simplification. 

For  those  engineers  not  interested  in  the  development,  the  final 
equations  of  motion.  Equation  2.8-3,  2.8-4,  and  2.8-5-  from  Section 
2  may  be  used  with  the  information  of  the  subsequent  sections  to 
obtain  the  desired  problem  formulation. 

Section  3  presents  the  definitions  of  and  the  formulae  applicable  to  the 
applied  loads,  which  are  shown  on  the  right  liand  side  of  the  equations  of 
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motion  (Equations  2.8-3,  2.8-4,  and  2*8-5).  The  importance  of  this  section 
of  the  report  lies  in  the  procedures  for  the  proper  introduction  of  these 
forces  into  the  system  of  equations.  The  specific  information  on  procedures 
for  calculating  applied  loads  is  based  on  the  best  currently  available  data. 
These  procedures,  howc/er,  may  not  always  exactly  fit  the  physical  condition 
being  simulated.  It  is  necessary,  then,  that  the  user  decide  upon  rhe  applic¬ 
ability  of  the  given  expressions,  and,  if  necessary,  modify  these  expressions 
or  completely  define  new  expressions  applicable  to  the  vehicle  being  inves¬ 
tigated.  In  the  latter  case,  the  information  of  Section  3  becomes  the  guide 
to  be  used  in  the  formulation,  and  provides  the  necessary  framework  for  prop¬ 
er  introduction  of  newly  defined  forces  into  the  equations  of  motion. 

The  procedures  whereby  the  results  of  Sections  2  and  3  may  be  manipulated 
for  use  in  the  landing  analyses  of  vehicles  of  current  interest  are  presented 
in  Section  4,  Vol.  II.  It  is  expected  that  Section  4,  Vol.  II  should  be 
sufficiently  complete  to  make  the  report  valuable  even  without  a  good  under¬ 
standing  of  the  Lagrangian  formulation  of  mechanics.  Section  4,  Part  II 
explains  the  proper  procedures  for  utilizing  the  method. 

Vol.  II  is  concluded  vrith  the  presentation  in  Section  5  of  several  nu¬ 
merical  solutions  of  the  example  problems  in  Section  4.  These  problems  are 
intended  to  exemplify  the  salient  features  of  the  method.  To  provide  ease  of 
understanding,  an  attempt  has  been  made  to  present  specific  aspects  of  the 
method  in  each  example  so  that  the  user  will  not  be  confronted  with  an  ex¬ 
cessive  number  of  new  concepts  in  a  single  problem.  After  an  investigation 
of  the  information  in  the  report,  (Vol.  I  and  II),  the  user  should  be  able  to 
proceed  to  the  solution  of  more  detailed  and  extensive  problems. 

Providing  a  complete  and  rational  method  for  predicting  alighting  gear 
loads  during  landing  impact  has  necessitated  the  presentation  of  many  diverse 
but  related  topics.  Some  of  these  are  basic  to  the  development;  for  example, 
distributed  mass  and  rotational  inertias  and  distributed  bending  and  torsional 
stiffness  of  all  parts  of  the  vehicle,  as  represented  by  the  various  A  and  K 
matrices  in  the  equations,  are  inherent  in  the  mavheroatical  model  used  for 
the  development,  and  thereby  vehicle  flexibility  is  included.  The  exact  con¬ 
figuration,  shape  and  dimensions  of  all  basic  components  of  the  vehicle  being 
studied,  as  well  as  those  of  the  composite  vehicle,  are  required  for  the  pur¬ 
pose  of  defining  the  various  forces.  The  mathematical  model  is  sufficiently 
general  to  allow  use  of  one,  two,  three,  four  or  more  separate  alighting  gears 
for  vehicle  support.  The  mathematical  model  also  provides  the  freedom  re¬ 
quired  for  establishment  of  variations  in  coefficients  of  other  parameters, 
and  of  any  desired  combinations  of  inital  conditions  on  displacements,  rates, 
or  accelerations;  for  example,  rigid  body  pitch  (attitude)  and  pitch  rate  may 
be  specified  simply  as  initial  conditions,  and  the  desired  initial  value  of 
pitch  acceleration  may  be  obtained  by  the  proper  unbalance  of  moments  in  the 
pitch  degree  of  freedom  equation  of  motion.  (This  same  statement  is  applicable 
to  vertical,  longitudinal  and  lateral  translations  and  to  roll  and  yaw  rota¬ 
tions.  ) 

To  provide  a  liandy  reference  for  inital  utilization  of  this  report, 
some  of  the  additional  topics,  which  are  not  general  throughout  the  report, 
but  which  are  worthy  of  special  note,  are  presented  in  Table  1,  page  seven. 

This  table  shows  the  page  location  for  some  of  the:general  topics  as  well 
as  for  such  special  topics  as  friction  forces  and  mechanical  springs. 
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While  it  is  not  feasible  to  cover  in  such  a  table  all  details  included  in  the 
report,  most  major  special  topics  are  enumerated. 
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TABLE  1 

LOCATION  OF  TOPICAL  DISCUSSIONS 


A.  CHARACTERISTICS  OF  THE  RIGID  VEHICLE 

1.  Development  of  Equations  of  Motion -  8 

2.  Mass  and  Rotational  Inertias  in  the  Rigid  Body  24 

Equations  of  Motion - 

B  CHARACTERISTICS  OF  FLEXIBLE  VEHICLE 

1.  Mathematical  Model -  8 

2.  Equations  of  Motion - — -  25 

C.  THE  INTERMEDIATE  STRUCTURE 


1.  General  Equation  of  Motion  Showing  Distributed  Mass 

and  Rotational  Inertias,  and  Bending  and  Torsional 
Stiffnesses -  23 

2.  Alighting  Gear  Restraints,  Including  Full  Cantilever 

with  One  and  Two  Restraints -  107 

3*  Angular  Orientation  of  Co or  with  Respect  to  Vehicle 

and  to  Alighting  Surface -  11 

4.  Articulated  Alighting  Gear  (Volume  II ) -  232 

5.  Semi-Articulated  Alighting  Gear  (Volume  II) -  24l 

6.  Characteristics  of  Shock  Absorbing  Devices - 

a.  Springs,  Linear  and  Non-Linear -  94 

b.  Energy  Absorption  Devices 

(1)  Hydraulic  Dampers -  87 

(2)  Crushable  or  Yielding  Materials -  105 

(3)  Release  of  Compressed  Fluids -  106 

c.  Strut  friction -  101 

d.  Bottoming  Forces -  105 

e.  Extension  Stroke  Forces -  92 


D.  SURFACE  CONTACTING  COMPONENTS 

1 .  Surface  Pads - 

2.  Tires,  Including  Load-Stroke  Characteristics  and 


Equivalent  Damping -  6l 

3.  Tire  Friction  Forces -  64 

4.  Spike-Soil  Forces -  70 

5.  Skid  and  Ski  Forces 

a.  Hard  Ground  Forces -  75 

b.  Soft  Ground  Forces -  78 

6.  Gas  Filled  Bag  Forces . 82 

7.  Equation  for  Bogie  Motion - - - - —  30 


8.  Consideration  o.f  Off -Set  or  In-Line,  Dual  or  Multiple 
Contacting  Elements - - — 


7 


SECTION  2 


FORMULATION  OF  THE  EQUATIONS  OF  MOTION 


2.1  GENERAL 


The  formulation  of  the  equations  of  motion  given  in  this  section,  al¬ 
though  somewhat  formidable  algebraically,  can  be  described  in  simple  terms. 

From  the  definition  of  inertial  and  body  coordinate  systems  the  kinetic  and 
potential  energies  for  the  mathematical  model  under  consideration  are  written 
in  the  body  coordinate  system.  These  energy  expressions  are  then  substituted 
into  a  modified  expression  of  Lagrange's  equations,  yielding  the  equations  of 
motion  of  the  vehicle  during  landing  impact.  The  desire  for  generality  and 
compactness  in  the  final  equations  has  dictated  that  the  variables  chosen  to 
define  the  vehicle  motion  be  those  expressed  in  the  accelerating,  or  body, 
reference  system.  This  choice  of  variables  requires  that  the  operations 
associated  with  the  usual  form  of  Lagrange's  equations  be  transformed  to 
operations  in  the  accelerating  reference  frame. 

In  order  that  continuity  may  be  maintained,  many  of  the  detailed  steps 
of  the  derivation  have  been  deferred  to  the  appendices. 

2.2  MATHEMATICAL  MODEL  OF  THE  VEHICLE 

The  vehicle  under  consideration  is  composed  of  n  components,  each  of 
which  may  be  elastic.  Each  component  is  assigned  a  subscript  i  ,  so  that 
i  a  1,  2,  — ,  N.  The  subscript  which  is  unity  refers  to  the  main  com¬ 
ponent  or  fuselage.  The  remaining  subscripts  refer  to  wings,  tails,  gears, 
etc.  Each  of  these  components  has  a  defined  volume,  a  mass  density,  and  an 
initial  orientation  in  space. 

The  motion  of  the  continuous  distribution  of  mass  requires  an  infinity  of 
variables  to  be  exactly  described.  It  is  approximated  by  the  motion  of  a  dis¬ 
crete  number  of  points  called  panel  points.  Mass  and  stiffness  properties  are 
assigned  to  the  panel  points  by  some  interpolation  scheme  such  that  the  energies 
of  the  continuous  system  and  the  panel  point  model  are  equivalent. 

The  motion  of  the  panel  points  will  be  defined  relative  to  the  body  axis, 
which  are  not  fuxed  in  the  inertial  frame  of  reference.  These  motions  must  then 
be  described  relative  mo  the  inertial  frame  in  order  to  formulate  the  kinetic 
energy  of  the  system. 

2 . j  NOTATION 


The  notation  in  this  report  is  made  consistent  with  that  in  the  rield  of 
aerodynamics  where  possible.  The  size  of  the  report  indicates  that  duplication 
of  symbols  may  occur.  The  forms  to  he  used  will  be  consistent  throughout  the 
report,  and  symbols  used  more  than  once  are  defined  for  each  usage. 
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General  notation  as  used  in  this  report  is  as  follows: 


Vectors  -  Vectors  are  denoted  by  a  double  bar  on  the  left  edge  of 
the  symbol;  IP ,  IR ,  IL . 


Matrices  -  Matrices  are  denoted  by  a  symbol  enclosed  by  braces  o: 

brackets,  the  distinction  between  types  are  listed: 


c  r 

[  ]'- 
u- 

t 

fij" 


A  rectangular  matrix;  the  number  of  rows  and  columns  depends 
on  the  particular  matrix 

A  rectangular  matrix;  the  transpose  of  the  matrix  indicated 
by  the  brackets  and  enclosed  symbol 

A  column  matrix;  a  single  column  of  elements  as  indicated 
by  an  enclosed  symbol 

A  row  matrix;  the  transpose  of  the  indicated  column  matrix 

The  identity  matrix;  a  matrix  with  unity  in  the  diagonal 
positions  and  zeroes  off  the  diagonal 


A  column  matrix  composed  of  the  indicated  column  matrices 


Time  derivatives  -  A  partial  derivative  of  a  variable  with  respect  to 

time  is  indicated  by  placing  a  dot  above  the  variable,  once  for 
each  time  it  is  differentiated; 


Primes  -  When  placed  on  a  matrix,  the  matrix  transpose  is  indicated. 

When  placed  on  a  variable  or  subscript  to  a  variable,  the 
quantity  is  considered  to  be  written  in  a  component  co¬ 
ordinate  system. 

Particular  cases  of  the  matrix  notation  involve  combinations  of  the  above 
notation.  These  cases  are  consister.*  in  notation  and  generally  will  not  be 
explained  in  detail. 

Symbols  and  notation  of  components  of  the  vehicle  are  considered  in  the 
definition  of  coordinate  systems  and  vectors. 

2.4  COORDINATE  SYSTEMS 


Three  distinct  coordinate  systems,  inertial,  body,  and  component,  are  used 
to  define  the  motion.  Each  of  the  coordinate  systems  is  defined  by  a  right  handed 
triad  of  mutually  orthogonal  unit  vectors  which  specify  the  directions  of  the 
coordinate  axes. 
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The  inertial  frame  or  ground  frame  of  reference  is  specified  by  the  unit 
vectors  E  ,  <3  ,  fl<  .  The  unit  vector  K  is  normal  to  the  horizontal 
ground  plane  and  is  positive  downward.  The  unit  vectors  H  >  are  in 
the  ground  plane. 

.  The  principal  or  body  coordinate  system  is  specified  by  the  unit  vectors 
fl  ,  d)  ,  IK  .  The  orientation  of  these  vectors  is  specified  initially  by 
the  principal  axes  of  inertia  of  the  undeflected  vehicle,  and  thereafter  by  the 
equations  of  motion.  The  unit  vector  IK  is  considered  positive  downward  and 
the  unit  vector  j  positive  forward,  where  those  directions  have  meaning. 

Component  coordinate  systems  are  specified  by  the  unit  vectors  ti  >  , 

IK  £  ,  The  subscript  refers  to  the  particular  component.  The  origins  of  these 
component  coordinate  systems  are  coincident  with  that  of  the  principal  or  body 
coordinate  system  and  are  fixed  relative  to  the  principal  or  body  coordinate 
system.  These  coordinate  systems  are  oriented  in  each  component  such  that  the 
motion  of  that  component  is  as  simple  to  describe  as  possible.  This  concept 
is  discussed  as  various  component  problems  arise. 

Transformations  between  the  coordinate  systems  are  given  by 


Ml 

m 

j> 

■  *  [n]  ■ 

4 

JK. 

UKj 

rn 

rn 

Inertial  to  vehicle  body 
coordinate  system 


Component  to  vehicle  body 
coordinate  system 


where  the  elements  of  the  transformations  are  the  direction  cosines  between  the 


coordinate  axes,  given  by 

[p]  = 

'hi 

a  •  3T 

a  •  K" 

<11  •  H 

v!)  'JT 

Jl  *  \K 

JK-  n 

IK' J 

JK  *  DC  _ 

1*1“ 

"a  •  i'i 

4  '4 

i  • 

Ji  •  it 

Jj  -4 

ji  -  ii i'i 

JK  *  i'i 

IK  •  K* 

iK'lK't 
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Since  this  method  is  primarily  concerned  with  the  cantilevered  type  of  landing 
gears,  it  is  assumed  for  the  general  formulation  that  the  component  transformation 
matrices  [Y]^  are  independent  of  time.  This  assumption  will  then  allow  simple 
descriptions  of  component  motions  if  the  orientation  of  the  components  are  fixed 
relative  to  the  body  cc'xdinate  system.  This  simplification  keeps  the  size  of 
the  formulation  within  reason  while  not  compromising  the  accuracy  fo'~  the  in¬ 
tended  application.  In  Section  Four  the  logical  extension  of  the  equations  of 
motion  to  include  articulated  landing  gears  will  be  demonstrated  by  example. 

The  matrix  [p]  is  the  Eulerian  transformation  from  the  ground  coordinate 
system  to  the  body  coordinate  system.  It  is  given  by  Equation  2.7-3,  and  is 
time  dependent. 

2.5  POSITION  AND  VELOCITY  VECTORS 

Any  formulation  of  a  dynamics  problem  requires  that  the  position  of  all 
masses  be  known  relative  to  an  inertial  frame.  With  reference  to  Figure  1, 
the  following  vectors  must  then  be  defined. 

IP  -  The  position  vector  of  the  origin  of  the  principal  or  body 
coordinate  system: 


lR“XII-*-YvP+ZG< 

IP  -  The  velocity  of  the  origin  of  the  body  coordinate  system: 

ir  -~xi  +  y jr  +  i  ik 
IR  =  /14 1  *  J|  +4 j^jk 


IL  -  The  position  vector  from  the  origin  of  the  body  coordinate 

system  to  the  unueflected  position  of  an  elemental  volume  dV: 


lL-xi  +  ^d)+/^!K 

IP  -  The  position  vector  oi  an  elemental  volume  dV  from  the  origin 
of  the  ground  or  inertial  reference  system. 
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-  The  angular  velocity  of  the  body  expressed  in  the  body 
coordinate  system 

a  =  nx  a  +  ik 

[D  -  The  total  displacement  of  an  elemental  volume  relative  to 
r  the  body  coordinate  system,  expressed  in  the  component  co¬ 
ordinate  system 

!P  -  Px.  S'  +  Pf-  dl'  +  f>,  IK; 

In  addition,  there  are  several  vectors  which  must  be  defined  for  use 
later  in  the  formulation. 


The  displacement  of  a  volume  element  dV  from  its  undeflected 
position  defining  component  rigid  body  displacement  with 
respect  to  the  body  axes,  expressed  in  the  component  coordinate 
system  in  which  the  motion  is  most  easily  described,  (referred 
to  as  "delta''  motion). 


n  -  The  displacement  of  a  volume  element  in  the  i  -th  component 
''  ijf  due  to  displacement  of  the  -v  -th  componert,  .  The  letter 

^  here  simply  refers  to  tne  particular  component  to  which  the 
i,  -th  component  is  attached. 

ip  ®  -  The  displacement  vector  of  an  elemental  volume  due  to  elastic 
*  deformations  of  the  {  -th  component  only, 


The  total  displacement  of  an  elemental  volume  is  then  the  sum  of  the  latter 
three  vectors. 

IP;  =  Iff  -  IPif  -  *'i 

For  the  main  component,  or  fuselage,  there  is  no  "delta"  motion  or  motion  caused 
by  other  displacements,  and 

IP,  Iff 

The  definition  of  the  principal  or  body  coordinate  system  requires  that  it 
be  initially  the  principal  axis  of  inertia  system,  so  that  by  defintion 


f  U-KmVKM-)  Jv  ‘  ° 

where  p  (x,  ij,  y)  is  the  mass  density  of  the  vehicle  and  V  is  the  total  volume. 
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2.6  KINETIC  AND  POTENTIAL  ENERGY  OF  THE  SYSTEM 


The  kinetic  and  potential  energy  of  the  system  can  now  be  written. 


The  kinetic  energy  of  any  system  may  be  defined  by 


(2.6-1) 


where  the  vector  /r  is  the  position  vector  of  an  elemental  mass 
relative  to  an  inertial  frame  of  reference.  From  Figure  1,  the  position 
vector  is  seen  to  be 


ir=  ir  +  il  +  ip 


(2.6-2) 


In  vector  calculus,  it  is  shown  that  the  total  time  derivative  of  a 
vector  expressed  in  an  accelerating  coordinate  system  is  found  from  the 
operation 


_d_ 

dt 


d>t 


+  fix 


(2.6-3) 


This  theorem  together  with  the  definitions  of  the  above  vectors  will  be 
used  to  find  the  velocity  of  the  elemental  mass  relative  to  the  inertial 
frame.  The  vector  JR  is  a  vector  in  the  inertial  frame  of  reference,  so 
that 


dJR  _  fR  (2-6-4) 

dt 


The  vector  IL  is  expressed  in  the  body  coordinate  system,  so  that  it  is 
dependent  on  body  orientation,  but  it  is  not  explicity  dependent  on  time. 
Then 


dlL 

dT  *  n  X  il 


(2.6-5) 


The  total  displacement  vector  of  an  elemental  volume  is  dependent  on  time 
and  body  orientation; 

=  P  +  0.  X  F  (2.6-6) 

dt 

The  velocity  of  an  elemental  volume  with  respect  to  the  ground  reference 
system  is  then 


=  jr  +  ft  x  jl  +  ip  t  n.  x  p 

dt 


.14 


The  square  of  this  quantity  is 

dir.  dr  =  iR1  '■Uh.flL?  +  f?2’  +  UUIPf 

dt  dt 

+  2IR«JL*1P  t  2  (R.  |P  +2 
+  2  ft*  1L^  IP  +  2  ft*  IP*  IP 

-f-  2.  (ift-x  IL)  -  (^L x  rp)  (2-6-8) 

This  form  is  substituted  into  the  expression  for  kinetic  energy  given  by 
Equation  2.6-1.  The  volume  integration  of  the  above  terms  containing  mnel 
point  displacements  or  velocities  is  broken  into  integration®  over  each 
volume  with  the  resulting  summation  indicated. 


1"  =  k'-pi'i  +  ^/WiMYJv 

♦niAv* 

i-=l  l  As  A. 


•h  fP'JJL*/  IP;  /*J|/  t  iPj/lPi/odV  +^/(lx  ip.^dv 


The  term 

IR  ll-^dv  (2.6-10) 

^  V 

does  not  appear  due  to  the  definition  of  IL  . 

The  desired  form  for  the  kinetic  energy  is  obtained  by  expanding  the  vector 
quantities  in  Equation  2.6-9  into  component  form  and  performing  the  indicated 
operations.  The  integrals  are  then  evaluated  to  give  the  form  of  the  kinetic 
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energy  consistent  with  the  mathematical  model  described  in  Section  2.2.  The 
detailed  manipulations  required  to  arrive  at  the  desired  final  expression  for 
the  kinetic  energy,  shown  in  Equation  2.6-11,  are  outlined  in  Appendix  A. 


(Summation  continues  on  next  page) 


JL.O 


(2.6-11) 


w  £“  <£• 

q?  <=»  CJ 


(2.6-11) 


IT 


It  should  be  noted  that  the  subscript  i  does  not  appear  on  each  term, 
but  rather  on  the  entire  bracket  of  terms  to  be  summed.  Thus  all  primed  sub¬ 
script  terms  should  have  the  proper  component  subscript  affixed,  as  should  each 
coordinate  transformation  CT]  .  The  vehicle  velocity  components  are  written  in 
the  body  coordinate  system  and  the  terms  involving  columns  of  panel  point  dis¬ 
placements,  {.  P}  ,  and  their  time  derivatives,  are  written  in  the  various 
component  coordinate  systems,  as  is  evidenced  by  the  primed  subscripts.  The 
transformation  of  the  component  coordinate  terms  into  the  body  coordinate 
system  is  evident  due  to  the  presence  of  the  coordinate  transformation  matrices. 
The  matrices  lA]  and  [nr]  are  defined  in  Equations  2.1-11,  12. 

The  matrices  [A]  with  various  subscripts  are  called  mass  matrices.  They 
have  different  forms  depending  on  the  type  of  integral  which  is  evaluated.  Each 
arises  from  some  numerical  scheme  which  relates  the  discrete  form  of  the  kine*ie 
energy  to  that  of  the  continuous  physical  system,  consistent  with  the  definition 
of  the  mathematical  model. 

It  is  assumed  that  the  potential  energy  in  the  vehicle  arises  only  from 
structural  deformations  in  the  linear  range.  This  requires  that  elastic  dis¬ 
placements  remain  small;  i.e.,  that  the  relative  displacements  of  adjoining 
panel  points  be  small  compared  to  the  distance  between  them.  The  potential 
energy  may  be  expanded  in  a  Taylor's  series  in  panel  point  elastic  displace¬ 
ments.  If  it  is  expanded  about  a  point  of  minimum  potential  (the  equilibrium 
position) ,  and  the  potential  at  that  point  is  arbitrarily  set  equal  to  zero, 
the  constant  and  linear  terms  in  elastic  displacements  will  be  zero.  If  all 
terms  higher  than  second  order  arc  neglected,,  a  quadratic  form  for  the  potential 
energy  results: 


[K*v] 

[K  *y] 

[K«yl 

m 

V  4LW 

[K  |'x'] 

'{$• 

M 

[K-J'x] 

iKvjd 

The  matrices  indicated  here  are  called  stiffness  matrices,  and  are  derived 
from  the  geometry  and  physical  characteristics  of  the  various  components  of  the 
systar  They  embody  not  only  the  distributed  stiffness  of  a  given  component,  but 
also  the  restraints  due  to  supporting  members  between  components.  Their  determina¬ 
tion  thus  involves  both  a  detailed  structural  analysis  and  some  sort  of  numerical 
technique,  sc  that  the  continuous  structure  is  properly  represented.  This  sub¬ 
ject  is  too  broad  to  be  included  in  this  report;  the  reader  unfamiliar  with  this 
aspect  of  the  analysis  may  consult  references  >,  5,  18  and  19. 

The  above  forms  for  the  kinetic  and  potential  energy  will  be  substituted 
into  Lagrange's  equations,  which  are  to  be  developed  next. 
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£.7  IAGRAilGB'S  EQUATIONS 


A  derivation  of  the  basic  form  for  Lagrange's  equations  may  be  found 
in  reference  1.  A  summary  of  the  concepts  involved  in  their  derivation  is 
necessary  in  order  to  understand  the  contents  of  this  section. 

Basic  to  the  derivation  is  the  requirement  that  the  position  and  orienta¬ 
tion  of  every  particle  be  specified  relative  to  a  fixed  frame  of  reference. 
This  is  accomplished  in  terms  of  any  set  of  variables;  as  these  may  be  dis¬ 
tances,  velocities,  angles,  or  perhaps  less  familiar  quantities,  they  are 
referred  to  as  generalized  coordinates.  The  appearance  of  Lagrange's  equa¬ 
tions  varies  among  authors.  A  convenient  form  for  use  here  is 


_d 

dt 


(2-7-1) 

L  =  i,Z,  -Mz 


The  quantity  L  is  known  as  the  Lagrangian,  the  difference  between  the 
kinetic  and  potential  energy:  L  =  T  -  U.  Here  the  are  the  generalized 
coordinates.  The  minimum  number  of  coordinates  which  may  be  chosen  is  dic¬ 
tated  by  the  number  of  degrees  of  freedom  allowed  in  the  motion,  and  the 
number  of  constraints  imposed.  The  quantities  Q ■  ?re  the  associated  gen¬ 
eralized  forces;  they  may  actually  be  forces,  moments,  etc.  These  gener¬ 
alized  forces  are  assumed  known  in  terms  of  .he  generalized  coordinates. 

The  quantities  Fj  are  a  set  of  constraint  relations  among  the  general¬ 
ized  coordinates/  If  there  are  M  -  of  these  constraints,  where  M,  is 
the  number  of  generalized  coordinates,  then  there  are  degrees  of  freedom 
for  the  system.  The  Lagrangian  undetermined  multipliers,  6~j  ,  are  the 

proper  functions  to  cause  each  term  on  the  right  hand  side  of  the  equations 
to  be  a  constraint  force  which  is  consistent  with  the  constraint  relation¬ 
ships  Fj  .  As  there  are  M,  generalized  coordinates  and  M,  -  14,  Lagrange 
multipliers,  and  only  M  Lagrange  equations,  it  is  necessary  to  solve  the 
set  of  M -  M,  constraint  relations  simultaneously  with  Lagrange's  equations. 
These  will  be  discussed  later. 


For  the  sake  of  algebraic  simplicity,  the  kinetic  energy  of  the  system 
was  written  in  terms  of  the  body  coordinate  system  variables.  Since  this 
system  moves  with  the  body,  it  does  not  specify  the  position  and  spatial 
orientation  of  every  particle  in  the  body,  contrary  to  the  assumptions  from 
which  the  above  statement  of  Lagrange ' s  equations  was  derived .  It  is  then 
necessary  to  modify  either  the  expression  for  the  kinetic  energy  so  as  to 
be  applicable  to  Lagrange's  equations  as  expressed,  or  to  modify  Lagrange's 
equations  to  utilize  the  kinetic  energy  in  its  present  form. 

An  investigation  of  these  alternatives  lias  shown  that  much  of  the 
cumbersome  algebra  can  be  avoided  if  Lagrange's  equations  are  modified  to 
accept  the  kinetic  energy  in  its  present  form. 
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Before  proceeding,  the  Eulerian  transformation  between  the  ground 
coordinate  system  and  the  body  coordinate  system  will  be  established.  The 
transformation  was  defined  such  that 


A 

'  jl  * 

UK. 


(2.7-2) 


The  form  for  the  individual  elements  in  the  transformation  matrix  will  depend 
on  the  definitions  of  the  Euler  angles  to  be  used,  that  is,  the  sequence  of 
rotations  which  take  the  vehicle  from  the  orientation  of  the  ground  reference 
system  to  its  instantaneous  space  orientation.  The  choice  of  these  rotations 
is  arbitrary.  A  convenient  sequence  of  rotations  for  landing  problems  is 
defined  as  follows: 


(1)  Rotate  the  system  through  the  angle  "'|r  about  the  Z  axis  in  the  right- 

handed  manner  (X  into  Y)  to  produce  the  system  (X  ,  Y  ,  Z  ) . 

Ill 

(2)  Rotate  the  system  (X  ,  Y  ,  Z  )  through  the  angle  9  about  the  Y^  axis 

in  the  right-handed  ^manner  (Z  into  Xt  )  to  produce  the  system 

(x^,  z2). 

(3)  Rotate  the  system  (X.,  Y  ,  Z  )  through  the  angle  0  about  the  X.  axis 

in  the  right-handed  ^manner  d  (Y  into  Z„)  to  produce  the  system^ 

(*^)- 

The  coordinate  transformation  is  given  then  by 
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(2.7-3) 


The  transformation  of  the  body  velocity  as  expressed  in  the  ground  reference 
coordinate  system  to  the  expression  in  the  body  coordinate  system  is  given  by 
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The  transformation  from  Eulerian  angle  time  derivatives  to  angular 
velocities  about  the  body  coordinate  axes  is  given  by 
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O 
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where  the  transformation  matrix  is 

-Sin  ■&  o 

[R]  =  Cos  (p 

CasCpCoS-O-  -  Sin.  Ip 

The  particular  property  of  this  choice  of  Eulerian  angles  useful  ii 
landing  problems  is  the  existence  of  the  inverse  of  the  matrix  [R]  when 
all  of  the  Euler  angles  are  zero.  Conventional  airpJaae  axes  are  thus 
chosen  to  represent  the  body  coordinate  system,  and  the  vehicle  orientation 
is  specified. 


->) 


A  set  of  generalised  coordinates  suitable  to  completely  specify  the 
motion  could  be  the  Euler  angles,  the  components  of  the  position  vector  to 
the  origin  of  the  body  coordinate  system,  and  the  panel  point  displacements. 
The  Lagrangian  then  would  be  giver*,  by 


L  -  L(x,nz,  V,-e,4>,f  px'?£  Ps'},{p?-})  U- 7-7) 

The  statement  of  Lagrange's  equations  using  i-hs^e  generalized  coordinates 
follows  immediately  from  the  general  form  of  Equation  2.7-1*  Thus 


the  last  set  of  equations  are 


written  once  for  each  component. 


a 


The  kinetic  energy  could  be  transformed  into  an  expression  relating  these 
variables,  and  the  indicated  operations  performed.  The  algebra  involved  in 
obtaining  the  desired  result,  however,  may  be  lessened  considerably  by  trans¬ 
forming  the  operations  in  these  equations  to  operations  on  the  variables  in 
which  the  kinetic  energy  is  already  expressed.  The  derivation  of  Lagrange's 
equations  in  "modified  form  is  shown  in  detail  in  Appendix  B.  In  terms  of 
the  matrices 
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the  modified  statement  of  Lagrange's  equations  is  shown  to  be 
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2.8  THE  EQUATIONS  OF  MOTION 


The  desired  equations  of  motion  may  now  be  obtained  by  the  substitution 
of  the  forms  for  kinetic  and  potential  energy  into  the  set  of  Lagrange  Equa¬ 
tions  2.7-13*  14,  15*  The  operations  indicated  there  are  somewhat  lengthy, 
and  are  relegated  to  Appendix  C. 

The  final  equations  of  motion  are  presented  as  three  sets  of  equations, 
each  defining  one  of  the  vectors  describing  the  motion.  The  components  of 
the  linear  velocity  vector  if?  are  governed  by  Equation  2.8-3,  and  those  of 
the  angular  velocity  ID.  by  Equation  2.8-4.  The  total  displacement  vector  IP 
is  governed  by  Equation  2.8-5-  These  equations  describe  the  motion  relative 
to  the  body  coordinate  axes,  which  constitute  an  accelerating  frame  cf  reference. 
The  motion  relative  lo  the  inertial  or  fixed  frame  of  reference  is  obtained  by 
transformation  of  the  body  velocities,  making  use  of  Equations  2.7-4,  5*  These 
velocities  are  then  integrated  with  respect  to  time,  yielding  the  components 
of  the  position  vector  [R  which  specifies  the  position  of  the  origin  of  the 
body  coordinate  system; 


(2-8- 

of  the  body  coordinate 


and  the  Euierian  angles  vhich  specify  the  orientation 
system  b>  Equation  2.7-2 


(2.8-2) 


Up  to  this  point,  the  number  of  restrictive  assumptions  imposed  on  the 
mathematical  model  have  been  small.  The  result  is  that  a  large  number  of  terms 
appear  in  the  equations  of  motion  in  the  appendix.  In  the  interest  of  presenting 
a  relatively  compact  set  of  equations  here,  some  of  the  terms  are  deleted  from  the 


Rigid  Body  Linear  Equations 
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Rigid  Body  Angular  Equations 


flexible  Body  (Panel  Point)  Equations 


■pud  Ip,}]  I"  [  0  ] -2^Axy]  2a^#x} 

[A^]  ^  [Pj  +  [0  ]  {fy} 

MwJ  NM  2aja^i  l  o  ]J({p^ 

A*}'  l({\) 

(^%/JCr/l^^/  {|} 

(/%* riyA^-/l^)[/\y]j[{  l)  ^ 


+  (Aj-  +Ax</ly)  |Ay j  “(A^  +  Uy)  [hytfl  ('Ax'^A^A^A^  {r-.} 

i*A^+A,%)CAyJ  (A,.r^Ay)lAy#3  “(A,-  +A^  [A^l  l^.} 


'-(A^  *Jly)  [AJ'  (V-M^A*]'  (Ay  ♦iUy)[A*.J|f(x'J 
+  (Ay  +A*%)  [a^]'  -(axa  *fl})  [Ay]'  (-AX-/1^)[A  J  {*] 
}-%  +AX'Aj')  [Ay]  (A/+A^[Ay]  -(A**  *Jlf)[hyX  WJ 


25 


complete  set  of  equations  in  Appendix  C.  The  Equations  2.8-3,  5  are  identical 
with  the  expressions  obtained  in  the  appendix.  The  set  of  Equations  2.8-4 
governing  the  body  angular  velocity  are  obtained  from  Equation  C-10  by  omitting 
a  number  of  terms.  The  omitted  terms  are  of  two  types.  It  will  be  .recalled 
that  the  body  coordinate  axes  are  initially  the  principal  axes  of  inertia  for  the 
undeflected  body.  If  this  remained  true  of  the  coordinate  axes  as  the  panel 
points  move  relative  to  one  another,  the  moments  of  inertia  of  the  v-nicle  would 
still  change  due  to  the  variation  in  position  of  the  various  distributed  masses. 
These  moments  of  inertia,  which  are  dependent  on  the  panel  point  motion,  are  de¬ 
fined  by  the  inertia  matrix  in  Equation  C-9,  and  the  manner  in  which  this  matrix 
enters  the  equation  of  motion  is  indicated  in  Equation  C-10.  It  will  be  assumed 
here  that  the  terms  involving  panel  point  motions  in  the  inertia  matrix  may  be 
neglected  in  comparison  to  the  constant  term.  In  addition,  since  the  body  co¬ 
ordinate  axes  are  not  instantaneously  the  principal  axes  of  inertia,  the  inertial 
forces  associated  with  the  panel  point  motions  will  exert  a  net  torque  about 
body  coordinate  axes.  Of  the  seven  types  of  terms  which  arise  because  of  this 
effect,  only  the  first  and  third  terms  listed  in  Equation  C-10  will  be  retained 
here . 


The  deletion  of  the  above  terms  is  justified  provided  that  the  principal 
axes  of  inertia  do  net  differ  appreciably  from  the  body  coordinate  axes  through¬ 
out  the  analysis  and  that  the  moments  of  inertia  remain  nearly  constant.  For 
landing  analyses  of  conventional  vehicles  this  assumption  is  always  true.  Con¬ 
figurations  which  do  not  satisfy  these  assumptions  will  require  the  additional 
terras  to  be  retained. 

To  be  complete,  the  system  of  equations  governing  the  vehicle  motion  must 
include  the  equations  of  constraint.  These  equations  are  written  in  the  form 

(2.8-6) 

These  constraint  equations  will  arise  if  the  set  of  panel  point  displacements 
used  to  define  the  motion  are  not  independent  variables.  If  this  occurs,  Equations 
2. 8-3, 4-, 5  do  not  provide  enough  relations  to  solve  the  problem.  The  remaining 
relations  necessary  are  the  set  of  constraint  equations  which  state  the  dependence 
among  the  variables.  The  various  constraints  that  arise  in  the  class  of  problems 
considered  in  this  report  are  discussed  in  Section  3*5* 

Equations  defining  the  rigid  body  motion  of  components  with  respect  to  the 

body  coordinate  axes  (A  -  motion).  ”  — 

Although  Equations  2. 8-3, 4, 5, 6  are  a  proper  description  of  the  motion,  they 
are  now  in  the  most  convenient  form  for  general  application.  It  vdlx  be  recalled 
that  the  total  displacement  IP  of  an  incremental  volume,  dv  was  composed  of  three 
displacements: 

|p£  -  The  displacement  of  an  elemental  volume  due  to  elastic  deforma¬ 
tions  of  the  i  -th  component  only. 

(P, .  -  The  displacement  of  an  elemental  volume  in  the  i  -th  component 
due  to  the  displacement  of  the  f  -th  component,  iPj  .  ..:ere  the 
subscript  -j  refers  to  the  component  to  which  t.ie  i  -th  com¬ 
ponent  is  attached. 
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fa .  -  The  displacement  of  a  volume  element  (IV  from  its  ur.defiected 
a  position  defining  component  rigid  body  displacement  with  res¬ 
pect  to  the  body  coordinate  axes  expressed  in  the  component 
coordinate  system. 

The  displacements  tP^e  and  ILi  define  displacements  for  separate  degrees  of 
freedom;  they  are  independent  variables.  Thus,  if  the  components  of  both  of 
these  displacements  along  one  of  the  coordinate  axes  are  retained,  there  will 
result  more  variables  than  equations.  The  resulting  restriction  is  therefore 
made  that  both  of  these  contributions  cannot  occur  for  a  given  component  in  the 
same  direction.  This  will  not  constitute  any  real  restriction  on  the  use  of  the 
equations  for  the  description  of  landing  impact,  sincethe  contribution  to  gear 
loads  from  elasticity  in  the  direction  in  which  a  component  may  move  relative 
to  the  vehicle  as  a  rigid  body  is  small. 

The  panel  point  coupling  terms  in  the  rigid  body  equations  required  con¬ 
sideration  of  the  totai  displacement  [P  ;  this  requirement  and  that  of  com¬ 
pactness  dictated  the  choice  of  the  total  panel  point  displacements  as  the 
generalized  coordinates  used  to  define  the  motion.  In  a  particular  application, 
however,  the  total  displacement  must  be  broken  into  its  separate  parts.  This 
paragraph  will  derive  explicitly  the  equations  governing  A  -motion  for  two  cases 
generally  required  in  landing  impact  analysis;  motion  along  a  line  such  as  gear 
stroking,  and  motion  about  a  line  such  as  rotation  of  a  bogie  about  an  axle  or 
wheel  spinup.  Other  displacements  are  discussed  in  Section  2.9* 


In  general,  the  equations  describing  -motion  of  a  component  may  be 
derived  from  the  panel  point  equations  2.8-5  directly.  They  are  obtained  by 
expanding  the  total  displacement  into  its  separate  parts,  and  seating  to  zero 
the  elastic  displacements  of  the  component  in  the  direction  in  which  component 
rigid  body  motion  is  to  be  allowed.  All  terms  are  retained,  however,  due  to 
elastic  displacements  in  the  remaining  direction(s)  and  displacements  die  to 
displacements  in  another  component.  For  convenience,  these  inertial  coupling 
terms  on  the  component  rigid  body  motion  will  be  included  with  the  applied  panel 
point  forces  (Q  }  until  the  algebraic  manipulations  are  completed.  Thus  (3} 
will  refer  to  both  applied  and  coupled  inertial  forces. 


Component  rigid  body  motion  along  a  line 


The  component  coordinate  system  for  the  component  which  may  move  as  a 
rigid  body  along  a  line  is  defined,  such  that  one  of  the  axes  lies  along  the 
line  of  motion.  This  will  be  designated  as  the  -j' -axis.  For  that  component 
there  can  be  no  \  P  ^  according  to  the  previously  mentioned  restriction. 
Also,  since  motion  along  the  ^  -axis  of  the  member  holding  this  component 


does  not  cause  motion  of  the  component,  then 
to  the  supporting  member.  Then 


zero. 


where 


refers 


l 


(2.8-7) 
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The  component  is  not  allowed  A  -motion  along  the  other  axes,  hence 


(Mi  =  Ip^-Ip*^ 

lpA  =  ip*'}i  +  lP»'W 


Thus,  the  component  may  have  elastic  deflections  along  the  other  axes,  and 
may  move  in  those  directions  due  to  motion  of  the  supporting  component.  The 
equation  for  the  -motion  along  a  line  is  then  derived  from  the  portion 

of  the  panel  point  equations  governing  motion  along  the  y'  -axis  of  the 
i  -th  component.  Use  of  Equations  2.8-7,  8,  9  yields 


Ma;.}  =  uy  (2.8,„, 

The  A  -motion  of  each  panel  point  is  the  same  for  motion  along  a  line,  and 
the  column  [  Av]  may  be  written  in  terms  of  any  one  of  the  (n)  panel  point 
displacements;  say  the  first,  Ay 


[Ay}  =  Ay,  10 


(2.8-11) 


This  may  be  rewritten  as  a  set  of  constraint  relations 


(Ay}-Ay({l]  -  [o]  (2.8-12) 

The  only  contributions  to  constraints  on  A  -motion  along  a  line  are  then 
those  which  make  the  component  move  as  a  rigid  body,  that  is, 

W  -  l Ay]  -  Ay[i]  =  {0}  (2.8-13) 

The  constraint  term  may  be  calculated  immediately; 


+®3  +  r 


The  equation  defining  the  motion  may  be  summed  by  premultiplication  by  a  row 
matrix  of  ones,  liV 

~  1*1  t  Q }  (2.8-15) 

The  internal  constraint  forces  thus  have  no  net  effect  on  the  rigid  body 
motion  of  the  whole  component,  as  expected.  Identifying  the  component  mass  as 

Mi  =  S'VCa^i]  (2.8-16) 

the  mass  moments  relative  to  the  origin  of  the  body  coordinate  system  are 

Mi?  -  IiVIa/WI 
Ml}'  =  Vlt Ayj'lfi 

Mif  =  l»VC Ay]\'y}  (2.8-17) 

Using  Equation  2.8-5  and  expanding  the  total  displacements  in  the  desired 
directions  under  consideration ,and  setting  to  zero  the  elastic  displacement 
of  the  component  in  the  direction  in  which  component  rigid  body  motion  is 
allowed,  the  final  equation  governing  A  -motion  along  a  line  is 


+  Ay) 

+ ~  (2.8-18) 
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The  right-hand  side  is  the  summation  of  all  the  external  forces  on  the  com¬ 
ponent  along  the  y  -axis .  The  parameters  Z'  ,  &>'  ,  are  the  coordinates 
of  the  component  center  of  mass  in  its  undeflected  position.  The  left  side 
of  the  equation  is  the  product  of  the  component  mass  with  its  instantaneous 
acceleration  along  the  ^  -axis  relative  to  the  fixed  or  ground  frame  of 
reference.  In  most  cases,  when  the  component  mass  is  small  compared  to  the 
total  mass,  only  the  first  significant  terms  of  the  equation  need  be  retained. 
The  remaining  terras  are  small. 

Component  rigid  body  motion  about  a  line 

The  motion  of  a  component  about  a  line  as  a  rigid  body  is  a  bit  more 
difficult  to  forumlate  than  that  along  a  line.  For  that  reason,  the  details 
of  the  formulation  will  be  presented.  This  problem  is  an  excellent  example, 
for  the  interested  reader,  in  the  use  of  constraints. 

The  coordinate  system  for  -motion  about  a  line  will  be  defined  so  that 
the  tf'  -axis  of  the  component  system  is  parallel  to  the  axis  of  rotation;  the 
motion  is  then  in  the  x'.y  plane.  This  implies  that  {  A  single  row 

of  panel  points  is  laid  along  a  line  in  the  component,  one  of  which  is  pic¬ 
tured  in  Figure  2.  The  angle  tl  measured  in  the  right-handed  sense  from  the 
x'  -axis  defines  the  position  of  the  line  of  panel  points.  The  distances  a- 
long  the  line  from  the  axis  of  rotation  to  each  panel  point  are  arranged  as  a 
column  matrix  U}  •  The  convention  will  be  used  that  the  distance  is 
positive  when  measured  along  the  end  defining  the  angle  \  ,  and  negative  if 
along  the  other  end. 
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Figure  2.  Component  Rigid  Body  Motion  About  a  Line. 
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From  Figure  2,  the  geometric  relations  between  the  angle  V.  and 
the  panel  point  displacements  are  seen  to  be 

[&*■}  ~  S.Jt]  O  ~  (2.8-19) 

{&■}'}  =  \J]  5iN^  (2.8-20) 

These  are  a  set  of  equations  relating  a  number  of  variables,  but  they  are 
not  suitable  as  constraint  equations.  Constraint  equations  can  only  relate 
variables  already  in  use.  However,  the  value  of  any  two  of  the  displacements. 


say  A*;  and  Ay  may  be  used  as 

~  ■£]  (.  I  C.os>^  (2.8-21) 

Ay  =  i,  siN'Jfl  (2.8-22) 

to  write  a  set  of  constraints 

lAx'l  =  (2.8-23) 

lAvt  -[j,}&y  (‘-s-w 

An  additional  constraint  may  be  derived  from  the  equations  relating  A*'  and 
Ay  to  the  angle  by  eliminating  that  angle; 

A*/*  ~2  J,  A*;  +  Ay*  s  O  (2.8-25) 

A  suitable  set  of  constraints  fox  the  motion  are  then 
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(2.8-26) 
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These  constraints  simply  require  that  all  of  the  panel  Joints  rotate  to¬ 
gether  by  an  amount  determined  either  by  AX|  or  A-j,  •  T‘iey  e  su 
stituted  into  the  set  of  equations  of  motion 


The  constraint  terms  are 


(2.8-29) 

(2.8-30) 
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Substitution  of  these  terras  and  the  constraint  relations  into  the  equations 
of  motion  yields 
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-Z*,) 
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(2.8-33) 
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(2.8-34) 


The  Lagrange  multipliers  may  be  eliminated  immediately.  A  linear  combination 
of  these  equations  is  found  by  preraultiplying  the  first  by  \A/JlxV  sinT^ 
and  the  second  by  1 4/J,\f  cos’?  so  that  '  * 


■h^Tnti  ~M,)  *>N 7? 


(2.8-35) 


!?,lM  l{}^  ^  -  ill'  %]  ^ 

+  ^Tn+l  U  )  CoS^l 


(2.8-36) 
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If  the  forms  defining  Ax', and  A-j'in  terms  of  7[  are  now  substituted 
into  these  two  equations  and  they  are  summed,  the  last  undetermined 
multiplier  is  eliminated.  The  equation  then  is 


ilN^C0S>t  ^  -SI HZ7t  7[) 

+  U]  (.sin?/  Cost/  rf  Z  -  cos1^  7[) 

~  6x'l  SIN'’/  MiVSflyl  cos?/  (2.8-37) 


This  may  be  written  in  terms  of  the  moment  of  inertia  about  the  -axis 

QS 


ij 

Expanding 


+  M 


-  ~\J\  \Qx'l  ■zin'tj  ~\J  cos^f  (2.8-38) 

the  inertial  terms  is  a  rather  long  process;  the  result  is 

t  r^ilv  COST/  -ily  Sin7/][Av' SIN?/  +  .fly  COST/)  -AyJ 
ti  +/^‘  cos 7/^  +Xly(/Uy  swTf  - AJ~x'  cos?/) 

i-/lTj'(ilx'  COST/  -ily  S(nV)  ■+  -tKx'llyyp^siN'l/ 

+  PyA  cosTi  +  (•fl**  +Ay)  cost/J 

-UVlQ^^-n  (2-8-39) 


where  is  the  component  mass  and  X  as  f'le  distance  from  the  axle  along 

the  line  to  the  center  of  mass  of  the  component.  The  inertial  moments  due  to 
displacement  of  the  axle  along  the  y  -axis,  PyA  ,  are  now  easily  seen. 

This  displacement  will  be  set  equal  to  piston  CU/in  the  case  where  a  bogie 
element  is  pivoted  on  an  axle  affixed  to  the  end 'of  a  gear  piston.  Moments  due 
to  fore  and  aft  and  lateral  deflections  of  the  axle  are  neglected. 

Particular  applications  of  interest  in  this  report  are  for  bogie  motion 
and  wheel  spinup.  For  bogie  motion,  there  are  only  three  panel  points  at 
which  moment  inducing  forces  are  applied;  those  at  the  front  and  rear  wheel 
axles,  and  one  at  the  point  at  which  the  bogie  rotational  spring-damper  is 
affixed.  The  distance  from  the  bogie  axle  to  front  wheel  axle  is  ,  l^at 
to  the  rear  wheel  axle  is  Jll  (negative) ,  and  w.e  bogie  equation  of  motion  is 
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IA^'  t  J/([(Ax'co3^  -Aysi^Ax'si u7i  +Hyco5->?)  -fly] 

+  lie  I  [(4c*  *lw7l  +  ^'  cos^)  +  Ay  (/tty  sinTJ  -4Tx /  cost?) 

cost?  -Ay-s/itfi)  +  (iy  f  Ay  Ay)  R^S(n7Z 

+  Py,  <«S^  +  (A*^  4- Ay)  PyA  COST?]  -  N^,^) 

~[Ui  ^x,- +  AQxj'jsiNTf  +^iQy,+-fz^^)cis?2]  , 


where  |M (or  "^)is  the  restraining  moment  of  the  rotational  spring-damper  of 
the  bogie  element. 

For  a  round  disc -like  object  rotating  through  large  angles  about  its 
central  axis,  such  as  a  tire  and  wheel  rotating  about  an  axle,  the  applied 
forces  do  not  act  directly  on  the  line  of  points  which  rotates  with  the 
wheel.  Therefore  to  define  properly  the  panel  point  forces  would  require 
additional  complexity.  It  is  noted,  however,  that  the  right  hand  side  of 
the  equations  must  be  the  external  moment  which  causes  the  wheel  to  rotate, 
and  it  can  be  simply  stated  as  such.  The  distance  X  is  zero  for  a  balanced 
wheel,  and  the  inertial  terms  in  body  angular  motion  are  small  compared  to 
wheel  angular  acceleration.  Hence  the  equation  for  wheel  angular  motion  is 

IA  "H  ~  EXTERNA!.  MOMENT  =  N  5U  (2.8-M) 

2.9  TRANSFORMATION  TO  MODAL  COORDIHATES 


The  formulation  of  the  equations  of  motion  as  expressed  in  2.8-3,  4,  5 
was  performed  with  panel  point  displacements  ($Px'l«^Py1  \P* ])  as  total  dis¬ 
placements  relative  to  the  ur.deflected  vehicle.  This  was  necessary  for 
clarity  and  brevity  of  the  formulation,  and  does  not  compromise  the  generality 
of  the  equations.  In  principle,  these  equations  may  be  solved  to  describe  the 
motion  during  landing  of  any  vehicle.  In  actual  practice,  it  is  desirable  to 
keep  the  number  of  variables  to  a  minimum  while  retaining  as  high  a  degree  of 
accux-’.c,,  as  possible.  It  is  not  generally  desirable  to  simply  reduce  the  num¬ 
ber  of  panel  points  until  the  total  number  of  variables  is  sufficiently  small. 
A  more  economic  approach  is  found  in  the  transformation  to  modal  coordinate',. 
The  elastic  motion  of  a  body  may  be  expressed  quite  accurately  by  retaining 
only  those  modes  whose  corresponding  frequencies  are  in  the  main  portion  of 
the  frequency  spectrum  of  the  forces  on  the  body.  Since  forces  in  landing 
problems  are  not  usually  composed  of  high  frequencies,  the  first  few  modes 
are  sufficient.  However,  the  mode  shape  of  an  entire  vehicle  is  nou  a  con¬ 
cept  too  useful  here,  rather,  the  mode  shapes  of  individual  components  are 
considered.  The  reason  for  this  is  fairly  simple.  If  the  airplane  ving 
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mass  is  small  compared  to  the  fuselage,  the  wing  elasticity  will  not 
contribute  to  landing  loads  and  may  be  omitted.  One  the  other  hand, 
the  wing  mass  may  be  the  greatest  part,  as  in  a  flying  delta,  and  fuse¬ 
lage  flexibility  may  be  omitted.  The  particular  problem  at  hand  will 
then  dictate  which  of  tne  effects  to  include.  With  these  ideas  in  mind, 
it  is  obvious  that  the  Equation  2.8-5  cannot  be  transformed  into  modal 
coordinates  immediately,  for  the  coordinates  will  depend  upon  the  par¬ 
ticular  problem. 

Mede  shapes  will  be  defined  for  the  following  set  of  constraints. 

(1)  Fuselage  modes  with  constrained  rigid  body  motion,  /£±  -motion, 
and  elastic  motion  of  all  but  the  fuselage. 

(2)  Wing  modes  with  constrained  rigid  body  motion,  ^  -motion,  and 
elastic  motion  of  all  but  the  wing,  unless  gears  are  on  the  wing. 

(3)  Modes  of  all  minor  elements  uith  all  motion  constained  except 
the  elastic  motion  of  that  element,  unless  another  element  is 
attached  to  that  element  and  not  the  main  body. 

Consider  the  problem  of  fuselage  modes.  The  displacement  of  all  panel  points 
on  the  wings  may  be  written  directly  in  terms  of  displacements  of  the  panel 
points  on  the  fuselage  where  the  wings  are  attached.  Normally, 


W  (Hij  rtR?n  riA„.y 

fR,'}  =  fU up,!  *  +  m*]’ 

iyJ  WF  Uyi  UMl 


(2.9-1) 


Since  wing  elastic  and  delta  motion  are  constrained,  the  constraint  equa¬ 
tions  are 


rpv)i  ftR.11 

W  -[TwdW-  = 

M„  ml 


(2.9-2) 


for  each  wing.  The  fuselage  panel  point  equation  is 


[feJ 
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*  to 


(2.^-3) 
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und  each  wing  panel  point  equation  is 


Suppose  that  all  other  masses  such  as  tails,  control  surfaces,  etc.  are 
small.  A  set  of  constraints  as  in  equation  2.9-2  is  written  for  boin  wings; 
the  equations  are  identical  if  the  wings  are  mirror  images.  The  right-hand 
side  of  2-9-b  is 


I  «i  =  VU 


7 


(2.9-6) 


and  the  right-nana  side  of  equation  2.9-3  is 


If  equations  2.9-5;  b,  7  are  substituted  into  equations  2.9-3;**  and  the 
contraints  are  used  to  write  wing  panel  point  displacements  in  terms  of 
fvsc'ege  panel  point  displacements,  the  equations  become 
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Premulti plication  of  equations  2 -9-9 >  10  hy  an<^  addition  with  Equa¬ 

tion  2. 9-8  eliminate  the  Lagrange  multipliers  and  yields 
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(2  9-11) 


Inclusion  of  the  inertial  effects  of  tail  masses,  etc.  is  now  obvious. 
Defining  the  total  body  mass  matrix  by 


(2.9-12) 


where  tails,  etc.  have  been  included,  the  panel  point  equations  defining 
fuselage  modes  are 
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The  transformation  separating  time  and  space  coordinates  is 


(2.9-13) 
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The  fuselage  modes  are  considered  to  be  normalized  on  the  total  mass  matrix 


[A}  ,  so  that 

(ii.))  m 

yw  -  w 

(©ip  (yf 

and 

(2.9-15) 

[A]U\  -  ^[K]>F]t 

(2.9-  16) 

where 

[01  - 

(2.9-17) 

define  the  mode  shapes  and  eigenvalues  for  the  fuselage. 


Wing  modes  are  defined  with  everything  constrained  except  elastic  wing 
motion.  Equations  2.0-5  for  the  wing  become 
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which  is  already  in  the  form  necessary  to  transform  to  modal 


W 

(2.9-18) 

coordinates. 


All  other  elements  are  handled  in  this  same  manner,  and  need  not  be 
written  explicitly.  The  transformation  to  modal  coordinates  is  then  quite 
simple;  substitutions  of  the  form  in  equations  2.9-1*+  arc  made  directly  into 
the  panel  point  equations.  It  is  the  definition  of  mode  shapes  for  the  main 
body  which  might  cause  some  debate;  the  choice  is  arbitrary,  but  it  is  felt 
that  the  above  definition  is  the  best  approximation  when  only  a  few  modes 
are  retained. 


In  the  case  of  a  gear  attached  to  a  wing,  the  wing  modes  are  not  de¬ 
fined  by  2.9-18.  Panel  point  displacement  of  the  gear  is  constrained  to 
motion  due  to  wing  motion,  constraints  are  written  defining  this  ge.°"  motion, 
and  a  term  such  as  [Twc3'[Ag]LTw6]  wil1  be  added  t0  £Aw]  in  the  final  equa¬ 
tion  defining  wing  mode  chapes. 
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The  transformation  to  modal  coordinates  for  a  problem  such  as  flutter 
analysis  is  very  useful  in  that  the  modal  coordinates  are  not  coupled  to¬ 
gether  in  the  final  equations,  and  the  simultaneous  equations  are  easily 
solved  on  a  digital  computer.  Were  the  problen  is  not  so  simple,  as  the 
rigid  body  motion  will  be  coupled  into  each  equation,  and  the  modal  co¬ 
ordinates  are  coupled  in  each  equation.  This  complexity  is  unavoidable. 

The  transformed  set  of  equations  are  considered  next.  The  total 
number  of  modal  coordinates  which  have  been  defined  are  not  generally 
used  due  to  the  resulting  large  number  of  variables  involved.  It  is 
desirable  to  truncate  the  square  matrix  [03  and  delete  the  column 
vector  to  the  desired  number  of  variables  for  each  component.  Then 
the  transformation 

-  L5>3 i  1  (2.9-19) 

is  made,  where  tfai  is  the  matrix  formed  by  deleting  the  higher  mode  shape 
columns,  and  is  not  square.  This  reduction  in  number  of  variables  may  be 
treated  in  one  of  two  ways.  The  set  of  equations  2.9-19  may  be  considered 
as  constraint  equations  on  the  system.  The  right-hand  side  of  the  panel 
point  equations  would  then  contain  total  panel  point  forces  minus  the  panel 
point  forces  which  constrain  the  higher  modes  of  deformation.  This  approach 
may  be  used  if  desired,  and  the  associated  Iagrange  multipliers  may  be  elim¬ 
inated  rather  easily.  A  somewhat  simpler  approach  may  be  used.  Direct  sub¬ 
stitution  of  equation  2.9-19  into  the  panel  point  equations  results  in  ob¬ 
taining  more  equations  than  unknowns,  and  the  solutions  are  not  unique. 
However,  a  linear  combination  of  these  equations  may  be  chosen  which  yields 
the  proper  result.  The  shorthand  notation  which  defines  this  operation  is 
the  preraultiplication  of  the  set  of  equations  by  the  transform  of  the  deleted 
modal  matrix,  [Wi  •  The  right-hand  side  of  the  equations  is  recognized 
as  the  contribution  of  the  panel  point  forces  to  the  modal  coordinates  re¬ 
tained,  which  is  exactly  the  result  obtained  by  the  first  approach. 

The  expansion  of  the  total  panel  point  displacement  for  the  i  -th 
component  is 


lP>i  *  MP);  +  SP'!;  Mai; 


(2.9-20) 


where  it  is  assumed  that  the  i  -th  component  is  attached  to  the  vehicle 
by  means  of  attaching  it  to  the  ^  -th  component.  As  it  is  possible  that  the 
y  rth  component  is  not  the  fuselage,  but  is  attached  to  the  fuselage,  then 

lp}j  =  [Tif]tPlr  Mp'Jj  *  Ulj 


4 1 


(2.9-21) 


An  example  is  that  of  a  gear  on  a  wing,  i  corresponding  to  the  gear  and  £ 
to  the  wing.  Using 


iP'h  -Wft 

expression  2.9-18  becomes 

<•  [TijUA'4 

+  U!; 


(2.5-22) 


(2.9-23) 


(2.9-24) 


One  may  well  appreciate  the  complexity  of  writing  the  equations  of  motion, 
2.8-3,  4,  5,  in  this  notation  for  a  general  body  whose  geometry  is  as  yet 
undefined.  Hence  the  general  sen  of  equations  using  modal  coordinates 
will  not  be  written;  rather,  the  transformation  will  be  performed  for  in¬ 
dividual  problems  as  necessary  to  reduce  the  number  of  degrees  of  freedom 
to  a  level  consistent  with  problem  requirements. 

2.10  POSITION  AND  VELOCITY  OF  A  BODY  POINT 


If  the  time  at  which  an  applied  force  is  imposed  on  the  body  is  to  be 
accurately  known,  as  well  as  its  magnitude,  which  may  depend  on  the  body 
kinematics,  it  is  necessary  that  the  position  of  the  point  to  which  it  is 
to  be  applied  is  also  known.  The  purpose  of  this  section  is  to  derive  that 
position  in  terms  of  the  variables  defining  the  motion  of  the  body.  This 
may  be  accomplished  by  the  integration  in  time  of  the  components  of  the 
velocity  of  the  point  relative  to  the  ground. 

Let  the  point  be  labeled  B.  Its  position  in  the  undeflected  body  relative 


to  the  origin  of  the  body  coordinate  system  is  ILa 
placement  relative  to  the  body  coordinate  system  is 
this  rxjint  may  be  expressed  by  either  of  the  forms 


V,  -  V„  i  *  V,,  Jl  <■  Vv  IK 
=  X,  1!  <-  Y„  Jl  *  2„  IK 


and  its  total  dis- 
IP„  •  The  velocity  of 


(2.10-1  ) 
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It  may  also  be  written  as 


VB  =  =  !R  +  SI*  Lg  +IFS  +  &*(PB 


(2.10-2) 


from  equation  2.6-7-  The  components  of  the  velocity  of  the  point,  expressed 
in  the  body  coordinate  system,  are  then 
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This  form  may  bt  expanded  into  components,  with  the  resulting  expressirn 
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As  the  point  B  may  be  in  a  component  other  than  the  fuselage,  the  term,  may 
be  bitten  in  the  component  coordinate  system.  Then 
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The  components  of  the  velocity  of  the  point  in  the  groj-.d  reference  system 
are  obtained  by  the  transfomation  of  this  expression. 
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(2.10-6^ 


The  position  of  the  point  in  the  ground  refe.ence  system  is  immediately 
obtained  as 
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X„  It)  =  X6(0)  +J*  XB  dt 
Yelt)  =  Y,W  Y0  <!t 
Z.ttY  =  zB(o\  +jf  Z,  it 


(2.10-7) 


where  X  e  (o'),  YB(o),  28(o)  are  the  coordinates  of  the  point  at  zero  time, 
the  beginning  of  the  problem.  They  may  be  found  from 
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(2.10-8) 


where  X(o),Y(o')  ,^.(o)  are  the  initial  components  of  the  position  vector 
to  the  body  coordinate  system  origin,  and  the  zero  subscript  on  the  Eulerian 
transformation  indicates  the  value  of  the  Euler  angles  at  zero  time. 
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SECTION  3 


APPLIED  FORCES 


3.1  GENERAL 

A  solution  of  the  equations  of  motion  written  in  Section  2.8  may  be 
obtained  when  the  applied  forces  have  been  defined.  As  the  number  of  equa¬ 
tions  is  the  same  as  the  number  of  variables  needed  to  define  the  motion,  the 
forces  must  be  completely  specified.  They  may  be  stated  as  known  constants, 
as  explicit  functions  of  time,  or  as  functional  relations  between  the  variables 
defining  the  motion. 

The  total  force  applied  to  the  vehicle  is  defined  by  the  vector  (Qx,Qy, 
Qy  ).  This  force  is  the  summation  of  all  the  externally  applied  forces  on 
the  vehicle,  expressed  in  the  body  coordinate  system.  The  contributions  to 
this  force  considered  in  Section  3*2  are 


Gravitational  Force 

-  © 

Thrust  Force 

-  Q 

Parachute  Force 

-  <Sl 

Aerodynamic  Force 

-  Q 

Ground  Force 

-  Q 

The  ground  forces  will  be  considered  for  a  variety  of  contacting  elements  on 
the  vehicle.  The  contribution  to  body  forces  is  developed  for  a  single  element 
of  each  type,  and  the  summation  over  the  elements  indicated.  The  total  exter¬ 
nal  moment  on  the  body  is  defined  by  the  vector  (Nx.N^.N’j.).  The  contribution 
from  each  of  the  above  forces  is  defined,  and  the  total  moment  obtained  by  a 
summation. 

The  definition  of  the  vectors  ({ Q=c'},{Q y/J  ,  {  Q  7/  )  )  is  dependent  both 
on  the  mathematical  model  to  be  used,  and  to  some  extent  the  manner  in  which 
the  problem  is  to  be  attacked.  The  panel  point  forces  in  each  set  of  panel 
point  equations  must  be  arranged  in  the  same  order  as  are  the  panel  point 
displacements.  They  may  include  the  distributed  forces  which  are  considered 
to  act  on  the  whole  vehicle,  the  interned  reactions  which  hold  the  vehicle 
together,  the  damping  forces  which  dissipate  energy  stored  in  elastic  defor¬ 
mations,  and  the  stroking  forces  between  components  which  move  as  rigid  bodies 
relative  to  one  another.  In  Section  3.3,  the  forces  on  the  component  rigid 
body  motions  are  defined  in  terms  of  ground  forces  and  stroking  forces.  In 
Section  3.4,  the  forces  on  the  panel  points  contributing  to  elastic  defor¬ 
mations  are  considered,  and  the  distinctions  made  as  to  which  of  these  are 
applied  forces  and  which  are  included  as  constraints  or  are  included  in  the 
definition  of  the  stiffness  matrix.  The  discussion  of  constraints  in  Section 
3.5  completes  the  definition  of  forces  in  the  system  of  equations. 

The  contents  of  this  section  are  intended  to  include  the  most  common 
forms  of  applied  forces  which  will  occur  in  landing  impact  problems.  M::*y  of 
them  will  hold  only  for  a  particular  configuration  and  are  not  meant  to  t; 
general  definitions  to  be  applied  for  any  problem.  The  individual  may  replace 
these  forms  with  others  more  suited  to  a  particular  application,  keeping  in 


mind  that  the  forces  must  be  completely  specified  either  in  terms  of  the 
variables  already  defined  or  as  explicit  functions  of  time. 

3.2  TOTAL  VEHICLE  APPLIED  FORCES 


In  this  section,,  the  forces  (Qx,  )  which  enter  the  rigid  body 
Equations  2.8-3  will  be  defined.  The  general  form  for  the  total  body  force 
is 


where  Q.  is  the  total  force,  given  by 


for  gravitational,  thrust,  ground,  parachute,  and  aerodynamic  forces, 
respectively. 

The  contributions  to  (Nx,Ny,,N^.)  entering  the  Equations  2.8-4  will  be 
calculated  in  each  section;  the  forms  vary  considerably  depending  on  the  type 
of  force.  The  contributions  from  each  section  may  then  be  summed  to  produce 
the  total  moment. 

3.2.1  Gravitational  Force,  Qw 

The  gravitational  force  acting  on  the  vehicle  is  directed  along  the 
positive  Z  axis;  the  magnitude  is  the  product  of  the  total  vehicle  mass 
and  the  local  acceleration  due  to  gravity  (  ^  ): 


Q„  =  MjlK 


(5-2-1) 


fhe  components  of  the  weight  vector  in  the  body  coordinate  system  are 
given  in  terms  of  the  Eulerian  transformation  as 


f~  SIN  9  ^ 

=  (  cos  0  sin  i  | 

jeos  0  cos 


(3.2-2) 
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Since  the  origin  of  the  tody  coordinate  system  will  not  be  coincident 
with  the  instantaneous  position  of  the  center  of  mass,  a  moment  about  the 
body  coordinate  axes  will  result  from  the  gravitational  force.  Let  |pw  be 
the  distance  from  the  origin  of  the  body  coordinate  system  to  the  center  of 
mass  of  the  total  vehicle.  The  moment  is  then  given  by 


fNjl  fipw  x  eu  *  a! 
<N*>  =  <  fPw 

In J  I Pw * * kj 

Vs.  ^  J 


(5.2-5) 


The  vector  Pw  is  determined  in  terms  of  the  displacement  from  initial  position 
of  all  the  vehicle  relative  to  the  body  coordinate  system; 


ip  _l 


(5.2-4) 


The  integral  may  be  broken  into  integrals  over  each  of  the  N  components  of 
the  vehicle; 


(5.2-5) 


In  terms  of  the  matrix  notation  which  arises  from  the  panel  point  concept 
and  the  interpolation  procedure,  this  form  is  rewritten.  The  components  of 
the  vector  Fw  are  then 


whereby] is  the  transformation  between  the  body  coordinate  system  and  the 
i  -th  component  coordinate  system.  The  moment  is  then  obtained  by  substi¬ 
tuting  Equation  3.2-6  into  Equation  3*2-3. 
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The  moment  calculated  here  Is  small  compared  to  moments  from  the  ground  . 
forces  and  will  normally  he  neglected  in  landing  impact  problems. 

3«2.2  Thrust  Force.  &  r 

It  is  assumed  ir  the  discussion  which  follows  that  the  magnitude  of  the 
thrust  force  is  known  explicitly  as  &  function  of  time,  or  is  a  krown  con¬ 
stant,  or  that  it  is  known  as  a  function  of  the  variables  defining  the  vehicle 
motion.  The  line  of  action  of  the  thrust  vector  is  defined  by  the  positions 
of  two  points  on  the  line,  labeled  (a,b).  The  positions  of  these  points  are 
defined  by  (L+P  )a  and  (  1L  +  P  )(,  as  in  Fig.  3.  The  thrust  vector  then  is 
collinear  with  the  vector  difference  of  these  two  quantities, 

U-T  =  (IL-IP)^  -(IL*lP)k 


which  may  be  expanded  to 


+(P»  -Px„  )i+ (P».  -P,k )  Jl ' -  (p,.-  P,„ )  IK 


The  direction  cosines  of  the  vector  1LT  in  the  body  coordinate  system  are  then 
variables,  given  by 


f(xft-*b)+(Pxa  “  Pxb  )"! 


(3-2-9) 


The  components  of  the  thrust  in  the  body  coordinate  system  are  then 


(3.2-10) 


J»8 


INSTANTANEOUS 
LINE  OF  THRUST 


Figure  3»  Line  of  Action  of  the  Thrust  Vector 


where  Qyis  the  thrust  magnitude.  If  there  are  several  independent  thrust 
sources,  a  sum  of  terms  like  this  one  will  give  the  total  force;  if  the  thrust 
magnitudes  are  the  same, 


cQ  as.  fV 

Q.)  =Qr  v 

a;jT  W 


(3.2-11) 


where  JC  runs  over  the  various  thrust  sources.  The  moment  due  to  any  one  of 
these  is 


fN.l  f(t+  ip),  x  <aT  •  il 
In?  j  k* ip),  x aT  •  iid 


(3-2-12) 


0  -C^-a+P/ja  )  (fa+Pip.)  fox 

=  (r*+p,J  o  -(ia+p£)  <q* 


~(fcL  +  Pyt  )  (Xa+PX0J  0 


2.2.3  Parachute  Force, 

Parachutes  are  used  in  conventional  aircraft  to  shorten  the  ground  run. 
For  vertical  alightment  vehicles,  the  parachute  may  he  the  primary  lift  device, 
or  it  may  only  be  used  to  ensure  that  vehicle  orientation  is  essentially 
vertical.  The  general  expression  for  the  parachute  force  is 


=  C'o  yzr.  \  ~~  z/°  SCD^r I  Vp 


(3.2-13) 
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where  the  various  parameters  are 


'  atmospheric  density 
C0S  effective  drag  area  of  the  parachute 

I!locity  of  the  parachute  attachment  point  with  respect  to 
the  atmosphere 


In  Section  2.10,  the  velocity  of  a  point  in  the  body  is  derived  in  detail.  If 
the  components  of  the  wind  velocity  expressed  in  the  body  coordinate  system 
are  given  as  vvx^ then  from  Eq.  2.10-5,  the  velocity  of  the 
parachute  attachment  point  relative  to  the  atmosphere  expressed  in  component 
form  in  the  body  coordinate  system  is 


In  terms  of  these  relative  velocity  components  and  the  general  parachute  force 
expression,  the  parachute  force  is 


where  the  positive  sign  of  the  radical  is  understood. 

If  the  parachute  is  to  open  during  the  time  period  in  which  the  landing 
is  analyzed,  the  drag  area  C0S  will  be  a  function  of  time.  Thus,  if  the 
parachute  is  initiated  at  time  t  j  ,  and  it  takes  tj  seconds  to  deploy  and 
seconds  to  open,  the  drag  area  is 
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0 

CDS  =<  CoS (t) 


t  <ti  +  td 


(CoS^ 


ti  +  td  ti  +  td+t0 

t£+td  +  tO  <  1 


(3.2-16) 


where  CoS(t)  is  the  growth  function  of  the  drag  area  during  the  opening 
period,  and  (CdS)*  is  its  final  value.  For  conventional  parachutes,  these 
parameters  may  be  found  in  reference  6 • 


The  moment  about  the  axes  of  the  body  coordinate  system  is  written  in 
component  form  as 


(ClL  +IP)p  x 
\(lL  +  7P)  x 


0 

(7r  +  P^p  ) 

~(#r 


“(7p  +  P'Jp) 

0 

Cxp  +  P*?) 


_(XF+P*p 

0 


(3.2-17) 


3.2.4  Aerodynamic  Force, 

The  aerodynamic  forces  acting  on  the  vehicle  are  those  forces  exerted  by 
the  surrounding  atmosphere  resisting  the  motion  of  the  vehicle  .  These  forces 
are  defined  here  in  terms  of  the  notation  used  in  describing  conventional  air¬ 
craft.  The  contribution  to  these  forces  due  to  clastic  motion  of  vehicle 
components  will  be  discussed  in  Section  3*4.  The  contribution  involving  rigid 
body  dynamics  and  control  surface  deflections  will  be  considered  here. 

The  general  forms  for  the  components  of  the  forces  and  moments  are 


=  J_  /oStzt' 


(3.2-18) 


52 


(3.2-19) 


where  the  indicated  parameters  are  defined  by 


/°  -  atmospheric  density 
S  =  wing  area 
c  =  mean  aerodynamic  chord 
Hr  -  wing  span 

HT  =  velocity  of  the  vehicle  relative  to  the  atmosphere 


c 

c 

c. 


non-dimensional  aerodynamic 
*  force  coefficients  referred 
to  the  body  axis 


Ct  1  non-dimensional  aerodynamic 
)  =  moment  coefficients  referred 
Q^j  to  the  body  axis 


Each  of  the  aerodynamic  coefficients  is  a  function  of  the  body  orientation 
and  velocities,  control  surface  deflections  and  velocities,  and  elastic  body 
motion.  For  the  vehicle  without;  elastic  moticis. 


(3-2-20) 


whereof,  the  angle  of  attack,  and/? ,  the  sideslip  angle,  define  the  orien¬ 
tation  of  the  free  stream  velocity  E/"  with  respect  to  the  body  x-axis.  The 
control  surface  deflection, 8 ,  and  the  corresponding  velocity,  j  ,  are 
considered  here  for  an  arbitrary  control  surface.  The  subscript  t  refers  to 
arty  of  the  above-mentioned  coefficients. 

In  the  theory  of  aerodynamics,  it  is  usually  adequate  to  assume  that  the 
aerodynamic  coefficients  are  linear  in  each  of  the  variables.  This  is  equiv¬ 
alent  to  expanding  the  coefficient  in  a  Taylor's  series  and  omitting  all 
non-linear  terms; 


Ci*Ct(0)  +  |fe(0)^+^(0)/ 


dCj 


( 3-2-21) 


The  higher  order  terms  are  assumed  negligible  ir.  the  linear  theory. 
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For  landing  response  problems,  not  all  of  the  linear  terms  are 
appreciable  in  a  given  coefficient.  In  the  following  list  of  coefficients, 
only  the  most  important  terms  have  been  retained. 


r< icx/9°c\ '  r«o 

cx  -  C3C()  +  <  <  $e  > 

[zCjiStj  L^J 

0 


^  C*  / 

ac;/^Jo 


(3.2-24) 


(3.2-23) 


fac~  /^vr 

+<  ac^/c)^EW 

(_d  C/m.  /  1 


(3-2-26) 


(3-2-27) 


C'm.  ~ 


j3  Cm.  I  R 


The  indicated  variables  are 

= pitch  control  surface  deflection,  such  as  elevator,  ailevator 
=flap  deflecton 
=  rudder  deflection 

$  =roll  control  surface  deflection,  such  as  aileron  or  spoiler 

The  values  C*e,  C^o>C*pt  „  for  zero  deflections  and  the  partial  derivatives 
evaluated  at  zero  deflections  are  obtained  from  wind  tunnel  tests. 

The  above  terms  are  obviously  not  all-inclusive,  but  should  provide  an 
adequate  definition  for  most  landing  problems.  A  complete  treatment  of 
aerodynamic  forces  may  be  found  in  reference  7  •  The  control  surface 
deflections  are  assumed  to  be  known  as  constants  or  explicit  functions  of  time. 

The  only  other  matter  to  be  covered  in  this  Section  is  the  calculation  of 
the  angles  of  attack  and  sideslip.  The  velocity  which  enters  Equations  3.2-18, 
19,  which  define  the  force  and  moment  on  the  body,  is  the  difference  between 
the  body  translational  velocity  and  wind  velocity,  which  yields 

2  w*  jv,  -rfO 

—  -/^r‘4W  /  \  /V~ u  >  (5.2-28) 

15*  -^yj  U?  -*VJ 

The  angle  of  attack,  oC,  is  defined  as  the  angle  between  the  velocity  vector 
jf  projected  into  the  x-  plane  anv  tne  x-axis.  This  requires 


1  ^2  ~ 

/£TX  -/ZCun 


(5.2-29) 


Similarly  the  sideslip  angle  is 
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(5-2-30) 


/3  -  TAN 


/(fy  —AfyYJ 


For  conventional  aircraft  landings,  the  forward  velocity ,/lT-x,  is  much 
larger  than  the  drift  or  sinking  speeds,  and  the  following  relations  are 
very  nearly  true. 


C'K  --n. ) 2 


(5.2-51) 


a C 


/^x-'^xtv 


(5.2-52) 


/&x~ 


(5.2-55) 


3.2.5  Ground  Forces ,  Qe 

In  this  Section,  forces  on  the  tody  due  to  interaction  with  the  ground 
will  he  developed.  The  section  will  he  divided  according  to  the  type  of 
contacting  element  under  coosidertfcion.  These  elements  will  he  of  se-veral 
types:  surface  pods,  tires,  skis  and  skids,  spikes,  and  gas -filled  hags. 

In  some  areas,  the  ground  will  he  considered  to  he  rigid.  In  these 
cases,  the  contacting  element  is  considered  to  develop  a  coefficient  of 
frlrt-’cn  with  respect  to  the  ground,  which  may  he  dependent  on  parameters 
of  the  motion.  In  other  areas,  the  ground  is  allowed  flexibility,  viscosity, 
or  compressibility  to  varying  extents. 

The  ground  coordinate  system  is  used  in  describing  most  of  these  forces. 
The  components  in  the  ground  coordinate  system  of  the  force  exerted  on  the 
body  by  the  ground,  called  ground  reactions,  are  (D^Sg.Vg  )m.  ,  where  the 
subscript  indicates  a  particular  contacting  element  on  the  vehicle.  The 
total  force  on  the  body  due  to  the  ground  is  then 


The  separate  contributions  to  force  on  the  body  are 


!Q, 


D, 


9  -Cr]s: 


a 


(3.2-35) 
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/H, 


Moments  about  the  axes  of  the  body  coordinate  system  are  found  from  the 
general  form 


(3.2-36) 


where  (IL^+IP^)  is  the  instantaneous,  lever  arm  of  the  force  Q*.  This  is 
written  in  component  form  as 


(tr-PJGu  --(?+ R>) 

(?+p’)q„  -o*p,;<2 
(*+p')q,;-(}.-p3)q 


0* 

xc. 


/ft 
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where  the  forces  and  lever  arms  for  each  contacting  element  are  indicated  in 
component  form. 

The  forces  as  written  in  the  following  pages  will  be  written  for  a 
single  contacting  element,  with  no  subscript  attached  to  identify  the  element. 
It  will  then  be  understood  that  these  forms  must  be  written  for  each  contact¬ 
ing  element  when  making  use  of  the  Equations  3.2-3U,  37  for  body  forces  and 
moments . 


A  surface  pad  is  a  mass  rigidly  attached  to  the  end  of  the  gear  piston. 
It  is  not  allowed  a  degree  of  freedom  for  motion,  hut  simply  moves  with  the 
piston.  It  has  sufficient  contacting  area  so  that  its  ground  penetration  is 
small  if  the  ground  is  considered  to  be  soft. 

The  ground  force  will  first  be  described  for  the  case  when  the  ground  is 
rigid.  The  surface  pad  will  develop  a  coefficient  of  friction  with  respect  to 
the  ground.  The  ground  force  will  be  developed  entirely  in  terms  of  the  body 
motion.  The  force  from  the  ground  along  the  stroking  axis  is  first  developed. 

The  drag  force  on  the  pad  is  opposite  to  the  direction  of  the  horizontal 
pad  velocity,  which  requires 


(3.2-38) 


•  « 

where  the  components  of  the  pad  velocity  in  the  ground  plane  are  (Xp,Y p  ),  the 
force  on  the  pad  parallel  to  theX-axis  isDg,  and  that  along  the  y  -axis  is  Sv 
The  magnitude  of  the  drag  force  is  the  product  of  vertical  force  on  the  ground 
with  the  coefficient  of  friction,  yielding 


(3.2-59) 


(3.2-40) 


where  whe  forms  in  velocity  components  divided  by  their  magnitudes  resolve  the 
arbitrariness  in  sign  in  solving  the  quadratic  Equation  3*2-39.  The  forces  on 
the  ground  are  related  to  the  forces  on  the  pad  by 
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(3.2-41) 


=  [r]'MP 


which  may  be  inverted  to  give 


(3-2-42) 


The  force Q^cp is  along  the  stroking  axis  of  the  gear.  After  contact,  there 
is  no  motion  of  the  piston  relative  to  the  ground  along  this  axis.  Then 
is  the  negative  of  the  piston  stroking  force,  Q>'s’  which  is  positive  downward 
along  the  stroking  axis.  The  stroking  force  is  the  sum  of  the  internal  gear 
axial  forces  to  be  developed  later.  The  Equations  3.2-40,  42  combine  with  the 
above  definition  to  relate  the  vertical  force  on  the  ground  to  the  piston  strok¬ 
ing  force: 


Va  = 


Q, 


0?  , 

o  JM  [r]|*(vp/l  *P00+  (*P/*f)T/z 


(3.2-43 


Combining  Equations  3.2-40,  43  yields  the  ground  force  for  a  single  contacting 
element, 


59 


,oV  ,  (>(xP/ixri)C !+(tf/xP)2]'^’ 
0  >[v-][r]^(Yp/|tpl)[l  +  (XP/YP)2]'/z 


V-  J 


From  this  expression  and  the  Equations  3*2-3^>  37  the  body  forces  and  moments 
may  readily  be  calculated. 

This  form  for  ground  forces  is  quite  complicated;  it  may,  however,  be 
simplified  for  most  applications.  One  of  the  examples  in  Section  4  will  make 
use  of  this  form,  and  the  simplifications  used  there  will  perhaps  enlighten 
the  reader. 

If  the  ground  is  soft,  the  above  statements  will  not  be  applicable.  It 
is  suggested  that  the  following  forms  for  ground  force  be  used.  The  vertical 
ground  force  may  be  represented  by  a  one-way  spring,  that  is. 


i?  >  0 
Zp,Mo 


(3.2-U5) 


where  the  velocity  of  the  surface  pad  normal  to  the  ground  is  obtained  from 
Equation  2.10-6.  The  spring  ratej^,  must  be  obtained  from  experiment,  as  it 
will  depend  on  the  geometry  of  the  pad  and  characteristics  of  the  ground. 

The  penetration  depth  d  is  obtained  by  integration  of  Zp  following  touchdown. 

The  ground  forces  in  the  ground  plane  may  be  represented  by 


D& 


d 


d>0 

dfCO 
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s 


G 


d^O 


where  the  velocities  are  obtained  from  Equation  2.10-6  and  the  constants*  2 
must  be  determined  from  experiment.  The  experimental  problems  might  be  avoided 
by  use  of  some  analytical  formulas  in  terms  of  pad  area  and  ground  character¬ 
istics;  a  form  of  this  sort  may  be  useful  in  prediction  of  peak  loads,  but  must 
be  used  with  caution  in  investigating  stability. 

3. 2. 5. 2  Tires 


Forces  generated  in  the  system  due  to  interaction  of  a  tire  and  the  ground 
have  never  been  very  accurately  described  on  a  theoretical  basis.  Generally, 
experimental  data  furnished  by  the  manufacturer  is  more  reliable  than  purely 
theoretical  data,  and  such  should  be  used  when  available.  These  data  may  in¬ 
clude  vertical  load-deflection  characteristics,  coefficient  of  sliding  friction, 
cornering  coefficient,  self-aligning  torque  coefficient,  or  other.  The  forms 
presented  here  are  in  general  taken  from  the  available  literature  and  are  felt 
to  be  consistent  with  the  state  of  the  art  concerning  tire  characteristics. 

The  vertical  load-deflection  characteristics  of  a  tire  will  be  discussed 
first.  The  forces  in  the  ground  plane  will  then  be  discussed  for  the  periods 
in  time  before  and  after  wheel  spinup. 

The  ground  is  assumed  to  be  rigid  throughout  this  development. 

Vertical  Ground  Force,  VG 

The  force  on  the  body  normal  to  the  ground  plane  is  written  here  in 
terms  of  physical  parameters  of  the  tire,  and  should  be  used  only  if  data 
from  an  experimental  program  are  not  available.  The  variation  in  ground  force 
with  tire  deflection  has  been  writ.rea  in  many  ways,  each  formula  representing 
an  approximate  fit  to  the  observed  data.  Hadekel  (Ref.  8  )  notes  tha-  the 
most  rational  approach  to  vertical  load-deflection  characteristics  produces  the 
form 

f  0  £<0 

Vo  =<  _ A  (P  +  Pc  +ap)  0<  S  <Sb  (5'2J,8> 

where  the  tire  contacting  area,  A  ,  is  given  by 
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and  the  various  other  quantities  are 


&  -  tire  vertical  deflection 
or  =  tire  width 
P  = undeflected  tire  pressure 
fit  = tire  wall  equivalent  pressure 
SP=  pressure  rise  on  deflection 
T5  =tire  outer  diameter 

=tire  deflection  on  initial  bottoming 

The  change  in  pressure  with  deflection  is  found  in  the  static  case  to  be 


This  approach  yields  results  negligibly  different  from  dynamic  loads  obtained 
from  drop  tests  for  low  deflections;  the  dynamic  loads  rise  more  rapidly  for 
large  deflections.  This  may  be  corrected  to  some  extent  by  using  a  polytropic 
compression  form  such  as  an  adiabatic  compression  would  yield.  The  above  forms 
are  inadequate  for  tire  bottoming  in  any  case.  It  is  then  necessary  to  replace 
Equation  3*2-48  by  the  form 

Vg  =  VGb  Sb)  (5'2'52) 

=  (Vob+4b  £b) -Jkb  8  £> 
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where  N^l  is  the  vertical  ground  force  vht  .  the  tire  first  hoc  terns  and  Jbt 
is  the  tire  bottcaing  spring  rate. 

There  sire  several  points  which  should  be  considered  here.  Hysteresis 
effects,  which  may  be  observed  in  standing  tires,  are  considered  by  Hadekel  to 
be  negligible  in  dynamics  problems  with  high  spinup  speeds;  hence  cquivaJ ent 
damping  in  tires  will  not  be  considered.  Tire  lateral  deflections  produce 
vertical  deflections  for  constant  vertical  load,  and  so  both  should  be  specified 
in  order  to  find  vertical  load.  The  omission  of  this  effect  will  not  apprecia¬ 
bly  change  the  time  history  of  the  load,  except  that  the  time  for  tire  bottom¬ 
ing  may  be  in  slight  error. 

It  has  been  assumed  here  that  the  vertical  tire  deflection  in  the  wheel 
plane  and  the  "average"  deflection  observed  in  a  tire  with  the  wheel  plane 
slightly  tilted  from  the  vertical  are  identical.  Thun,  the  load  deflection 
characteristics  are  assumed  independent  of  wheel  plane  orientation  for  small 
angles  from  the  vertical. 

The  vertical  deflection  is  found  as  follows.  The  distance  to  the  axle 
from  the  ground,  which  will  be  negative  according  to  the  sign  convention,  is 


ZA  -  (iR  +  ILa  +iPA  j  *  IfcL  (5.2-53) 

where  are  the  undeflectcd  position  and  the  deflection  of  the  axle.  This 

may  be  written  in  component  form  as 


where  1  is  the  position  of  the  origin  of  the  body  coordinate  system,  and  the 
subscript  refers  to  the  axle  and  tne  coordinate  system  of  the  gear  to  >hich  it 
belongs .  The  inste  ntaneous  angle  which  the  axle  makes  with  the  ground,  a, 
defines  the  projection  of  the  undeflected  tire  radius  on  the  normal  to  the 
ground,  which  yields  the  deflection  as 


S  =  H 


0 

M  [rlUo  \  +  r  cos  ocf 


(5-2-55) 


63 


For  normal  aircraft  landings,  the  angle is  small,  so  thatCOS^^is  very 
nearly  unity.  The  error  here  is  of  the  same  order  as  that  in  using  the  deflec¬ 
tion  obtained  in  the  formula  presented  by  Hadekel. 

For  most  aircraft  landings,  the  Euler  angles  remain  small,  and  the  Euler ian 
transformation  may  be  linearized.  The  resulting  expression  for  tire  deflection 
is  then 


fx'+p^V  ,f-e 

S  =  H+<*-'  +  p,.  >0]  <  * 

l 1 


+  r 
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The  vertical  ground  force  as  written  here  is  valid  both  before  and  after  wheel 
spinup.  The  forces  on  the  ground  in  the  ground  plane,  however,  are  dependent 
on  whether  or  not  the  wheel  is  spun  up.  In  these  areas,  the  load-deflection 
characteristics  lateral  and  tangential  to  the  tire  will  not  be  used  in  the  sense 
of  Equation  3.2-48.  Fre-spinup  forces  are  derived  from  the  sliding  coefficient 
of  friction,  the  vertical  force,  and  the  geometry.  Post-spinup  force  along  the 
wheel  plane  is  considered  negligible,  and  that  normal  to  the  wheel  plane  is  found 
from  the  cornering  characteristics  of  the  tire. 

Ground  Forces  in  the  Ground  Plane,  D& .  && 

Forces  from  the  ground  along  the  X  -axis  and  Y  -axis,  D^and  Sg  ,  are 
usually  defined  in  two  separate  regions:  pre-spinup  and  post-spinup.  In  the 
region  of  transition  between  the  two,  a  form  may  be  picked  which  is  intuitive¬ 
ly  satisfying  in  order  to  obtain  continuous  forcing  functions.  A  search  of  the 
literature  shows  that  no  analytical  procedures  derivable  from  physical  laws 
are  in  existence  for  this  region,  and  no  experimental  data  is  of  sufficient 
accuracy  to  define  any  variation  with  tire  parameters. 

It  is  also  convenient  to  present  in  this  section  the  moment  about  the 
axle  during  spinup. 


Pre-Spinup  forces  in  the  Ground  Plane 


The  forms  derived  here  are  based  on  the  assumption  that  the  force  from 
the  ground,  in  the  ground  plane,  is  in  the  opposite  direction  of  the  axle 
velc.ity  vector  parallel  to  the  ground,'^,  given  by 


x  =  xA  u  +  ya  j 


(3-2-57) 


The  magnitude  of  the  force  is  assumed  to  be  the  product  of  vertical  force  with 
the  sliding  coefficient  of  friction  between  the  tire  and  the  ground.  The 
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components  of  the  force,  as  pictured  in  Fig.  h,  are  then 


t  tt  (5.2-58) 


t  <  H  (5-2-59) 


where t-^  is  the  spinup  time.  The  components  of  the  axle  velocity  may  be  found 
from  Station  2.10-6.  For  conventional  aircraft  landings,  the  forward  velocity 
is  larger  than  any  other,  and  the  approximation  |VaMXaI  is  valid.  Other¬ 
wise,  the  root-sum-square  representing  the  velocity  magnitude  must  be  used. 

Some  discussion  of  the  variation  of  the  sliding  coefficient  of  friction 
w.'th  velocity,  pressure,  and  temperature  is  worthwhile  here.  Hample  (Ref.  9) 
presents  experimental  data  for  pressure,  load,  and  temperature  variations  of 
the  static  coefficient  of  friction  for  tire  materials.  Velocity  variations 
are  presented  by  Luthman  (Ref .10)  and  Gough  (P?f .  ll),  et  al.  The  data  by 
Hample  show  that  the  static  coefficient  decreases  sharply  from  room  temper¬ 
ature  to  300°  F,  then  gradually  to  500°  F,  and  again  sharply  to  the  melting 
point.  Material  was  taken  from  a  B-29  nose  wheel  tire.  He  also  notes  that 
the  static  coefficient  decreases  with  increasing  normal  pressure,  the  decrease 
becoming  sharper  with  increasing  temperature.  This  variation  is  substantiated 
by  Luthman,  although  his  main  interest  is  the  variation  with  velocity.  This 
is  characterized  by  a  general  decrease  as  velocity  increases,  but  with  oscilla¬ 
tions  superimposed.  Luthman  states  that  this  phenomenon  has  been  noted  prev¬ 
iously,  but  that  no  explanation  has  been  found.  Gough,  et  al.,  present  similar 
although  less  extensive  data  points;  the  oscillations  seem  jess  dominant,  and 
they  present  a  smooth  curve  through  the  data  points. 

These  statements  describe  the  variations  in  the  coefficient  of  friction 
in  a  general  manner.  Theoretical  work  on  the  subject  is  generally  avoided  due 
to  the  many  coupled  variables  involved.  Since  the  relative  velocity  of  the  tire 
footprint  and  the  ground  varies  from  air  speed  to  zero  during  spinup,  it  is  ad¬ 
visable  <,0  include  the  variation  of  the  coefficient  with  velocity.  A  linear 
term  in  contact  pressure  may  also  be  included.  Denoting  the  relative  ve’ccity 
by/#^  and  the  contact  pressure  by  PCT  ,  the  variation  may  be  approximated  by 


SG  =/^s  Vg 


Ya 

IVaI 
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(3.2-60) 


q^R 

where  B  is  the  value  for  low  pressure  and  high  velocity,  C  is  the  ratio  of 
values  at  high  to  low  velocities,  and Q.  W  must  be  determined  from  experiment. 
Contact  pressure  as  a  function  of  vertical  deflection  is  derived  from  Equation 
3-2-48,  the  ratio  of  vertical  load  to  footprint  area.  To  the  extent  to  which 
the  variation  in  coefficient  of  friction  with  velocity  is  known,  the  relative 
velocity  may  be  approximated  by 


•  P„  )  =  B(  I  -A  P„/l  +  f; 


AT*.  *  |VAI  -  r  7C 


(3.2-61) 


It  will  be  noted  later  that  these  forms  may  be  too  complicated  for 
ordinary  use,  but  simplified  forms  can  be  used. 


The  above  equations  are  all  valid  for  t  <  tfc ,  where 
time,  is  found  from  solution  of  the  equations  of  motion, 
time  may  be  approximated  with  use  of  a  formula  by  Flugge 


,  the  spinup 
If  desired,  spinup 
(Ref.  12): 


,  .  I  Ill/lTx.  2l~ 

1  ~  r  J  Ap  A  (3.2-62) 

where  I  is  the  wheel  moment  of  inertia,  is  the  gear  forward  spring  rate,  and 
F  is  the  vertical  load  time  derivative  (assumed  constant). 

It  is  useful  to  insert  at  this  point  the  forces  along  and  normal  to  the 
line  of  intersection  of  the  wheel  plane  and  the  ground  plane,  although  they  are 
not  used  in  determining  body  forces. 

The  force  along  the  line  is  the  force  used  in  defining  the  spinup  moment, 
and  is  designated  Wsu .  The  force  at  right  angles  is  the  side  force,  designated 
Vv's  .  These  are  forces  on  the  tire,  given  by 


W£„ 

=  VQ  cos  a 

t  <  tt 

(3.2-63) 

Ws 

-  Vg  sin  /) 

t  <  it 

(3.2-64) 
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where  the  angle  A  ,  shown  in  Pig.  4,  is  given  by 

(fy)  0.2-< 

The  angle/#  is  that  through  which  the  axle  is  rotated  about  the  gear  center- 
line  due  to  torsional  elasticity,  and  ‘p’  is  defined  in  the  Eulerian  transfor¬ 
mation.  The  spinup  moment  is  given  by 


Nsu  *(T-S)W„ 

•  • 

Note  that  for  conventional  airplane  landings,  Y  «X  ,  hence 


-i 

SIN 


(3.2-67) 


Post-Spinup  Forces  in  the  Ground  Plane 

The  force  from  the  ground  in  the  ground  plane  following  spinup  is  usually 
considered  to  be  normal  to  the  line  of  intersection  of  ground  plane  and  wheel 
plane.  This  assumes  rolling  friction  is  negligible  and  wheel  braking  is  not 
present.  The  side  force  for  conventional  aircraft  landings  is 


Ws  =  -  K  A  t  >  tt  (3.2-68) 

where  K  is  the  tire  cornering  coefficient  and  A  is  the  slip  angle,  defined  in 
Equations  3.2-65,  67.  If  the  slip  angle  and  vertical  force  are  small,  this 
form  is  valid  with  the  cornering  coefficient  constant.  For  high  vertical  loads, 
a  method  (Ref.  13)  has  been  devised  to  yield  a  cornering  coefficient  dependent 
on  tiro  vertical  deflection, 


v  _  if  (X  PCot  +  Cj) 

G 


(3.2-69) 
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vhert 


(5-2-70) 


The  remaining  various  parameter  s  are 

Q  =  tire  section  radius 
P  =  inflation  pressure 
d  =  head  seat  radius 
7  =  tire  undeflected  radius 

G  =  r  ~d 

£  s  tire  vertical  deflection 

This  approach  is  valid  for  high  vertical  loads  prior  to  tire  bottoming,  but  will 
still  be  restricted  to  small  slip  angles.  As  the  slip  angle  increases,  the  side 
load  reaches  a  maximum,  then  decreases  'co  the  value  which  would  occur  if  the  tire 
were  skidding  laterally.  Kadekel  presents  some  data  on  this  effect;  no  theoret¬ 
ical  forms  for  high  slip  angles  have  been  developed. 

Forces  along  the  ground  coordinate  axes  for  the  case  where  I  are  then 
given  by 


D.  =WS  (A-Y#/XJ  t>tt 

(5.2-71) 

s5  =  ws  t  >  t, 

(5-2-72) 

Ref.  12  also  gives  the  tire  self -aligning  torque  M 
angle  A  as 

arising  from  slip 

M  =  /m 

(5.2-75) 

in  the  range  <  5  ,  where 

nm,  =  12  a  (l.57  G  - a)(ajr  +G)p 

(5.2-71*) 
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3. 2. 5. 3  Spike-Soil  forces 


Penetration  on  impact  is  a  problem  not  easily  handled,  and  is  not  discussed 
in  the  literature  which  is  readily  available.  However,  an  approximate  form  may 
be  developed  which  is  useful  for  impact  on  sand  or  hard  soils. 

Penetration  forces  may  be  broken  into  two  groups:  compression  forces  and 
friction  forces.  Compression  forces  are  described  as  follows. 

The  bulk  modulus  of  a  medium  is  defined  as  the  change  in  pressure  per  unit 
volume  on  compression  of  the  medium: 


B  =  V 


dP 

dv 


(5-2-75) 


This  definition  is  valid  only  for  static  pressures  and  confined  volumes.  Ref. 

14  notes  that  dynamic  values  of  the  bulk  modulus  of  sandy  soils  are  200  to 
300  percent  of  the  static  value.  The  static  form  will  be  used  to  formulate 
the  compression  forces  on  impacting;  the  volume,  volume  change,  and  pressure 
change  will  be  considered  time  dependent,  and  the  value  of  the  dynamic  bulk 
modulus  will  be  used  for  calculations. 

Consider  a  body  with  velocity  /V  and  cross-sectional  area  A  normal  to  the 
velocity  vector.  As  this  area  moves  through  the  medium  an  amount  ds  ,  a  com¬ 
pression  wave  travels  outward  radially  from  the  area.  The  volume  encompassed 
by  the  wave  is 


V  = 


3 


(5.2-76) 


where  £  is  the  distance  from  the  area  to  the  wave  front.  The  pressure  change 
due  to  expansion  of  the  wave  front  an  amount  d  $  is 

dP  =  Fd(-^)=FayF)  =  -F^  <«.„ 


so  that 


VdP  =  --§■  FdJ 


(5  2-78) 
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This  form  is  independent  of  the  radius  of  the  compression  vave.  The  radius  may 
then  be  considered  small,  and  the  force  F  as  that  force  on  the  soil  produced 
by  the  notion  of  the  area  A  .  The  change  in  volume  of  the  soil  is  related  to 
the  cross-sectional  area  A  by  the  penetration  dS  vhich  occurs  in  the  time  for 
the  pressure  change  dP  to  occur: 


dV  =  -AdS 


(3.2-79) 


The  change  in  pressure  vave  radius  may  be  written 


(3-2-30) 


vhere  V  is  the  velocity  of  the  compression  vave  in  the  medium,  since  S  may 
be  considered  implicitly  dependent  on  time.  The  bulk  modulus  is  then 


o_  2.  _F  V 

D"  3  A  ^  (3-2-81) 


vhere  the  force  F  is  now  that  on  the  area  A  ,  and  not  on  the  soil,  since  a  minus 
sign  has  been  added.  The  force  opposing  the  velocity  /2T  is  then 


This  form  will  be  used  to  calculate  forces  and  moments  on  a  spike.  In  the 
derivation  the  assurption  has  been  made  that  E>  and  V"  are  independent  of  pene¬ 
tration  depth,  and  that  no  local  contracting  will  occur.  The  omission  of  these 
factors  will  not  affect  the  portion  of  a  landing  during  vhich  peak  loads  occur. 
However,  the  stability  of  a  body  landing  on  a  single  spike  is  determined  from  re¬ 
straining  moments  on  the  spike  alone,  and  this  area  will  not  be  described  too 
accurately. 

Spike  forces  will  be  calculated  for  a  single  spike,  symmetrically  located 
in  a  symmetrical  body,  from  the  form  in  Equation  3.2-82.  For  the  spike  problem 
to  be  formulated,  the  body  5-axis  will  remain  in  a  plane  normal  to  the  ground. 
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Figure  5«  Spike  Geometry 
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and  the  body  nay  rotate  about  the  M/  -axis  only.  There  will  then  be  two  com¬ 
ponents  of  force  on  the  body  and  one  component  of  moment.  The  extension  from 
this  problem  to  one  in  which  a  vehicle  is  supported  by  several  spikes  will  not 
be  explicitly  performed  were.  The  spike  will  be  considered  rigid  in  this 
derivation.  Consideration  of  spike  flexibility  will  require  more  rigorous 
definition  of  the  ground  forces. 

The  geometry  for  the  problem  is  depicted  in  Fig.  5*  The  spike  forces  and 
moments  will  be  calculated  directly  in  the  body  coordinate  system,  as  they  are 
dependent  only  on  relative  velocities. 

•The  component  of  force  along  the  spike  axis,  or  body  y  -axis,  is  immediately 
written  as  0 


"  “("I  y)  A(d)(/^r  +A^)  (3.2-83) 


where  is  the  stroking  velocity  of  the  piston  relative  to  the  body,  andA(d) 
is  the  cross-sectional  area  of  the  spike  normal  to  the  axis,  as  a  function  of 
penetration  depth  d  .  The  leading  edge  of  the  spike  will  be  pointed;  hence 


for  a  conical  point,  where  A<>  is  the  maximum  cross-sectional  area,  /£  is  the 
apex  angle,  and  the  penetration  depth  is  given  by 


(3.2-85) 


Time  is  considered  to  start  on  Impact. 

Neglecting  the  sharpened  point,  the  cross-sectional  area  normal  to  the 
body  X  -axis  isdA-Dd^o,  where  D  is  the  spike  diameter  and dy,o  is  an  incre¬ 
mental  distance  measured  along  the  spike, being  measured  from  the  tip.  The 
lateral  force  on  this  area  is  * 
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(3.2-86) 


(>o)-  2  y  D  sir x  (?tf)  d/^0 


(3.2-87) 


The  lateral  velocity /^(^is  given  by 

(^o)  =  ^x  +  +  ^  +  l  -  Jo) 

and  the  total  lateral  force  on  the  body  is 

Qxg  =  ”  i  ^  £>y'/irx  (*}o)dy0 

=-4AD[(i  +  0^+A,  +i))d-%£] 

The  distance  p.  is  from  the  origin  of  the  body  coordinate  system,  or  center  of 
mass  of  the  undeflected  body,  to  the  center  of  mass  of  the  piston. 

The  moment  about  the  body  ^  -axis  due  to  dQxC  is 

dN<JG  “?°)dQXG 

which  integrates  to 


2~}  (3*2-8 £ 


(3.2-89) 


=(^A>+i-cj4)Q„ 


(3.2-90) 


The  moment  about  the  piston  center  of  mass  is 


t\ro  -Ut-d/l) Q 


XG 


BD  O.d  3 
V  8 


(3. 2-91) 


3. 2. 5.4  Skis  and  skids 

The  elements  referred  to  as  skis  or  skids  will  not  be  distinguished  here, 
as  they  have  essentially  the  same  properties.  Either  will  be  a  device  of  some 
length  attached  to  the  lower  end  of  a  gear,  oriented  such  that  the  rearward  end 
contacts  the  ground  first.  The  element  will  generally  be  allowed  to  rotate 
relative  to  the  gear,  and  will  have  some  sort  of  rotational  spring  damper  to 
control  this  mcMon. 

In  the  case  of  rigid  ground,  the  ground  forces  will  be  derived  from 
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elimination  of  a  degree  of  freedom  for  motion  of  the  contacting  element,  as  in 
the  case  of  the  surface  pads  on  a  rigid  ground.  If  the  ground  is  not  rigid, 
ground  forces  may  be  developed  using  the  ground  penetration  development  ex¬ 
pressed  in  Paragraph  3.2. 5.3. 

Skis  and  Skids,  Rigid  Ground 

If  the  ground  is  to  be  considered  rigid,  the  ground  forces  cannot  be 
derived  from  properties  of  the  ground,  and  must  be  determined  from  the  motion 
of  the  contacting  element,  as  was  done  in  Section  3.2.5.I  for  surface  pads  on 
a  rigid  ground. 


The  rotational  motion  of  the  element  about  the  axle  as  a  rigid  body  is 
governed  by  Equation  2.8-36.  However,  if  the  moment  of  inertia  of  the  elemert 
about  the  axle  is  sufficiently  small,  the  inertial  terms  may  all  be  ignored, 
and  the  restraining  moment  about  the  axle  may  be  set  equal  to  the  moment  about 
the  axle  from  the  ground  forces.  The  degree  of  freedom  of  the  rotation  is  then 
eliminated  as  a  variable  in  order  to  yield  the  ground  force.  If  the  resv.rair.ing 
moment  is  only  for  the  purpose  of  opposing  rebound  and  not  to  absorb  laniing 
impact,  the  forces  involved  before  the  clement  flattens  on  the  ground  are 
actually  small  at  any  rate,  The  element  must,  of  course,  be  considered  rigid 
in  this  case. 


The  geometry  is  depicted  in  Fig.  6.  He  force  Wm  is  the  component  of  the 
total  ground  force  acting  to  rotate  the  element  about  its  axle.  This  component 
is  not  generally  in  a  plane  normal  to  the  ground  plane,  as  the  figure  might 
indicate.  It  is  first  necessary  to  develop  thr  form  for  this  component.  The 
axle  is  assumed  to  be  parallel  to  the  body  -axis  initially.  The  element  may 
rotate  about  the  gear  center  line,  or  the  <y-axis,  through  an  angle /3 
(assumed  small)  due  to  torsional  elasticity.  Its  angle  relative  to  the  x  -axis 
of  the  gear  coordinate  system  is  *  ,  consistent  with  the  definition  used  in 
forming  Eq.  2.8-39  governing  component  rigid  body  motion  about  a  line.  The 
coordinate  system  in  which,  instantaneously,  the  X  -axis  lies  parallel  to  the 
element  and  the  -y,  -axis  lies  parallel  to  the  axle,  is  derived  by  the  product 
of  transformations  due  to  the  rotations/^  ,  i £  from  the  gear  coordinate  system; 


COS* 
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’  1 
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(3.2-92) 


The  unit  vectors  in  the  gear  coordinate  system  may  be  written  in  terms  of  those 
in  the  ground  coordinate  system  by 


j'i'l  r  fi 

<J|'  >  =  |y]  tr]<J 

bo  b<, 


The  components  of  the  ground  force  in  the  contacting  element  coordinate  system 
are  then  related  to  those  in  the  ground  coordinate  system  by 
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1 


(5*2-94) 


[v]  [rjfsA 

LVJ 

The  component  Wq  xies  along  the  element,  the  others  at  right  angles  to  the 
element.  The  component W(,  may  be  used  to  find  the  applied  moment  for  the 
torsional  motion  in  the  gear.  Only  the  component Wn  contributes  to  the  elanent 
rotation  about  the  axle.  It  is  assumed  that,  the  moment  due  to  this  force, J?  W#, 
is  balanced  by  the  restraining  moment 


COStf  /?COS  It  —  SIN  71 

v*  \  0 

(wL  )= 

|wj 

sin  n  ft  cos  n 

(5.2-95) 


It  is  assumed  that  the  drag  force  on  the  trailing  end  is  opposing  the  horizontal 
velocity  of  the  end,  and  that  its  magnitude  is  the  product  of  the  vertical  ground 
reaction  with  the  coefficient;  of  friction  between  the  element  and  the  ground. 

The  derivation  producing  Eq.  3.2-44  is  used  here,  with  the  pad  velocity  components 
replaced  by  the  components  of  velocity  of  the  element  trailing  edge,  so  that 


The  Eqs.  3.2-94, 


96  may  be  combined  to  yield  Wjg  in  terms  of  V&  : 


W„=VG  i 


(5.2-97) 


where 
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This  result  may  be  substituted  into  Eq.  3*2-95  and  the  vertical  ground  reaction 
found  in  terms  of  the  restraining  moment.  Use  of  Ea-.  3*2-98  then  yields  the 
remaining  components  of  the  ground  reaction.  Thus,  for  E  A  <0, 
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This  result  is  quite  similar  to  that  obtained  in  the  paragraph  on  surface  pads 
for  the  case  of  rigid  ground.  The  form  is  valid  only  until  the  element  flattens 
out  on  the  ground.  Thereafter,  the  forces  are  derived  exactly  in  the  manner  of 
the  section  on  surface  pads,  except  that  the  coefficient  of  friction  will  generally 
he  different,  and  the  components  of  the  pad  horizontal  velocity  become  those  of 
the  element  axle  velocity.  The  axle  height  2/,  may  be  obtained  by  integrating  the 
axle  velocity  normal  to  the  ground,  given  by  Eq.  2.10-6.  The  ground  forces  on 
the  body,  expressed  in  the  body  coordinate  system,  are  then  given  by  Eq.  3 .2-35  • 

The  element  must  be  considered  rigid  in  the  previous  derivation.  If 
flexibility  is  to  be  included  in  the  same  direction  as  the  component  rigid  body 
motion  about  a  line,  the  problem  may  be  approached  in  the  same  manner  as  that 
in  Paragraph  4. 2. 5. 2, which  discusses  articulated  gears.  The  ground  force  is  then 
defined  in  the  following  paragraphs. 

Skis  and  Skids.  Soft,  (ire-mud 

Prior  to  the  time  when  the  element  flattens  onto  the  ground,  tn.  ground 
forces  may  be  derived  from  the  general  form  for  ground  penetration  forces 
expressed  by  Eq.  3.2-82,  which  is  written  in  differential  form  as 


/ird  A 


(3.2-100) 


vhere  the  velocity  or  is  that  relative  to  the  ground  and  normal  to  the  differential 
area  dA.  The  geemetry  is  depicted  in  Fig.  7. 


The  parameter^  is  measured  along  the  element  from  the  trailing  end;  its  value 
at  the  surface  of  the  ground  is  D  .  This  "planing"  length  is  related  to  the 
penetration  depth,  d  ,  the  axle  height,  and  the  distance  from  the  axle  to  the 
trailing  end  by 


D 


(3.2-10J.) 


The  element  is  assumed  to  have  a  "planing"  width  W,  ,  whioh  is  independent  of 

,  and  a  constant  lateral  vidth  ar~.  The  forces  normal  to  these  areas  are 
those  in  the  contacting  element  coordinate  system,  defined  in  the  paragraph  on 
rigid  ground.  The  differential  force  normal  to  an  element  of  r-rea  ar,ctjT  is 
then 


*-(f  |^v.(y)«n  d? 

and  that  normal  to  an  area  oJ^dJis 

dWL  =-(f  (?)“•* 


(3-2-102) 


(3.2-103) 
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She  differential  of  friction  force  aue  to  these  forces  is  parallel  to  the  element; 


d  WD  -s*  (dWM  +  dWL) 


(3.2-104) 


These  three  terms  define  the  force  on  an  elemental  lengthd  f.  The  velocities 
at  that  point  are 


QQ*/3rNA+*(i|-2') 

(y)  =  /^LA‘V^COS^(^''cr) 


The  transformation  from  the  gear  coordinate  system  to  that  of  the  contacting 
element,  expressed  by  Eq.  3*2-92,  yields  the  components  of  axle  velocity /^hA 
/Z'la-  Tt0  axle  velocity  in  the  eround  coordinate  system  is  derived  from  Eq.  J 
2.10-6,  and  may  be  transformed  to  the  gear  coordinate  system.  The  desired  axle 
velocities  are  given  by 


in  terms  of  the  panel  point  displacement  and  velocity  of  che  axle  expressed  in 
the  gear  coordinate  system.  These  velocities  ire  not  dependent  or,  the  pa"  jueter 
JT  ;  hence  the  forces  are 
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(3.2-109) 


w.  -  -(?¥)“->  (f  »  +  y(£<-  -r)]D 

WL=  +/^cos ?(i,-  -^]d  (3-2-110) 


w„  =^(wN+w,.) 


(3.2-111) 


in  terms  of  the  planing  length  D  .  From  Eq.  3.2-101,  the  planing  length  is 
determined  by  the  penetration  depth.  This  may  be  obtained  by  integration  of  the 
velocity  of  the  trailing  end  of  the  element  normal  to  the  ground.  The  addition 
of  the  velocity  of  the  trailing  end  relative  to  the  axle  to  the  form  in  Eq.  2. 
10-6  yields  the  desired  quantity.  Thus, 


(3.2-112) 
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COS# 

ft  COSH 

-  SIN  # 


-/3 

J 

0 


•  SIN  ft 
/^SIN  ft 
COS# 


(3.2-114) 


These  forms  are  valid  only  until  the  element  becomes  horizontal,  or  parallel 
to  the  ground.  It  has  been  assumed  that  the  element  rotates  about  an  axle,  th? 
rotation  being  described  by  the  angle  f[  .  If  this  angle  is  allowed  to  vary, 
the  forms  above  for  ground  forces  will  generally  be  valid  until  the  axle  height 
becomes  zero.  It  is  suggested  that  when  this  occurs,  the  forces  may  then  be 
obtained  by  a  form  which  assumes  that  the  axle  vertical  velocity  is  zero.  This 
is  the  form  used  in  deriving  the  forces  on  a  surface  pad  in  the  case  of  a  rigid 
ground,  which  should  hold  in  the  case  of  a  ski  or  skid  of  large  surface  area. 

If  the  element  is  actually  rigidly  affixed  to  the  gear,  or  the  rotational 
spring-damper  is  sufficiently  strong,  the  element  will  not  become  parallel  to 
the  ground  and  the  above  forms  will  remain  valid. 

3.S.5.5  Gas-Filled  bags 

Although  this  formulation  of  the  landing  impact  problem  does  not  describe 
gas-filledbag  inertial  properties,  the  bag  and  vehicle  mass  may  be  combined  as 
a  single  rigid  body  to  describe  the  inertia.  At  present,  the  only  analytical 
work  of  any  rigorous  type  lias  been  restricted’  to  vertical  alightment.  Ref.  15 
is  an  example  of  this  type  of  effort.  In  the  notation  of  this  report,  the 
equation  for  vertical  motion  used  there  is 


M  (3.2-115) 

The  external  forces  considered  in  the  reference  are  gravitational,  atmospherical 
drag,  and  ground  force  from  the  gas-filled  bag.  The  latter  is  given  by 


Q>0  ="A(P-P0)  (3.2-116) 


where  p  is  the  bag  total  pressure  and  Po  is  the  atmospheric  pressure.  The 
pressure  variation  in  the  bag  is  assumed  adiabatic,  for  the  case  of  no 
bleeding,  so  that 


(3.2-117) 


82 


Then  for  a  cylindrical,  non-bulging  bag,  the  volume  is  related  linearly  to  the 
body  deflection  after  touchdovn; 


(5.2-118) 


vhere  h  is  the  undeflected  cylinder  height  and  d  is  the  deflection  of  the  bag, 
obtained  by  integration  from  touchdovn  time. 

Additional  complications  in  the  form  used  for  the  pressure  in  the  bag  vil’ 
result  if  the  bag  shape  varies  or  if  a  gas  orifice  is  used,  likavise,  the 
contacting  area  may  be  a  variable.  These  points  are  discussed  in  the  reference. 

It  is  pointed  out  that  for  stability  requirements,  a  multiple  beg  system  is  usually 
helpful. 

Generally,  then,  it  is  necessary  to  define  the  variation  in  the  bag  shape, 
contact  area,  and  bag  pressure  in  terms  of  the  geometry  of  the  bag,  the  bag 
deflection,  and  some  polytropic  compression  form.  The  simplest  approach  is  to 
assume  that  the  bag  shape  remains  constant.  The  contact  a'rea  is  then  defined  by 
the  geometry  and  the  bag  vertical  deflection.  If  the  bag  shape  varies  with  pres¬ 
sure,  the  elastic  properties  of  the  bag  mur-t  be  included. 

If  the  gas  bleeding  orifice  is  controlled  srch  that  the  bag  pressure  remains 
constant,  then  the  form  for  the  vertical  force  is  fairly  simple.  If  the  rate  of 
bleeding  does  not  produce  this  result,  then  the  pressure  variation  must  be 
described  in  some  manner  in  terms  of  the  bag  deflection  and  deflection  rate.  The 
gas  flow  characteristics  necessary  to  define  the  pressure  variations  will  not  be 
written  here. 

Suppose  a  set  of  independently  operating  bags  are  located  on  a  plane  surface 
on  the  vehicle.  For  small  angles  from  the  vertical,  the  restraining  force  from 
each  bag  may  be  considered  to  be  nearly  parallel  to  the  body  ^  -axis;  each  is 
given  by 


Vs  =  -A(P-Po) 


(3.2-119) 


after  toucudown,  and  is  zero  beforehand.  The  vertical  height  of  the  bags  defines 
the  touchdown  times,  as  in  the  section  on  surface  pods.  The  drag  force  on  °".ch 
bag  is  assumed  to  have  conponents  given  by 


Dg,  =  S#  Vq 


Xb 

IVJ 


(5.2-12G) 
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(3.2-121) 


where /t  Is  the  coefficient  of  friction  between  the  bag  and  the  ground.  Body 
forces  are  then  given  by  Equation  3-2—3^  and  body  moments  by  Equation  3.2-37. 
The  instantaneous  positions  of  the  contacting  points  should  be  interpreted  as 
the  positions  of  the  center  points  of  the  bag  ends.  The  deflection  Py 
of  each  bag  would  be  the  stroking  deflection,  and  the  lateral  deflections  Pac  , 
Py  may  be  derived  from  the  drag  forces  and  bag  lateral  spring  rates.  These 
spring  rates  will  probably  be  dependent  on  vertical  deflection.  The  vertical 
deflection  of  each  bag  will  be  found  by  integrating  the  velocity  normal  to  *r.z 
ground  of  the  center  point  of  the  bag  attachment  area.  The  above  statements 
are  sufficient  to  work  a  stability  problem  for  landing  on  a  set  of  gas-filled 
bags,  provided  that  the  variation  with  all  parameters  of  the  bag  area  and 
pressure  in  Equation  3*2-116  arc  defined. 

3*3  APPLIED  FORCES  Oh  COMPONENT  RIGID  BODY  MOTION 


3*3.1  General 


The  purpose  of  this  section  is  to  define  the  applied  forces  which  enter 
Equations  2.8-18,  40  for  component  rigid  body  motion  for  the  particular  cases 
or  gear  piston  stroking  and  bogie  rotational  motion.  The  forces  are  of  two 
types,  ground  forces  and  stroking  forces. 

For  gear  stroking,  the  total  force  along  the  piston  axis  is  given  by 


(3.3-1) 


where  Q,'fi  is  the  component  of  the  ground  force  along  the  stroking  axis, 
is'the  total  stroking  force. 


and 


For  bogie  rotation,  the  total  moment  is  due  to  the  ground  forces  and  the 
rotational  damper  forces.  The  ground  forces  are  derived  later.  The  damper 
moment  may  be  expressed  by 


N  =  N  (#,*) 


'5.3-2) 


where  it  is  the  bogie  rotational  angle..  The  relation  between  this  angle  and 
its  time  derivative  to  the  actual  damping  mechanism  is  considered  briefly  here. 
If  the  damping  mechanism  is  a  stroking  element  affixed  to  the  bogie  element  and 
the  gear  piston,  the  exact  relation  will  be  non-linear  .  The  gear  coordinate  axes 


e>: 


have  been  translated 


y 


to  the  bogie  axle  for  convenience.  The  bogie  is  originally  at  an  angle  from  the 
DC  -axis,  defining  the  position  of  zero  stroke  for  the  rotational  spring  damper. 

If  the  bogie  rotates  only  through  a  small  angle,  the  restraining  moment  is  re¬ 
lated  to  the  spring  damper  force  by 

N~X-s\h/3  Fso 


The  force  may  be  derived  from  the  paragraphs  or  springs  and  damping  mechanisms  by 
replacing  the  element  stroke  by 


for  small  rotations.  Hence 


X  ~{Ir  siN/^ )  -q 


For  large  rotations,  these  expressions  must  be  replaced  by  forms  nonlinear  in  the 
angle  ^  . 

3.3-2  Ground  Forces 


The  Equation  3.2-35  defines  the  force  on  the  body  due  to  interaction  with 
the  ground  of  a  single  contacting  element.  This  force  may  be  transformed  into 
the  gear  coordinate  system  to  yield  the  conponents  needed  for  gear  stroking  and 
bogie  rotation.  The  force  on  a  single  contacting  element  is  written  in  uhe  body 
coordinate  system  as 
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This  is  written  in  the  gear  coordinate  system  as 


M  O’] 


(5.3-4) 


If  there  are  several  contacting  elements  on  a  single  gear,  this  form  is  summed 
over  those  elements; 


1  /  V  fD*l 

=M  Cr]>  s’ 

J  -  K  J 

•'g 


This  form  yields  the  ground  forces  in  the  gear  coordinate  system.  The  forces 

v'z/y,  are  used  in  the  calculation  of  bogie  moments,  and  the  force 
Q?S  is  used?  to  define  the  ground  force  for  piston  stroking. 

3.3.3  Stroking  Forces 

The  purpose  of  this  paragraph  is  to  define  the  stroking  force  This 

force  is  the  summation  of  all  the  forces  on  the  piston,  along  the  strolling  axis, 
except  for  ground  forces.  It  is  important  to  realize  that  the  formulation  is 
sufficiently  general  that  the  effects  of  any  ’mown  type  of  energy-absorbing 
me c ha:.. Ism  may  be  included  simply  by  defining  the  stroking  force  properly. 

The  remainder  of  Section  3.3  will  be  devoted  xo  the  stroking  forces. 
Obviously,  not  all  the  various  types  of  shock-absorbing  devices  can  be  considered. 
Those  to  be  included  are  a3  follows.  Paragraph  3*3»3.1  covers  hydraulic  forces. 
Variations  with  oil  compressibility,  metering  pins,  and  relief  valves  are  included. 
In  Paragraph  3>3*3.2,  springs  of  the  mechanical,  pneumatic,  and  liquid  types  are 
considered.  Paragraph  3. 3. 2. 3  discusses  bearing  friction  forces  for  a  particular 
configuration.  The  piston  bottoming  or  restraining  forces  are  defined  in  Paragraph 
3. 3. 3. 4.  Crushable  materials  are  discussed  in  Paragraph  3*3.3»5>  <ukL  gas  expulsion 


devices  in  Paragraph  3«3-3.6  complete  the  types  of  stroking  forces  to  be  develop¬ 
ed  in  the  report. 

In  order  to  write  the  relatione  defining  these  forces  in  compact  form,  the 
displacement  of  the  piston  relative  to  the  cylinder  vill  be  used.  This  dis¬ 
placement,  referred  to  as  stroke  and  given  the  symbol  -d  ,  is  measured  from 
the  fully  extended  position  and  is  positive  as  the  piston  enters  the  cylinder. 

The  equation  governing  the  piston  motion  is  written  in  terms  of  the  piston 
displacement  relative  to  the  body  coordinate  system,  A  V  *  which  is  positive 
along  thepf  -axis  of  the  gear  coordinate  system.  Let  t'ne  point  at  which  the 
cylinder  is  affixed  to  the  vehicle  be  labeled  B  ,  and  its  total  displacement 
IPft  .  These  three  variables  are  then  related  by 

-  -Ay  •  IK/  =  -Ay  +R^a  (3.3.6) 

This  result  indicates  that  stroking  may  occur  due  to  elastic  deformations  in  the 
vehicle.  This  result  must  be  used  in  all  the  forcing  functions  in  this  section, 
if  the  vehicle  is  considered  to  be  flexible.  The  stroking  velocity  may  be  obtained 
by  taking  the  time  derivative  of  Equation  3*3-6. 

3«3*3.1  Hydraulic  farce.  pH 


When  the  piston  shown  in  the  figure  is  displaced 
with  a  velocity  A  ,  a  pressure  field  is  created  in 
the  lower  chamber  wh^ch  resists  this  motion.  In  a 
rigorous  definition  of  this  pressure  field,  where 
consideration  is  given  to  viscosity  and  compressibility 
of  the  fluid  and  the  unsteady  aspects  of  the  flow, 
one  is  faced  with  a  formidable  and  in  most  cases 
unsolvable  formulation.  Satisfactory  definitions 
of  loads  in  the  piston  have  been  achieved,  however, 
by  semi-empirical  means  where  the  form  for  the 
pressure  in  the  piston  is  defined  by  Bernoulli's 
principle  for  an  ideal  fluid.  Neglecting  the 
static  bead,  the  velocity  of  the  Jet  stream  at 
a  point  outside  the  orifice  at  which  the  stream¬ 
lines  are  parallel  is  given  by: 


2  Pa 

/°H 


(3.3-8)  * 


where 


• 

xi  :  stroking  velocity  of  the  piston 

p  3  pressure  upstream  of  the  orifice  where  the  streamlines  are  paral.1el 
fS  s  pressure  in  the  air  chamber 
density  of  the  fluid 

•Equations  on  this  page  misnurmered. 
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In  actual  practice  the  velocity  /ir  is  never  attained  due  to  the  dissipation  of 
energy  in  overcoming  the  resistance  to  flow.  The  ratio  of  the  actual  velocity 
to  the  velocity  given  by  the  Bernoulli's  relation  is  defined  by  the  coefficient 
of  velocity,  Ov  *  Equating  the  rate  at  which  the  volume  of  fluid  is  displaced 
by  the  piston  to  the  rate  at  which  it  is  discharged  through  the  orifice. 


AH-^  -  Ay  0 \/ sir 


(3.5-9) 


where  is  the  hydraulic  area  of  the  piston,  and  A^.  Is  the  area  of  the  jet 
stream  ac  a  point  where  the  streamlines  are  parallel.  The  ratio  of  to  the 
orifice  area  A#is  given  by  the  coefficient  of  contraction,  Cc  •  For  most 
hydraulic  dampers,  is  small  compared  to  /V  ,  and  combining  Eqs.  3.3-7, 

8  yields 


2  •  2 


2(CcCv)2Aw2 

/*H 


(P-Po) 


(3.3-10) 


The  hydraulic  force  FH  is  then  given  by: 


Fh 


=  (P~Pa)A„ 


/oH  AH3i.  2 

=i(cX)* 


(5.3-11) 


where  the  orifice  coefficient  Cd'Is  the  product  of  Cc  and  C.v  . 

The  value  of  the  coefficient  Cd  is  dependent  upon  the  size  of  the  orifice, 
the  shape  and  finish  of  the  orifice  face,  the  kinematic  viscosity  of  the  fluid, 
the  velocity,  and  the  motion  in  the  fluid  approaching  the  orifice  which  causes  a 
dependence  on  the  length  of  the  oil  column  remaining  in  the  strut.  The  relation¬ 
ship  between  size,  velocity  and  kinematic  viscosity  can  be  expressed  as  Reynold's 
number.  Thus,  the  /alue  of  the  orifice  coefficient  for  a  given  orifice  can  then 
be  said  to  be  dependent  upon  Reynold's  number  and  piston  stroke. 

The  results  of  an  experimental  study  of  orifice  coefficients  in  a  small 
oloo-piidumatic  strut  with  a  constant  orifice  are  contained  in  Ref.  16.  This  test 
investigated  the  effect  on  orifice  coefficient  of  variations  of  Reynold's  number 
in  the  range  from  9*500  to  66,500.  The  results  of  the  test  were  sunm»ri,.od  by 
an  empirical  relation  between  orifice  coefficient,  stroke  and  stroke  velocity. 

The  variation  between  the  minimum  and  maximum  values  of  the  coefficient  for  the 
tests  ranged  from  0.86  to  0.93*  The  final  conclusion  from  the  experiment  was  that 
an  average  value  of  orifice  coefficient  could  be  used  as  a  constant  to  lei  ermine 
strut  loads. 
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Effect  of  Oil  Compressibility 


When  the  contact  velocity  of  the  strut  is  high  (generally  above  15  feet  per 
second), the  compressibility  of  the  hydraulic  oil  and  in  certain  cases  the  volumet¬ 
ric  expansion  of  the  strut  will  have  an  effect  on  the  loan  time  history'  • 

The  following  treatment  of  this  effect  has  improved  substantially  the  correlation 
between  analysis  and  drop  test  results  at  Chance  Vought. 

The  change  in  pressure  d  P  in  the  fluid  due  to  compressing  the  volume  by 
an  amount  dv  is  given  by 


dP=K^ 


(3.3-12) 


where  |6  would  be  an  equivalent  bulk  modulus  expressing  both  the  compressibility 
of  the  fluid  and  the  elasticity  of  the  strut.  Since  the  volume  change  is  small 
compared  to  the  total  volume,  the  differentials  can  be  replaced  by  finite 
differences 


AP 


(5.3-13) 


The  term  Av  is  the  difference  between  the  volume  swept  by  the  hydraulic  area 
and  the  volume  expelled  through  the  orifice: 


dl'AV) 

dt 


=  Ah  A-Cd  An  /v 


(3-3-14) 


For  a  constant  Ah, 


AV  =  Ah  'ds  ~  f~ Cd  An  w 


(3.5-15) 


The  volume  V  is  the  volume  of  the  oil  chamber  at  any  time,  hence 


V  =  V0-A„^ 


.3-16) 


where  Vo  is  the  volume  of  the  oil  chamber  when  the  strut  is  fully  e:  tended. 
Thus,  AP  becomes 
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Ap-U  Ah^»  Jq  Ah  /V  dt 
Vc  -  Ah  ^ 


15-5-17) 


The  term  A  P  represents  the  total  change  in  pressure  from  time  t  =  0  to  time 
"t,"  so  that 


AP  »  P-Pac 


(5.3-18) 


where  P  is  the.pressure  in  the  oil  chamber  and^P^  is  the  initial  air  pressure. 
The  pressure  AP  is  due  to  hydraulic  pressure  and  air  pressure  ; 


AP  =  P„+-PA”PA0 


(3.3-19) 


thus  the  expression  for  becomes 


E5  -  J/  CpAn  /IT  dt  = 

P«~K  Vo-A„^  ~p 


Pa-^Pao 


(3.3-20) 


The  velocity /V  is  defined  as  the  Bernouilian  velocity. 


/2Ph 


(3.3-21) 


The  hydraulic  force  F h  is  given  by 


Ph  -  Ah  Ph 


Orifice  and  Relief  Valve  Combination 


(3.3-22) 


The  relief  valve  operates  on  the  principle  that  when  the  hydraulic  pressure 
reaches  a  predetermined  value,  a  valve  cracks,  introducing  additional  orifice 
area.  Theoretically,  this  method  of  controlling  the  load  in  the  strut  >s  superior 
to  the  metering  pin  approach  since  it  functions  from  load  level  rather  than  stroke 
There  are,  however,  sufficient  design,  qualification  and  manufacturing  difficulxie 
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to  detract  from  its  theoretical  superiority. 

Writing  a  general  expression  for  a  hydraulic  strut  employing  a  relief 
valve  is  difficult  since  designs  vary  widely.  The  following  is  offered  as  an 
example  to  illustrate  the  concept.  The  geometry  is  depicted  in  Fig.  8. 

The  piston  of  area  strokes  with  a  velocity  Jl,  ,  metering  fluid 
through  an  orifice  with  area  A 4  •  When  the  pressure  acting  on  the  area  A, 
overcomes  the  downward  force  pRo  on  the  valve  due  to  the  spring,  the  valve 
moves  upward,  uncovering  an  area  A3(xr)  .  The  pressure  acting  on  the  relief  valve 
before  cracking  is  PH  .  The  force  required  for  cracking  is  therefore  PH  A.  . 


Figure  8.  Hydraulic  Damper  with  Relief  Valve. 

The  force  acting  on  the  valve  due  to  the  spring  is  given  by 

Fr  -Fr0  +  K.rXr  (3.5-23) 

where  Kg.  is  the  spring  rate  of  the  relief  valve  spring. 

If  the  compressibility  of  the  fluid  and  the  mass  of  the  relief  valve  can 
be  neglected, 
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(3-3-24) 


n  -  Sl  Ah  i  2. 

Ph  ~ 2C02 A42  ^ 


pH  A.  <  Fr, 


A3  (xr)  limit  value  =  A 5 


where 


-PhA,-FRc 

kr 


(3.3-26) 


The  hydraulic  force  Fh  is  then 

Fh  =  Ah  Ph  (3.3-27) 

Oil  compressibility  effects  can  be  included  in  the  same  manner  as  discussed 
previously. 

Extension  Stroke 


In  general  the  damping  characteristics  of  a  self-positioning  strut  made  up 
of  a  spring  and  hydraulic  damper  ii.  series  will  be  different  during  the  extension 
stroke  tlian  during  compression.  The  extension  characteristics  are  governed  by 
two  requirements:  oil  must  be  returned  to  lower  chamber  fairly  rapidly  for 
another  energy  absorption  cycle;  and  adequate  damping  must  be  provided  to  reduce 
bottoming  loads  at  the  zero  stroke  position  during  rapid  extension. 

The  energy  available  for  the  extension  stroke  is  that  stored  in  the  gear 
spring  which  may  be  pneumatic,  liquid  or  mechanical. 


The  form  for  the  hydraulic  force  is  the  same  as  for  the  compression  stroke 
except  the  constants  are  different 


Fh  =•  ~ 


/*h  Ah53 

2(Co  An)J 


(3.3-28) 
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where  the  parameters  are 

Ahs  »  hydraulic  area  associated  with  the  extension  stroke 
(£0  Ah)s  =  effective  orifice  area  associated  with  the  extension  stroke 

The  relations  for  the  hydraulic  force  previously  derived  are  summarized 
below.  The  relation  defining  the  stroke,  ,  from  Eq,.  3.3-6  shov'vi  be  kept  in 
mind. 

Incompressible  oil,  metering  pin 


_/VAhj^2 

2(Cq  An)  2 


■A*  ^  0 


(3.5-29) 


c  -  ^Ah,U2 
2(Cd  An)s2 


0 


Incompressible  oil,  relief  valve 

.  *  3  1  2 

r  _  Ah  -A 

Fh“  z(c„a4)z 


X  >  0 
Ph  Aj  >Fr0 


r-  /^h  Ah  A, 

2Cd2[A4+-A5(x4)]2 


A  >  0 

A|  >  FRo 


(3.3-30) 


(3.3-31) 

(3.3-32) 


Ph  A 2  ~  Fro 


F 


H  = 


a  3:2 
Ahs  _ 

2(Cd  An)s2 


(3.3-33) 


^  <  0 


(3.3-34) 


Compressible  oil,  metering  pin 


Fh  -  Ph  Ah 

A»  Co  An  /IT dt 


f 


Ph  =  1C 


V0-Ah^ 


-Pa+Pao  >  ^  >  0 


(3.3-35) 


Fh  = 


-v^H  ^HS  ^ 

2(C0  An)s2 


<  0 


(3.3-35) 


The  parameter  An  is  in  each  case  the  metering  function,  which  may  be  a  function 
of  stroke,  and  may  be  of  a  form  for  the  return  stroke  different  from  that  of  the 
compression  stroke. 


3.3*3»2  Spring  forces 

A  spring  is  by  definition  an  elastic  body  or  device  that  returns  to  its 
original  shape  after  being  distorted.  A  spring  when  distorted  will  generate 
a  restoring  force  that  is  functionally  related  to  the  displacement  of  the  spring 
along  a  specified  axis.  If  energy  is  dissipated  during  the  distortion-recovery 
cycle  the  functional  relationship  between  restoring  force  and  deflection  will 
be  double  valued,  as  shown  below,  where  the  shaded  portion  represents  the  energy 
dissipated: 


The  general  expression  for  a  spring  exhibiting  the  above  load-deflection 
character la  lie  is  given  by 


C*4)  X  }  0 

Fs  s2  ('&)  X  <  0 


(3.3-57) 


Since  the  spring  itself  has  inertia,  the  load-deflection  characteristics 
are  dependent  somewhat  on  the  rate  of  loading  of  the  spring.  These  effects, 
in  general,  will  be  small  for  the  anticipated  applications  of  this  report  and 
not  been  considered* 
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Three  types  of  springs  are  considered  in  the  succeeding  paragraphs: 
mechanical,  pneumatic,  and  liquid. 


Mechanical  Springs 


Mechanical  springs  occur  in  almost  endless  variety,  from  simple  helical 
springs,  vhich  possess  linear  load-deflection  characteristics,  to  the  mere 
complicated  Belleville  and  ring  springs  vhich  exhibit  non-linear  characteristics 
and  significant  energy  dissipation.  Mechanical  springs  lend  themselves  to 
fairly  accurate  analytical  description  since  they  vill  he  deformed,  in  general, 
only  in  the  elastic  region  of  the  material.  Detailed  formulas  of  spring  load- 
deflection  characteristics,  however,  vill  not  he  derived  in  this  report  since 
these  relations  are  veil  documented  in  design  manuals,  texts  and  manufacturers  * 
literature. 

For  the  purposes  of  this  report  it  vill  he  assumed  xhat  the  force-deflection 
characteristics  for  a  mechanical  spring  can  he  expressed  as: 


F,  -/*<  U) 
Fs  (A) 


^>0 

^<0 


(3-5-38) 


For  springs  in  vhich  hysteresis 
hy  a  single  function. 


is  negligible 


and  the  cycle  is  defined 


Assuming  the  functional  relation  for  is  known, yd  must  still  be  ex¬ 
pressed  in  terms  of  the  variables  of  the  analysis,  as  in  Eq.  3.3-6. 


Pneumatic  Spring 


A  pneumatic  spring  is  one  that  depends  upon  the  compressibility  of  a  gas 
to  generate  a  restoring  force.  The  lav  governing  the  compressibility  of  a  gas 
in  a  closed  container  is  given  by 


CONSTANT 


(3.3-39) 
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where 


P/so  /-  initial  pressure  in  cylinder 

V0  ~  initial  volume 

P/s  '-pressure  -vt  stroke 

V  —volume  at  stroke  11  A."-  y^-  ^  -4- 

VI  /-exponent  that  indicates  the  exact  polytropic  nature  of  tne  compression 
/-pneumatic  area 

Substituting  the  relation  for  V  into  the  above  equation  the  expression  for 
becomes 

Pa.v.0 

and  the  force  Fa  acting  downward  on  the  piston  is 


Pa 


Fa=Aa  (Pa  -Pj 


(3.3-41) 


The  atmospheric  pressure?^  may  usually  be  neglected  in  comparison  with  the 
cylinder  pressure. 


The  value  of  the  exponent  fb  to  be  used  depends  on  how  much  heat  is 
transferred  to  and  from  the  gas.  If  the  compression  or  expansion  takes  place 
rapidly,  such  that  little  heat  is  transferred  from  and  to  the  gas,  the  process 
can  be  assumed  adiabatic  and  rv  becomes  the  ratio  of  the  specific  heat  of  the 
gas  at  constant  pressure  to  the  specific  heat  at  constant  volume.  For  dry  air 
IT_>  =  1.4o6.  If  the  compression  or  expansion  process  is  such  that  the  temper¬ 
ature  of  the  gas  is  unchanged  (isothermal), lbs  1.  In  the  general,  polytropic 
case,  IT_>  must  be  determined  from  an  analysis  of  the  thermodynamic  process. 

When  the  pneumatic  spring  discussed  above  is  part  of  an  oleo-pneumatic 
strut,  the  thermodynamic  process  is  further  complicated  by  the  cooling  action  and 
vaporization  of  the  oil  spray.  The  net  effect  of  this  spray  is  to  cause  the 
thermodynamic  process  to  approach  isothermal.  Experiments  to  evaluate  the 
exponent  rv  for  an  oleo-pneumatic  strut,  reported  in  Ref.  17,  indicate  an  average 
value  of  rv»  =  1.06  would  adequately  represent  the  compression  process  for  the  impacts 
investigated.  For  most  practical  analyses  it  is  sufficiently  accurate  to  choose 
fU  =  1,  for  which  case 


h 


PaoVo 
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(5-3-42) 


Here  again  the  variable  "A>"  must  be  written  in  terms  of  the  variable  of  the 
analysis,  as  given  by  Eq.  3.3-6. 

Liquid  Spring 

In  recent  years  a  number  of  aircraft  landing  gears  have  been  designed 
solely  on  hydraulic  principles,  where  the  function  of  the  pneumaxic  spring  is 
replaced  by  compressing  the  hydraulic  fluid.  Usually  some  mechanical  advantage 
is  involved  between  wheel  stroke  and  shock  absorber  stroke  since  the  liquid 
spring-shock  absorber  is  inherently  a  short  stroke  device.  Two  spring  and 
damper  concepts  were  examined.  Although  these  differ  materially  from  a  design 
and  utility  standpoint,  they  both  submit  to  the  seine  analytical  treatment. 


A  simple  geometry  for  a  liquid  spring  is  pictured  below.  Initially,  the 
volume  V|0  to  the  left  of  the  piston  and  the  volume  V^o  to  the  right  of  the 
piston  are  filled  with  a  compressible  fluid  to  some  initial  pressure.  (The 
mechanical  stops  are  not  pictured . ) 


The  force  on  the  piston  is  made  up  of  the  hydraulic  force  resulting  from 
the  difference  in  hydraulic  pressure  and  area  on  either  side  of  the  orifice,  the 
frictional  force  resulting  from  the  normal  pressure  of  the  seals,  and  a  force 
Q  p  which  represents  any  other  forces  external  to  the  liquid  spring  system. 

The  piston  force  Fp  is  therefore: 

FP  »  Fw  +  Ff  +Qa 

**  f?  AP  -  Pz  (Ap  “Ar) 

-Ap(F?-Pj^ARPa 

where  Ap  is  the  area  of  the  piston  and  AR  is  the  area  of  the  rod. 

The  frictional  force  Fc  will  generally  be  a  significant  contribution  to 
the  net  force  in  xhe  piston  due  to  the  severe  sealing  requirements. 

The  analytical  expressions  for  the  normal  forces  on  the  pisten  due  to  the 
sealing  pressure  will  be  functions  of  the  particular  design.  In  general  they 
can  be  expressed  as: 
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FNU  * ^(Pi  jF2.) 

Fjy/c  L  (  P,  ,  P-  ) 


\ 


The  frictional  force  on  the  piston  is  then: 

Ff  =(^u  Fnu  Fnl)  jfj 

FOR  Fw  +  /#L  Fnl<Q^F„ 

and 

Ff=-(Q^r„) 

FOE  Fnu  +  F„_  >  Q a  +  Fh  (3-3-41.) 


whereto  and  are  the  iriction  coefficients  of  the  upper  and  lower  seal 
respectively. 

The  hulk  moduli  of  oils  used  in  liquid  springs  in  general  vary  linearly 
with  pressure  throughout  the  range  of  interest.  Thus,  in  the  two  regions. 


B. 


dP, 

dV. 


(5-3-45) 


d  Fz 

dV2 


(3.3-46) 


Asoum<-  for  the  moment  that  the  contribution  to  the  change  in  total  volume  from 
the  cylinder  and  seal  elasticity  can  be  neglected.  The  instantaneous  volumes 
of  oil  on  either  side  of  the  piston  are  dependent  on  both  stroke  and  The  amount 
of  fluid  which  has  been  metered  through  the  orifice.  These  will  be  given  by 
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(3.3-47) 


V|- V/0  A p^  +  vM 
VV2o+<Var)a-Vm 


and  their  differentials  by 


dy=-Apds  +  dVM  (3.3-49) 

dVa  -  (Ap  -  AR)d^.  -  dVM  (3.3.50) 

The  volume  of  oil  metered  through  the  orifice  is  related  to  the  orifice 
coefficient  Cq  ,  the  orifice  area  ,  and  the  velocity  of  the  oil  in 
the  orifice  /V  ,  by  ?V* 


dVM  =  C0  An  SV  dt 


(3.3-51) 


(3.3-52) 


is  used.  With  these  expressions,  Eq.  3*3-35  becomes 


which  ..ay  be  written  as 


f>  _  (a  +JrFjX~ Ap^  +  C0  An  /I/") 

4  VJ0  -  Ap  ^  +  /'  CD  A„/ydt 


(3.3-53) 


(3.3-54) 


A  similar  form  may  be  developed  from  Etp.  3*3-46.  The  resulting  expressions 
which  define  the  force  acting  on  the  piston  rod  are  recapitulated  below. 
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The  inclusion  of  the  volumetric  expansion  of  each  region  due  to  cylinder  ex¬ 
pansion  and  seal  compression  is  quite  complicated.  The  change  in  volume  from 
each  of  these  effects  must  be  written  in  terms  of  the  pressure  and  the  piston 
stroke.  In  terms  of  the  geometrical  volume  change  d  Vic  due  to  cylinder  ex¬ 
pansion  in  region  one  and  dV|$  due  to  compression  of  any  seal  in  region  one, 
and  similar  terms  in  region  two,  the  above  relations  may  be  restated  as  follows. 


Liquid 


Elastic  Cylinder  and  Seals 


(3.3-56) 


Thus,  if  the  cylinder  is  considered  inelastic,  the  corresponding  partial  deriv¬ 
atives  may  be  sec  equal  to  zero.  Generally,  each  of  the  partial  derivatives 
will  be  a  function  of  pressure,  stroke,  or  both.  The  terms  in  seal  comprc.'sioa 
are  constant  if  the  seals  are  compressed  in  the  linear  range;  otherwise  they 
will  vary  with  pressure. 
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The  geometry  of  a  particular  liquid  spring  may  be  somewhat  different  from 
that  presented  here,  but  the  approach  will  remain  the  same.  The  metering 
function  An  is  retained  underneath  the  integral  signs  in  the  previous  forms,  as 
it  may  vary  with  stroke  and  hence  implicitly  with  time.  Relief  valves  for  rapid 
return  strokes  have  not  been  included  here.  They  may  be  incorporated  by  simply 
stating  that  the  press'ires  in  the  two  regions  are  identical. 


When  the  piston  strokes  within  the  cylinder  an  amount  of  kinetic  energy 
will  be  dissipated  by  friction  at  the  bearing  surfaces  between  the  piston  and 
the  cylinder.  The  frictional  force,  opposing  the  motion  of  the  piston,  is  ex¬ 
pressed  as 


(3.3-57) 


where 

Fg^bearing  forces  at  the  upper  and  lower  bearings  respectively  required 
to  balance  the  lateral  loading  on  the  piston. 

41.V*  coefficients  of  sliding  friction  at  the  upper  and  lower  bearings 
'  respectively. 

Fp0  frictional  force  at  zero  lateral  loading.  This  force  is  assumed  to 
be  known. 

The  bearing  forces  and  are  the  reactions  that  put  the  piston  in 
equilibrium  with  the  external  and  inertial  forces  in  the  piston. 

Consider  first  the  case  when  the  inertia  forces  in  the  piston  can  be 
neglected  when  determining  bearing  forces. 
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From  the  equilibrium  conditions  on  the  piston; 

ffx'G,  4 

F&LX  =  +  A. 


(3.3-58) 


-%g,  4  * 

"  4  *  ^ 


(3.3-59) 


where  4  is  the  bearing  separation  at  zero  stroke. 
The  resultant  bearing  forces  are  therefore, 


(3.  3-60) 


3-fc) 
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The  dimension  cL  is  the  length  of  the  piston,  JlD  ,  plus  the  tire  radius,  T  , 
minus  the  tire  deflection  along  the  IK.  axis,  fico&K,  or 

~^p  +  r  -  £  COS  ^4  (3.5-62) 

where  J"  an are  defined  by  Equation  3.2-55. 

In  the  general  case  when  the  inertial  loads  are  to  be  considered,  the  panel 
point  loads  on  the  piston  are  defined  from  the  equation  of  motion  of  the  piston 


OWO-C'Vi'Kp/) 


(3.3-63) 


(3.3-64) 


where  OWL  are  the  stiffness  matrices  of  the  piston  and  0*}, 

<Tp  ,  £  y  are  tue  elastic  displacements  of  the  piston.  x 

if 


■IK 


Given  the  loads  at  tne  above-indicated  panel  points  the  bearing  forces  are 
determined  from  the  equations  of  statics  as  before: 


^BLX'  +^~BUX.'  OXFp*'}  ^ 

^Blx'  (J 5  +xd')  =  &p)  {fj pX'}  ^  ^ 

where  is  the  distance  of  each  panel  point  along  JK 
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3-3-3-^  Bottoming  force,  p 

O 

The  piston  of  the  hydraulic  strut  will  be  required  to  stay  within  certain 
values  of  stroke,  consistent  with  the  construction  of  the  strut.  For  strokes 
less  than  zero  the  piston  will  contact  the  lower  bearing.  For  strokes  greater 
than the  piston  will  contact  some  mechanical  stop  at  the  upper  end  of  the 
cylinder.  To  ensure  that  the  piston  stays  within  the  required  range  of  strokes 
the  following  functions  are  introduced. 


fb  *  ^  ^ 


Bl_ 


^4  0 


(3.3-73) 


F  =  «  ^ 


^  ^  MAX 


(3. 3-74) 


where  KR..  ,  \{_  are  the  spring  rates  of  the  upper  and  lower  ends  of  the 
cylinder. 


The  values  of  the  spring  rates  are  quite  high,  and  may  arbitrarily  be  assigned 
such  that  the  displacement  past  the  mechanical  stop  is  extremely  small.  Generally, 
damping  terms  are  included  so  that  the  piston  will  cease  to  oscillate  from  the 
boUoiuiug  and  stroking  spring  forces.  This  is  a  mathematical  artifice,  and  the 
damper  rates  are  arbitrary. 


3*3-3 * 5  Crushable  materials 


The  use  of  crushable  materials  in  shock-mi oigating  devices  is  fairly 
recent.  There  is,  however,  a  wide  range  of  types  of  these  devices.  The  design 
of  such  a  device  will  depend  greatly  on  the  vehicle  mass  and  the  limits  to  be 
placed  on  its  deceleration  rate. 

Almost  all  of  the  shock-absorbing  mechanisms  using  crushable  materials  are 
designed  so  that  the  force  on  the  main  component  never  exceeds  a  certain  value 
but  remains  very  near  that  value  throughout  the  gear  stroke.  This  is  possible 
since  the  primary  characteristic  of  most  crushable  materials  is  that  of  a  constant 
load-stroke  curve.  This  generally  holds  until  the  volume  of  the  material  is 
reduced  to  one-fifth  of  its  initial  valu-.,,  at  which  time  the  crushing  character¬ 
istics  become  nonlinear. 

Some  shock-absorbing  devices  must  have  the  characteristic  that  the  rate  of 
loading  does  not  exceed  a  certain  value,  rather  than  the  load  itself.  This  is 
accomplished  by  shaping  the  leading  end  of  the  device  which  penetrates  the 
material.  Most  of  the  crushable  materials  will  exhibit  a  high  onset  force,  which 
quickly  reduces  to  the  constant  force  for  which  it  is  designed.  This  is  the 
reason  for  the  shaping  of  the  initially  crushed  surface  or  the  device  wh does 
the  crushing. 

It  has  beer,  assumed  h°re  that  the  crushable  material  is  interior  to  a 
stroking  device.  If  it  is  a  shaped  piece  of  material  simply  affixed  to  the 
underside  of  a  vehicle,  the  force  should  be  considered  as  an  exterior  appli°d 
force  on  the  body  rather  than  an  interior  force  applied  to  an  unspring  mass. 
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If  the  material  is  interior-,  the  contribution  would  be  positive; 


^  Jy  < 
fd  ^ 

and  if  it  is  exterior,  the  force  is  a  ground 
touchdown  time  t0  ,  and  is  negative, 

v  =^°  t<to 

30  l-f.  t>t. 

If  it  is  exterior,  the  contribution  to  \A  from  ground  flexibility  is  usually 
negligible.  The  drag  force  would  depend  on  the  vertical  load  and  a  coefficient 
of  friction  dependent  on  the  surface  area  of  the  material  on  the  ground. 

3. 3*3*6  Gas  compression  and  relief  valves 

This  paragraph  is  concerned  with  the  energy  absorption  device  which  makes 
use  of  gas  compression  and  release  through  an  orifice. 

A  device  of  this  type  is  usually  designed  with  a  pressure-sensitive  bleeding 
mechanism.  The  orifice  area  opened  by  this  instrument  is  designed  such  that  the 
interior  pressure  is  maintained  as  constant  as  possible.  Consider  the  figure 
below. 


0 

0 


(3-5-75) 


force  applied  to  the  vehicle  at 


The  spring  holding  the  relief  valve  is  preloaded  so  that  a  critical  pressure 
P0  is  required  to  open  the  valve.  Prior  to  the  opening  of  the  valve,  the 
pressure  inside  the  cylinder  is  found  from  the  usual  polytropic  form  in  terms 
of  the  initial  pressure  and  volume; 


PVn=  p0v0" 


(3.3-77) 


The  volume  of  the  gas  during  this  interval  is  linearly  related  to  the  stroke, so 


that 


(3-3-78) 


In  terms  of  the  area  Ay  of  the  valve  prior  to  opening  and  the  spring  preloading 
force  f  ,  the  critical  pressure  at  which  the  valve  will  open  is 


Pc=Avf 


(3.3-79) 


The  corresponding  critical  stroke 


is 


(3.3-80) 


After  the  relief  valve  is  opened,  it  is  assumed  that  the  area  of  the  opening  may 
vary  sufficiently  rapidly  that  the  cylinder  pressure  remains  at  the  critical  value. 
The  force  on  the  piston  is  then 


F- 


IAP0(I  -*/i) 
A  Pc 


-n 


^  >-^c 


(3.3-81) 


in  terms  of  the  area  A  of  the  piston. 

The  assumption  has  been  made  that  the  piston  stroking  velocity  is  suffi¬ 
ciently  low  that  the  pressure  is  uniform  throughout  the  cylinder. 

3.4  ELASTIC  BODY  FORCES 

3.4.1  General 

The  forces  and  internal  reactions  giving  rise  to  elastic  deformations  in  the 
vehicle  are  discussed  in  this  paragraph.  Perhaps  the  most  difficult  point  to 
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understand  in  the  formulation  of  the  equations  of  motion  for  an  elastic  body  is 
in  this  area.  Several  different  concepts  may  be  used  in  the  interpretation  of 
the  panel  point  equations.  These  concepts  vary  according  to  the  manner  in  which 
the  internal  reactions  in  the  vehicle  are  entered  into  the  equations  governing 
component  elastic  motions.  The  reactions  may  be  entered  entirely  as  constraints 
on  the  elastic  motions,  or  they  may  be  entered  partially  as  applied  loads.  The 
former  of  these  methods  is  perhaps  the  most  straightforward.  It  does  not,  how¬ 
ever,  lend  itself  readily  to  approximate  solutions,  as  the  transformation  to 
modal  coordinates  is  not  easily  accomplished  with  that  method.  The  latter  will 
then  be  used.  The  method  will  be  elaborated  here,  with  several  examples  which 
exhibit  the  concepts. 

Consider  a  wing  attached  in  a  cantilevered  manner  to  a  fuselage;  that  is, 
it  does  not  rotate  relative  to  the  fuselage.  The  reac cions  at  the  wing  root 
which  hold  the  wing  to  the  luselage  are  not  considered  to  be  applied  forces  or 
constraints  on  the  wing.  Their  effect  is  entered  into  the  stiffness  matrix  for 
the  wing,  so  that  the  stiffness  matrix  is  that  of  a  cantilevered  wing.  Thus,  the 
equations  defining  the  wing  modal  coordinates  in  the  paragraph  on  modal  trans¬ 
formations  will  then  produce  cantilevered  mode  shapes  for  the  wing.  Consider 
the  symbolic  fonn  for  the  wing  panel  point  equations; 

[av/](pw}+  >  (3.4-u 


where  the  stiffness  matrix  is  that  of  a  cantilevered  wing,  and  the  right-hand 
side  includes  both  the  applied  forces  and  terms  coupling  rigid  body  motion  and 
elastic  motion.  There  are  no  constraints  on  the  wing.  If  the  fuselage  elastic 
displacements  are  set  equal  to  zero,  the  wing  total  displacement  becomes  the 
wing  elastic  displacement,  and  the  left-hand  side  of  the  equation  becomes  the 
form  which  defines  the  wing  modal  coordinates.  The  equation  defining  the  wing 
panel  point  total  displacements  is 


The  matrix[T^p  ]  relates  the  displacements  of  the  wing  panel  points  -  with 
the  wing  considered  as  a  rigid  body  -  due  to  displacements  of  the  fuselage 
panel  points.  This  matrix  then  picks  out  only  the  fuselage  panel  points  located 
at  the  wing  root,  and  geometrically  defines  the  wing  displacements  due  to 
fuselage  displacements.  The  displacements  of  a  rigid  body  are  completely 
defined  by  the  displacements  of  any  three  points  in  the  body,  so  that  three 
panel  points  at  the  wing  root  are  sufficient  for  the  general  case. 


This  geometric  relationship  has  a  useful  property.  Suppose  Eq.  3.4-2  is 
substituted  into  Eq.  3.4-1,  and  the  latter  is  p’-emultiplied  by  the  transform 
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The  latter  term  on  the  right-hand  side  is  just  the  inertial  load  on  the  fuselage 
panel  points  at  the  wing  root  due  to  their  acceleration  of  the  "rigid"  wing.  The 
other  terms  on  the  right-hand  side  arc  loads  on  those  panel  points  due  to  wing 
elastic  accelerations  and  wing  externally  applied  forces.  Consider  the  symbolic 
form  for  the  fuselage  elastic  motion; 
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where  again  the  right-hand  side  includes  applied  forces  and  forces  due  to  coup¬ 
ling  of  the  rigid  body,  and  fuselage  elastic  motions.  One  of  the  applied  forces 
on  the  fuselage  elastic  motion  is  at  the  wing  root.  This  force  is  represented 
in  the  column  matrix  form  by  either  side  of  Eq.  3.4-3*  Then  Ea.  3*4-4  may  be 
written 
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where  the  starred  (*)  quantity  does  net  include  the  wing  root  applied  forces, 
and  the  right-hand  side  of  Eq.  3*4-3  has  been  used.  The  reason  that  the  left 
side  was  not  used  is  now  obvious,  for  transposition  of  the  last  tern  on  the  right 
to  the  left-hand  side  introduces  the  inertial  effects  of  the  "rigid"  wing  into 
the  total  mass  matrix  needed  in  the  definition  of  the  fuselage  modes  from 
Eq.  2.9-12; 

[A]  =  [Af]+  z  Cv ][Aw][Twf ]+ .  (3.4.6) 

In  a  like  manner,  the  applied  forces  on  the  fuselage  panel  points  from  the 
various  appended  components  may  be  evaluated.  Thus,  it  is  seen  in  this  example 
that  the  wing  root  forces  are  applied  forces  on  the  fuselage  elastic  motion,  but 
their  effect  on  the  wing  elastic  motion  is  included  in  the  definition  of  the 
wing  stiffness  matrix. 

This  approacu  will  normally  be  used  throughout  the  analysis;  the  reaction 
forces  are  applied  forces  on  the  more  central  component  and  are  included  in  the 
stiffness  matrix  for  the  more  remote  component.  Thus,  for  a  gear  affixed  to  a 
wing,  the  trunnion  loads  are  applied  forces  on  the  wing  elastic  motion,  but 
their  effect  on  the  gear  is  included  in  the  gear  stiffness  matrix.  It  is 
apparent  that  if  a  chain  of  components  affixed  to  one  another  are  all  elastic, 
this  may  become  a  lengthy  process. 

In  Section  4  it  is  pointed  out  that,  for  practical  applications,  some  of 
the  components  would  be  considered  to  be  rigid.  Suppose  that  for  a  particular 
application,  the  fuselage  may  be  considered  to  be  rigid.  There  would  *heu  be 
no  fuselage  panel  point  displacements.  Consider  an  element  which  is  pin-~:pported 
to  the  fuselage,  and  has  a  pin-supported  strut  connected  to  the  fuselage.  This 
configuration  is  pictured  below. 
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Assume  that  the  strut  is  rigid,  and  that  the  fuselage  is  rigid.  The  elastic 
motion  of  the  element  is  defined  by  its  distributed  stiffness  properties  as 
the  two  restraints.  For  this  geometry,  the  elastic  deflection  of  the  element 
must  be  zero  at  the  two  support  points.  This  effect  will  be  included  in  the 
stiffness  matrix  of  the  element;  the  stiffness  matrix  is  that  of  an  element 
simply  supported  at  two  points.  In  this  case  there  are  no  applied  loads  on 
fuselage  elastic  motion. 

Let  the  fuselage  now  be  elastic,  so  that  the  support  points  may  move. 

The  strut  is  still  assumed  to  be  rigid.  The  motion  of  the  panel  points  on 
the  fuselage  where  the  element  is  supported  obviously  determines  the  displace¬ 
ment  of  the  element  as  a  rigid  body.  Thus,  the  problem  becomes  exactly  the  same 
in  concept  <*3  the  first  case  of  a  fuselage  with  cantilevered  wing.  The  stiff¬ 
ness  matrices  of  this  element  and  of  the  wing  differ  due  to  the  change  in  support 
methods,  and  the  geometric  relations  involved  in  the  matrix  £TWfJ,  are  different 
from  those  in  this  case.  The  applied  forces  on  the  fuselage  panel  points  are 
found  in  the  identical,  manner. 

These  examples  show  the  manner  in  which  the  internal  reactions  will  be 
handled.  Cases  have  been  discussed  in  which  both  components  are  elastic,  and 
in  which  one  component  is  rigid.  If  both  components  are  assumed  to  be  rigid, 
but  are  attached  to  a  third  elastic  component,  the  same  methods  apply  in 
obtaining  applied  forces  on  the  elastic  component.  The  panel  point  elastic 
displacements  of  the  components  assumed  to  be  rigid  are  simply  set  equal  to 
zero.  The  proper  geometric  properties  are  then  used  to  enter  the  applied  forces 
on  the  elastic  component  due  xo  tne  inertia  of,  and  applied  forces  on,  the 
rigid  components. 

An  example  of  the  general  method  is  shown  for  the  simplified  geometry 
below.  The  horizontal  beam  is  simply  supported  at  two  points,  and  the  vertical 
beam  cantilevered  to  the  horizontal  beam.  Motion  is  considered  only  in  the 
plane  of  the  page,  and  each  beam  is  assumed  to  be  incompressible. 
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The  horizontal  (E)  beam  may  have  elastic  displacements  parallel  to  the  y  axis. 
The  vertical  (V)  beam  may  have  elastic  displacements  parallel  to  the  X  axis} 
and  motion  as  a  rigid  body  parallel  to  both  axes  due  to  elastic  displacements 
of  the  horizontal  beam.  The  total  displacements  on  the  vertical  beam  are 


Assume  that  the  interpolation  scheme  used  to  calculate  the  mass  matrices  is  the 
trapezoidal  scheme;  i.e.,  the  displacement  betveen  panel  points  is  a  straight 
line  ending  at  the  panel  point  displacements.  The  panel  points  labeled 
on  the  horizontal  beam  then  determine  entirely  the  motion  of  the  vertical  beam 
as  a  rigid  body.  The  displacements  of  the  panel  points  on  the  vertical  beam 
are  then  geometrically  related  to  those  of  the  panel  points  ,Jl+\  by  the 
matrix  [TVhL  which  yields  Eq.  3.4-8.  It  is  assumed  that  these  two  points  are 
spaced  a  distance  d/z  to  either  side  of  the  connection  point.  The  distances 

>'■ .  >  are  measured  from  the  connection  point  down  to  each 

panel  point  on  the  irertical -beam. 
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The  stiffness  matrix  for  the  horizontal  heam  is  that  of  a  beam  simply 
supported  at  two  points.  The  stiffness  matrix  for  the  vertical  beam  is  that 
of  a  beam  cantilevered  at  one  end.  The  motion  of  the  entire  system  as  a  rigid 
body  dees  not  occur  due  to  the  fixed  simple  supports  on  the  horizontal  beam. 

The  corresponding  rigid  body  coupling  terms  from  the  panel  point  Eqs.  2.8-5  may 
then  be  omitted.  Ho  constraints  are  formally  entered  into  the  motion,  so  that 
the  resulting  panel  point  equations  may  be  written  as 


^A? >  V *  [4  J ^ fa 
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The  quantity  00  is  discussed  in  a  moment.  For  simplicity,  assume  that  the 
only  external  applied  forces  on  the  system  occur  at  the  bottom  panel  point  on 
the  vertical  beam.  Let  this  force  have  the  components  Fx  ,  Fy  .  The  applied 
force  on  the  vertical  beam  is  then  * 
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The  applied  force  on  the  horizontal  team  is  due  only  to  the  vertical  beam 
connection.  This  is  given  by 


,Q  ,  |y  VrR,r,  6U  1  f<Px>t\ 


(3-4-12) 


The  quantity  {R^is  seen  to  be  the  difference  in  the  applied  force  and  inertial 
reaction  of  the  vertical  beam  along  its  axis.  This  quantity  cannot  be  determined 
from  stiffness  properties,  as  it  vas  assumed  that  the  beam  is  rigid  along  the 
axis.  The  stiffness  in  that  direction  is  then  infinite,  but  the  elastic 
deformation  is  zero,  and  the  product  is  indeterminant.  Thus,  f p^>is  found 
in  terms  of  the  applied  forces  and  inertial  reactions.  ^  J 


The  first  term  on  the  right-hand  side  of  Eq.  3.4-12  is  easily  evaluated; 


This  form  enters  the  effect  of  the  auplied  forces  on  the  vertical  beam  into  the 
horizontal  beam  panel  points  Jt  ,  >+r  in  the  proper  manner.  The  other  term  on 
the  right-hand  side  of  Eq.  3-4-12  essentially  subtracts  the  inertial  effects  of 
the  vertical  beam  from  the  effects  of  the  applied  forces  on  the  vertical  beam. 

In  that  term,  the  Eq.  3*4-7  will  be  substituted.  The  two  resulting  terms  are 


The  latter  will  be  transposed  to  the  left-hand  side  of  Eq.  3-';**9  and  included 
with  the  first  term  there.  This  term  enters  the  inertial  effects  of  the  vertical 
beam  as  a  rigid  body  into  the  proper  place  in  the  mass  matrix  of  the  horizontal 
b9sm=  The  remaining  term  enters  the  inertial  effects  of  the  lateral  elastic 
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The  equation  governing  motion  of  the  vertical  beam  parallel  to  the  ^--axis  could 
be  included,  but  no  additional  information  is  obtained.  This  motion  is  defined 
entirely  by  the  matrix  [Tvh]  and  the  displacements  of  the  panel  points  JC,  , 

-C+  I  .  If  this  is  kept  in  mind,  the  equation  may  be  omitted.  The  form  of 
Eq.  3-4-15  is  now  proper  for  the  transformation  to  modal  coordinates,  as  the 
mass  matrix  includes  the  rigid  body  inertia  of  the  vertical  beam.  The  Eq.  3 .4— 16 
is  in  the  proper  form  when  the  total  displacement  is  expanded  according  to 
Eq.  3.4-7. 


This  completes  the  discussion  of  the  applied  forces  on  component  motions. 
The  forces  on  the  whole  vehicle,  which  were  discussed  in  Paragraph  3*2,  must  be 
distributed  to  the  various  panel  points.  This  process  will  not  be  disc  seed 
here,  as  it  depends  only  on  the  particular  configuration.  Several  examples  in 
the  next  section  exhibit  this  process. 

Two  types  of  applied  forces  on  panel  point  elastic  motion  are  explicitly 
written  next.  One  is  the  aerodynamic  force  on  e.  wing  due  to  changes  in  the 
local  surface  angle  of  attack  caused  by  wing  elastic  deformations.  The  ether 
is  the  structural  damping  force,  which  is  physically  not  an  external  force,  but 
is  considered  to  be  as  a  convenience. 
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The  contribution  to  individual  panel  point  loads  here  will  be  considered 
only  on  airplane  wings,  and  only  that  contribution  due  to  elastic  motions.  The 
body  and  wing  ^  axes  will  bo  assumed  coincident.  The  pressure  distribution  due 
to  elastic  deformations  is 

PC(*y,t)  =  LR(x'^0  + ir  LC  fcfO  (3.4-17) 


,  /  .  1  2Cp(x,v)  ^p,'(xV,t) 


(3.4-18) 


Lc  (x' lj/,  t)  =  \/2'rtf  2  (x  Y» fc)  (5.4-19) 


=  atmospheric  density 

/IT  =  velocity  of  the  vehicle  relative  to  the  atmosphere 
<=«£  s  airplane  angle  of  attack 
Cp  =  local  pressure  coefficient 

The  virtual  work  entailed  in  a  cliange  of  displacement  5  Py  is 

JW  PC(x'^,t)  £P.  dxelij/ 

where  the  integration  is  over  both  wings. 

This  may  be  written  in  terms  of  an  interpolation  scheme  as 


(5.4-20) 


(5.4-21) 


2  2' 
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Panel  point  loads  for  either  "ins:  are  then 


KjAj  */  \ir*K?j  w  ir&j l'?j/  (3-4_22) 

It  has  been  assumed  tbrt^Q^and^Q^^are  unchanged  from  their  rigid  values. 
The  contributions  to  body  forces  and  moments  from  these  panel  point  forces  are 


(5A-23) 

Nxa  ^»{Q^  }w,+  \z 

(3.4-24) 

(3.4-25) 

These  may  be  of  significant  value  for  large  wings  with  considerable  bending. 


The  inclusion  of  contributions  from  tails  or  other  surfaces  may  be  per¬ 
formed  in  a  similar  manner  but  wiil  not  have  any  first-order  effect  on  landing 
loads. 

3.^-3  Structural  Damping 

Damping  in  a  built-up  structure  is  not  a  problem  to  be  approached  analytically. 
The  results  of  damping  are  easily  incorporated  for  systems  of  defined  normal 
modes,  however. 


In  the  test  laboratory,  the  definition  of  a  mode  of  vibration  includes  the 
requirement  that  the  shape  of  the  deformation  remain  constant  as  the  amplitude 
decreases  to  zero.  This  may  be  used  to  define  a  damping  parameter  in  the  modal 
equations.  Consider  the  equation 


[A]{P}  +[kn]{P)  ={0) 


where  [k„J  is  the  stiffness  matrix  which  defines  the  normal  modes  of  vibration. 
The  transformation  to  normalized  modal  coordinates  yields 


a  separated  set  of  equations.  The  frequency  of  the  undamped  motion  defined  by 
each  equation  iso)=:=  the  equation  governing  that  coordinate  is 

XA 
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(3.4-28) 


+  0)  2  Cl/  =  0 


Suppose  that  damping  were  included  in  some  manner  such  that 

[a]{p}+[d]{p}  +[kJ{p)  ={o> 


(3 ->*-29) 


If  transformation  to  the  (normalized)  normal  modes  is  made,  the  equation  becomes 


(3.4-30) 


The  requirement  that  the  shape  of  the  deformation  remain  unchanged  as  the  motion 
dies  out  requires  that  the  equations  for  the  modal  coordinates  be  uncoupled  even 
with  damping,  which  means  that 


mm^ 


(3.4-51) 


where  the  elements  are  the  modal  damping  coefficients.  These  elements  may 
be  assigned  values  in  terms  of  critical  damping,  but  are  not  derivable  from  basic 
properties  of  the  structure  in  most  cases.  Consider  the  equation 


(jr  +  d<^  +  (JO  =  0 


(3.4-32) 


where  OJ  is  the  frequency  of  vibration  of  the  mode  corresponding  to  as  for  the 
undamped  structure.  If  the  substitution  ® 


(3.4-33) 


is  made,  the  equation  for  f(t )  is 

'/  +  (*) 2 =  0 


The  solution  is  harmonic,  of  frequency 


(3.‘.~34) 


(3.4-35) 
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Critical  damping  is  defined,  as  that  value  for  which  h->  -  0  ,  that  is. 


«»e 
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If  the  fraction  of  critical  damping  for  the  mode  is 


or 


d  -Xtu  dc  (3.4-57) 

and  I,  the  frequency  of  the  damped  motion  is  very  nearly  that  of  the 

undamped  motion.  Thus  the  damping  term  may  be  written 


d=2uj/w  (3.4-38) 

For  many  applications,  the  fraction  of  critical  damping  of  the  basic  or 
fundamental  mode  of  vibration  may  be  taken  as  ten  percent. 

3.5  CONSTRAINTS 


The  reader  who  has  progressed  through  the  contents  of  the  report  to  this 
point  is  well  aware  of  the  important  part  the  concept  of  constraints  may  have 
in  the  formulation  of  a  complex  problem.  The  purpose  here  is  to  discuss  the 
general  concept  of  a  constraint  and  to  present  some  examples  which  may  aid  the 
reader. 

In  the  formulation  of  a  complex  problem,  it  is  often  quite  difficult  to 
retain  compactness  and  at  the  same  time  choose  a  set  of  independent  variables 
which  completely  define  the  motion.  It  is  more  often  convenient  to  choose  a 
set  of  variables  which  least  complicate  the  entire  problem  and  use  the  method 
of  constraints  to  eliminate  the  dependencies  which  have  been  included.  There 
are  also,  in  some  cases,  problems  in  which  the  applied  forces  corresponding  to 
the  independent  variables  cannot  be  defined,  but  may  be  found  as  constraint 
forces  using  a  set  of  dependent  variables.  These  concepts  are  discussed  in 
Ref.  1.  In  that  reference,  it  is  shown  that  the  constraints  must  be  of 

a  particular  form  in  order  that  they  may  be  handled  by  the  methods  presented 
here. 

On  the  right-hand  side  of  Eq.  2.8-5  is  a  set  of  terms 


(5.5-1) 


The  functions^,  are  the  constraint  relations;  they  are  algebraic  relations 
between  some  of  the  panel  point  displacements.  These  algebraic  relations  must 
be  of  a  particular  form.  They  must  relate  only  the  variables  already  defint  *. 
in  the  set  of  equations  of  motion;  they  cannot  introduce  new  variables.  They 
must  be  written  as 


> . »  >0  "  0  (5*5-2) 

where  jk  ,Jl>  •  *  *  •  ,  Al  indicate  the  particular  panel  point  displacements  that  are 
not  independent.  Note  that  this  form  does  not  include  inequalities. 

The  elements  (TV  in  Eq.  3.5-1  are  called  Lagrange's  undetermined  multipliers. 
They  are  the  factors'which  convert  each  partial  derivative  into  a  constraint 
force  consistent  with  the  constraint  relations.  These  elements  are  determined 
only  by  solution  of  the  equations  of  motion,  and  they  are  generally  complicated 
functions  of  the  panel  point  motion  throughout  the  system.  They  may,  however, 
be  eliminated  before  solving  the  equations  for  many  cases  of  interest  in  this 
report. 

A  single  problem  to  illustrate  the  constraint  concept  is  pictured  in 
Fig.  9. 


5  fe _ 7 

0  o  o 

E 


O  I 


E 


Fig.  9.  A  Constraint  Example 

In  the  left  part,  two  beams  are  simply  supported  at  two  points.  Each  beam  is 
assigned  four  panel  points,  and  each  beam  is  then  allowed  four  dog-ees.  of  lree- 
dom  for  elastic  notion.  The  equations  of  motion  are 
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(3-5-3) 


CA]a{P}a+CK]a{P}a={Q}a 

LK  '¥-%^  m/r  (?-5J0 


Sr?."®  freedom  altogether .  Suppose  now  a  rigid,  weightless 

bar  is  connected  between  tne  beams  at  panel  points  three  and  seven.  These  two 

^“strained  to  move  together  or  have  the  same  displacements. 
This  constraint  relation  is  written  as 


F,  =P3  -P7  =0 


(3-5-6) 


which  is  of  the  general  fom  of  Eq.  3.5.2.  The  partial  derivatives  of  this 
constraint  relation  with  respect  to  the  panel  point  displacements  are 


The  Eqs.  3.5-3,  b  are  then  no  longer  valid;  the  constraint  term  must  be  added. 
The  equations  of  motion  vhich  must  be  solved  simultaneously  for  the  coupled 
elastic  motion  are  then 


P3  -P7  =0  (3.5-10) 


There  are  now  the  eight  panel  point  displacements  and  one  Lagrangian  undetermined 
multiplier  or  nine  unknown  variables  and  nine  equations  of  motion.  On  the  right- 
hand  side  are  the  eight  panel  point  applied  forces  and  the  constraint  forces  at 
panel  points  three  and  seven  which  cause  those  panel  points  to  move  together . 

The  reader  whose  experience  with  these  concepts  is  limited  should  now 
examine  in  detail  the  Sections  2.8,  9  in  order  to  obtain  some  working  knowledge 
on  the  subject.  The  section  or.  component  rigid  body  motion  alorg  a  line  follows 
directly  from  the  simple  example,  except  that  the  number  of  constraints  involved 
is  large.  Fortunately,  in  that  and  the  following  cases  of  use  of  the  constraint 
methods,  the  undetermined  multipliers  may  be  eliminated.  The  constraint  forces 
for  each  case  are  simply  the  forces  necessary  to  bold  the  component  together  to 
move  as  rigid  body  under  the  applied  and  inertial  forces. 


121 


APPENDIX  A 


DERIVATION  OF  TOE  KINETIC  ENERGY 


The  details  of  the  derivation  of  the  kinetic  energy  as  presented  in 
Section  2  are  written  here.  Several  vector  and  matrix  identities  are  established 
for  later  use. 


A  vector  in  three-space  rnay  be  written  in  matrix  form  as 


(A 


-1) 


(A-2) 
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In  the  same  manner,  the  vector  product  is  written  as 


Ax  IB  = 


and  the  triple  scalar  product  as 


(A-5) 


A-IB*ID= 


A  coordinate  transformation  may  he  written  in  this  manner: 


(A-U) 


TOO 


Then 


l/J- 


m 

~l-V  l-f  A  •  IK  1 

o'>= 

j  *  a'  j.jt  j.ik 

llK'J 

.IK-T  IK*  J*  !K •  IK  .. 

(A-5) 


The  elements  of  the  latter  matrix  may  be  recognized  as  the  direction  cosines  of 
the  unit  vectors  of  one  coordinate  system  in  the  other  system.  One  may  then  write 


With  these  identities,  the  kinetic  energy  from  Eq.  2.6-9  may  be  put  in  matrix  form. 
Beginning  with 


(A-7) 


(A-8) 


In  the  second  term,  the  identity 


(JTL*  H)L  =  *  (lL*.&) 


may  be  written  using  Eqs.  A-2,4  as 
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(A-10) 


The  coordinate  system  of  the  body  is  defined  as  the  principal  axis  of  inertia 
system,  so  that  the  off-diagonal  elements  integrate  to  zero,  and  the  second  term 
becomes 


{]'v(nxiL)VdV=l 


(A-ll) 


In  the  remaining  terms  of  the  kinetic  energy,  the  displacement  vector  of  the 
elemental  volume  occurs.  If  these  terms  ere  to  be  evaluated,  some  scheme  must 
be  devised  by  which  the  continuous  displacements  may  be  represented  by  the 
displacements  of  the  discrete  set  of  panel  points.  This  is  accomplished  by  a 
numerical  procedure  known  as  an  interpolation  scheme.  Interpolation  schemes  re¬ 
late  the  displacements  in  the  regions  between  the  panel  points  to  the  panel  point 
displacements.  Consider  the  next  term  in  the  kinetic  energy.  It  may  be  written 
as 


12? 


since 


and 

K  Mi  =  M 


(A-13) 


(A-14) 


In  terms  of  some  interpolation  scheme,  each  component  of  the  continuous  dis¬ 
placement  (or  velocity)  is  related  to  the  values  at  the  panel  points.  Thus, 


(A-15) 


The  resulting  form  is  simply  a  quadratic  expansion  in  the  panel  point  velocities. 
The  values  of  the  elements  in  the  mass  matrix  [Ax' x’j  will  depend  on  the  inter¬ 
polation  scheme  used,  which  in  turn  dictates  the  accuracy  of  the  representation 
of  the  continuous  system.  The  third  term  in  the  expansion  of  the  kinetic  energy 
is  then 


1 

l 


1 


(A-16) 


In  the  remaining  terms,  the  s\xbscr.ipt  i  is  omitted  until  the  final  form  for  the 
kinetic  energy  is  vritten,  it  being  understood  that  all  elements  must  be  written 
in  the  proper  coordinate  system. 


In  the  fourth  term. 


(SL*  !P)1; 


■Pl) 
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-p  p 

*  t 

-p  p 

*  7 


-P  P 

'x  ri 


-p?pJ , 


(A-17) 


From  the  properties  of  a  similarity  transformation  (see  Ref.  l),  this  may  be 
written 
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.  fM 

//*..  rr\V‘ _  ) r\  V  fm") 

\jfL-in  -  Uy  / 1*  j 

luj 


-Px'V 

H'  {P^P/}  "p<Va  ! 
"V *V 
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■Px'P| 
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M'ysf 
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Y 


\H-XO) 


This  may  be  verified  by  direct  expansion.  Finally, 


£[*(&*  IP)Vdv  = 
ISL 


LA. 


({f^[»,vKP^<Pj'>'f»jyKPj'» 

-{rJK-AirJ 

-{tj'bJfKfy 

-{Pj>'(ijrKPy> 

-  <Py>'  Dw  ]<V> 

-<P/>'[Ay>.]  <P?-> 

Thus,  since  the  quadratic  form  expressed  in  Eq.  A-l8  has  coupling  terms  between 
the  components  of  the  displacement,  the  form  in  Eq.  A- 19  will  have  mass  matrices 
with  the  subscripts  indicating  the  coupling.  These  mass  matrices  will  differ  from 
those  on  the  main  diagonal  in  general,  due  to  the  variation  in  mass  density,  panel 
point  spacings,  and  possible  variations  in  the  interpolation  schemes  used  in  the 
various  directions  in  the  component  coordinate  system. 


W 

(A-19) 


The  fifth  term  may  be  written  in  either  of  two  forms: 


[R  •  Jl  *  IP 


-  -ft.*  ik  *  ip 


(A -20) 


Defining 


[flj=  -ii? 
■■n"t 


f  °  /ir^  -  /ir^ 

[/v'l  =  -ATf  o  a rx 

./IT  -ATX  0  . 

then 


(A-21) 


(A-22) 


**•■**- ^  ) [si] { p„  y =  (o  } m  { p, 
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(A-23) 


(A-24) 


fA*v  fe  1  rM  ]  f{p*>] 

<n>(/ir)M  {i>  (M  <{p„->> 

1  IaJ  i  (i>J  L  u?.i j  [{pj)| 

The  new  type  of  mass  matrix  introduced  here  arises  from  forms  of  the  type 


f  P%'/»dV  ={l>'lAx-]<P>.->  (*-85) 

JVl 

The  same  interpolation  scheme  as  in  Eq.  A- 15  would  be  used  to  obtain  \fcis  result, 
but  the  elements  of  the  single  subscript  type  of  mass  matrix  will  In  general  be 
different  from  those  of  the  double  subscript  type.  The  following  shorthand 
notation  has  been  adapted: 
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(A-26) 


The  sizes  of  the  indicated  column  vectors  are  dictated  by  the  number  of  panel 
points  whose  displacements  are  indicated  in  the  adjacent  matrices. 


The  sixth  term  is 
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In  the  seventh  term, 
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(A -28) 


(A-29) 


The  mass  matrices  containing  components  of  the  position  vector  IL-  inside  ..he 
matrix  brackets  to  the  left  of  the  mass  matrix  symbol  arise  from  integrations 
of  the  form 


fvl  =  {'  y  [>a-x'3  (Px'J 


(A- 


7* 
J  * 


\ 
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In  the  limit  ns  the  number  of  panel  points  becomes  large,  and  the  interpolation 
scheme  becomes  increasingly  accurate. 


y.  •*»  r  .  .  1  CA  \ 

v;  u'vjw 


c  \ 1  r. 

IAX 


/•  \ 

'J  \JV/ 


(A-31) 


where  the  column  matrix  ,<v>  indicates  the  positions  of  the  panel  points  in  the 
component  coordinate  system,  along  the'}'  -axis.  Although  the  original  forms 
will  be  retained  in  defining  the  kinetic  energy,  the  equations  of  motion  will  be 
derived  using  the  approximation  that  Eq.  A-31  is  true  without  reservations. 


Likewise  the  eighth  term  is 

J7L'  f  .  pxliVdv  - 

Jvt 

•  (o}'  ^>'  M  f<U 

{P,>'M  <0>'  0,) 
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In  the  ninth  term, 
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so  that 
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JUL*s4-U-*iP)/,dv  - 


M 


-<ihx'.vKp,-> 
-  Cl  >'  [x-A  -]  CP  > 


-{|>'U'AX-J<PX'> 


(A-5U) 


-<>}'  [vAx-K?x'> 

-{Ob’A«']{Pu) 


U  J  n»'  J  V  <4  i 


Each  of  these  terms  is  written  for  i  -th  component  and  the  sum  of  these  is  the 
kinetic  energy.  These  terms  are  collected  next.  The  summation  on  i  goes  from 
i  *  1  to  i  a  K.  Instead  of  distributing  subscripts  throughout,  the  entire  set 
of  terms  to  be  summed  over  is  indexed.  The  index  will  fall  on  all  elements  except 
the  rigid  body  velocities. 


131 


The  desired,  form  for  the  kinetic  energy  is  written  in  Eq.  A-35.  It  will 
be  substituted  into  the  modified  set  of  Lagrange's  equations,  which  are  derived 
in  Appendix  B. 
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(A-55) 
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End  Suramation 


APPENDIX  B 


TRANSFORMATION  TO  THE  MODIFIED  FORM  OF  LAGRANGE'S  EQUATIONS 


The  details  of  the  derivation  of  the  modified  form  of  Lagrange's  equations 
as  presented  in  Section  2  are  written  here. 


Lagrange's  equations  are  first  written  as 
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(B-2) 


M 


where  (X,  Y,  2),  (fad ,  ),  )  and  their  time  derivatives  com¬ 

pletely  specify  the  position  and  velocity  or  every  particle  in  the  system  relative 
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to  the  inertial  frame  of  reference.  The  Lagrangian,  L  =  T  -  U,  is  expressed  in 
terms  of  the  body  linear  ana  angular  velocities,  (4r>.  ,4C,flQ  and  )  > 

and  the  panel  point  displacements  and  velocities.  These*’ are  not  a  suitable  riot 
of  variables  for  use  in  the  ordinary  form  of  Lagrange's  equations,  as  they  do  not 
specify  the  motion  relative  to  an  inertial  frame  of  reference.  The  tas.  P-1, 

2,3  are  valid,  however,  and  may  be  transformed  so  that  all  operations  in  the 
equations  act  on  the  variables  used  to  define  the  Lagrangian. 

The  transformation  matrices  [p]  and  [R.]  and  the  matrices  and 

may  be  recalled  from  Eqs.  2.7-^>  5>  11,  12.  Their  use  as  a  brief  notation  is 
valuable  here. 

The  first  operation  indicated  in  Eq.  B-l  may  be  transformed  by  use  ot  Vie 
chain  rule  for  differentiation  in  calculus  and  the  definition  of  the  matrix  fp]  ; 
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Since  the  kinetic  energy  is  independent  of  the  inertial  displacements  (X,  Y,  Z),  -then 
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If  no  kinematic  constraints  are  imposed  on  the  "rigid-body"  motion  of  the  vehicle, 
then  (X,  Y,  Z)  are  independent  coordinates,  and 
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The  time  derivative  of  Eq.  B-5  is 
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(B-8) 
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Since  the  potential  energy  is  not  dependent  on  the  body  coordinates, 
B-l  is  written  using  Eqs.  B-6,  7>  8  as 
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The  generalized  forces  are  transformed  by 


into  the  body  coordinate  system.  Premultiplication  by  [r]  of  Eq.  B-! 
with  the  identity  O']  OT  and 


then  Eq. 


(B-9) 


(B-10) 


together 


(B-ll) 
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yields  one  of  the  desired  transformed  equations 


pT  ] 

r  n 

Qv 

✓V. 

{'3  * 

3T 

d/tfy. 

>=< 

s 

3T 

ld/%  J 

a. 

The  transformation  of  Eq.  B-2  is  a  bit  more  tedious, 
differentiation  and  the  definition  of  w , 
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Using  the  chain  rule  oi 
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also,  the  chain  rule  yields 
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which  may  be  shown  to  produce 


r,.-\  (B-I5) 

using  the  definition  of  [/v  J. 


Differentiation  of  Eq.  B-13  with  respect  xo  time  yields 
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If  no  kinematic  constraints  are  imposed  on  the  "rigid-body"  motion  of  the  vehicle, 
then  are  independent  coordinates  and 
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Substitution  of  Eas.  B-15,  16,  17  into  Eq.  B-2  yields 


(B-l8) 


since  the  terms  explicitly  written  in  Eq.  3»15  and  Eq.  B-l6,  when  subtracted, 
yieid  -  ORJM 

The  components  (N^.Ns  )  of  Eq.  B-l8  are  the  generalized  forces  associated 
with  the  coordinates  {ifr,  0,  <f>  ) .  They  are,  physically,  moments  about  the  line 
of  nodes,  and  are  not  mutually  orthogonal  components  of  a  vector.  They  trans¬ 
form  into  the  body  coordinate  system  by 
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Premultiplication  of  Eq.  B-l8  by  [rJ 


(B-19) 

then  yields  the  modified  form 


l4o 


The  desired  transformations  are  thus  obtained, 
modified  form  are 


Lagrange's  equations  in 
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APPENDIX  C 


DERIVATION  OF  THE  EQUATIONS  OF  MOTION 

Some  of  the  details  in  the  derivation  of  the  equations  of  motion  will  he 
presented  here. 

The  forms  for  the  kinetic  and  potential  energy  of  the  system  will  be  sub¬ 
stituted  into  the  modified  set  of  Lagrange's  equations.  A  number  of  useful 
identities  will  first  be  established  in  order  that  the  operations  indicated  in 
The  modified  equations  may  be  performed  in  matrix  form  rather  than  individually. 

If  an  arbitrary  function,  Q,  ,  has  the  form 


where  the  indicated  variables  are  also  arbitrary,  then  it  is  easily  shewn  that 
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also,  in  ohe  same  sense,  if 
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With  the  use  of  these  identities  and  the  alternrte  form  of  Eq.  A-24,  the  follow¬ 
ing  operation  on  the  kinetic  energy  is  obvious: 
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The  equations  of  motion  governing  the  body  linear  velocities  follow  immediately 
from  Eq.  B-21,  and  are 


The  derivation  of  the  equations  governing  the  angular  velocities  of  the  body 
are  more  .Lengthy  due  to  the  number  of  terms  in  the  kinetic  energy  involving 
the  angular  velocities.  Using  the  previous  identities,  the  remaining  operation 
in  Eq.  B-22  is  written  in  Eq.  C-8. 
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This  form  cogether  with  that  in  Eq,  0-6  is  substituted  into  The  second  of  the 
modified  equations.  A  time- dependent  moment  of  inertia  is  defined  by  collect¬ 
ing  together  all  the  resulting  terms  which  multiply  the  angular  accelerations. 
The  moment  of  inertia  matrix  is  given  by  Eq.  C-9.  In  terms  of  this  matrix,  the 
equations  governing  the  angular  velocities  of  the  body  are  given  by  Eq.  C-10. 
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The  panel  point  equations  are  more  easily  derived  if  use  is  made  of  the 
identities 


The  rigid  body  velocities  are  then  written  in  each  component  coordinate  system. 
It  will  be  understood  that  primes  on  their  subscripts  moan  that  they  must  be 
written  in  the  i-th  system.  This  avoids  considerable  complication  in  the 
equations,  for  many  of  the  terms  cannot  be  broken  down  into  simple  product  form. 

The  same  identities  for  column  operators  may  be  used  to  advantage  where 
panel  point  displacements  appear  in  ordered  column  form.  Where  they  appear  in 
complicated  matrix  form,  the  matrix  may  be  expanded  and  the  operations  performed 
separately.  Noting  that 
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(note  also  the  primes  on  the  it  subscripts) ,  the  panel  point  equations  may  be 
written  in  the  form  of  Eq.  C-13. 
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APPENDIX  D 


INTERPOLATION  SCHEMES 


The  method  of  analysis  ir  this  report  requires  that  the  elastic  displace¬ 
ments  of  a  continuous  system  be  approximated  by  displacements  at  a  finite  number 
of  discrete  points.  This  numerical  approach  is  general  in  that  it  may  be  applied 
to  bodies  of  any  shape.  In  order  that  the  forms  remain  simple,  the  illustration 
here  will  be  confined  to  one-dimensional  motion  of  a  thin  beam.  The  purpose 
of  the  scheme  is  to  make  use  of  an  interpolation  formula  to  relate  the  continuous 
displacement  at  all  points  to  the  discrete  displacement  at  neighboring  points. 


The  region  between  panel  points^,  andi-f-l  is  referred  xo  as  the  -th  bay. 
Local  bay  coordinates  are  defined  by  * 

%  --  xi  ~  x-  xi*>  o-n 

Hence  Z  varies  between  zero  and  one. 

Various  interpolation  schemes  are  as  follows. 

Lumped  Mass 

This  method  is  not  properly  an  interpolation  scheme  but  is  included  since 
it  is  the  most  easily  used  method.  A  row  of  panel  points  is  laid  alon^  the  center- 
line  of  the  beam.  In  the  region  between  each  pair  of  panel  points,  the  mass  and 
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center  of  mass  are  calculated.  This  mass  is  then  beamed  out  to  the  pair  of  panel 
points.  When  this  process  is  completed,  the  mass  of  the  beam  is  represented  by 
a  series  of  lumped  masses  at  the  panel  points.  For  this  case,  the  mass  matrix 
is  diagonal;  the  elements  along  the  diagonal  are  the  lumped  masses. 

Trapezoidal  Buie 

This  rule  is  called  a  two-point  rule  because  the  displacements  in  a  bay  are 
defined  in  terms  of  the  displacements  at  the  two  neighboring  panel  points: 


V'  i 


In  terms  of  local  bay  coordinates,  this  may  be  written 


This  rule  is  useful  in  calculating  mass  matrices  which  do  not  have  to  be  highly 
accurate  but  are  wished  to  be  more  representative  of  the  system  than  would  be 
a  lumped  mass  approximation.  This  rule  is  also  somewhat  limited;  it  cannot  be 
used  to  calculate  the  beam  curvature  as  there  is  no  second  derivative  of  P  . 

The  kinetic  energy  of  a  beam  may  be  written 


I  rL  °  z 

T=— J  P(x)  /7Cx)dx 


where  /°(yc)  is  the  mass  per  unit  length  along  the  beam.  This  may  be  written  in  the 
form 


o  jy  o  ^  ^ 

4 

wher°//  is  the  length  of  the  ^  -th  bay  between  panel  points  ^  ,  ^.-tl  .  The 
displacement  may  be  written  in  the  form 


P(e)  = 


MM  O 


uiPi-, 


which  yields 


I  '1 
0  ! 


The  kinetic  energy  in  the  -j  -th  bay  is  then 


(D-7) 


(D-8) 


(0-9) 


The  displacements  ,  Rj+'(  may  be  picked  out  of  all  the  displacements  by  an 
operation  of  the  form 


(D-10) 


where  the  matrix  is  composed  of  zeroes  except  for  unity  in  two  elements.  The 
kinetic  energy  is  then 


(D-ll) 


or 


T=  ~  {P}'[A]{P} 


(D-12) 


The  mass  matrix  obtained  by  this  process  will  have  non-zero  elements  on  the  diagonal 
and  in  the  first  positions  off  the  diagonal. 
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Many  interpolation  schemes  have  been  devised  which  make  use  of  an  increased 
number  of  neighboring  panel  point  displacements.  Cubic  rules  make  use  of  three 
of  these,  and  four-point  rules  use  two  on  each  side  of  the  local  bay.  Of  the 
four-point  rules  in  existence,  one  will  be  chosen  for  description.  This  is  referred 
to  as  the  diparabolic  formula.  It.  is  felt  to  be  superior  to  the  ethers  since 
continuity  of  slope  is  preserved  between  bays. 

Diparabolic  Scheme 

If  -f-  (?)  is  the 
and*-  f  and f; , ,  (?) 

Xj*,1  ’ 

o<g<i  (D.15) 

is  the  formula  for  interpolation  in  the  ^  -th  bay.  Thus  the  formula  is  a  weighted 
average  of  the  two  parabolas. 


parabola  drawn  through  the  set  of  points  at  Yj-t  >  "XV 
is  the  parabola  drawn  through  the  set  of  points  at  Y- 

i 


then 


This  may  be  described  in  terms  of  the  displacements  at  the  panel  points.  The 
parabolas 

(D-14) 

,  .  01  CZ1  (Si  - 

(g)=a-pt  £  +  a/H  ?" 

(D-15) 

are  defined  by 

W=Pj- 

(D-16) 

fjrt  ~  =  ’  ^j+>  Pj+z 

(D-17; 
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The  constants,  ,  may  be  evaluated  in  terms  of  the  panel  point  deflections  by 
combining  the  above  forms.  They  may  be  written 


15* 


Formulas  for  the  parabolas  are  then 
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The  diparabolic  formula  is  then 

r*. 


P05)*{l  5  §*}((l-5)W  Jpj 


iS*y 
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which  is  of  the  form 
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Numerically,  for  eq.ua!  intervals. 
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The  kinetic  energy  of  a  beam  is  written 


L  •  2 

P'x)  V  C-JOdx 

0 


(D-32) 


where  /o(x)  is  the  mass  per  unit  length. 
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In  the  ^  ~ih  bay. 
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The  kinetic  energy  of  the  j  -th  bay  is  then 
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If  the  displacements  in  the 
by  an  operation  of  the  form 


< 


Pi 

PjH 


CP} 


1+zj 


-th  bay  are  picked  from  all  the  displacements 


(d-3T) 


vhsrG 


[&]; 


is  a  matrix  whose  elements  are  one  or 


zero, 


then  one  may  write 


T  =  i{P>'[A]{P} 
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